INTELLIGENT EARLY WARNING METHOD OF MEMBRANE FOULING

Abstract

To solve problems of frequent occurrence and great harm of membrane fouling during MBR wastewater treatment process, the invention proposes a membrane fouling intelligent early warning method to realize online and accurate early warning of membrane fouling. This early warning method achieves accurate prediction of water permeability by constructing soft-computing model based on recurrent fuzzy neural network. The intelligent early warning of membrane fouling is achieved by the comprehensive evaluation method, which solves the problem that membrane fouling is difficult to be early warning in the MBR wastewater treatment process, improves the pretreatment ability of membrane fouling, reduces the damage caused by membrane fouling, ensures the safe operation of MBR wastewater treatment process, and promotes efficient and stable operation of MBR wastewater treatment plant.

Claims

1. Membrane bioreactor-MBR membrane fouling intelligent early warning method, including data acquisition of the running process, data pretreatment of the running process, intelligent prediction of membrane fouling, and intelligent early warning of membrane fouling, comprising the following steps: (1) data acquisition of the running process: data are collected by the acquisition instrument installed on the process site, including: water flow, water pressure, chemical oxygen demand, pH, biological oxygen demand, total phosphorus, oxidation-reduction potential in anaerobic zone, oxidation-reduction potential ORP in anoxic zone, dissolved oxygen in aerobic zone, nitrate in aerobic zone, aeration; the acquired data is transmitted to the Programmable Logic Controller through Modbus communication protocol, and Programmable Logic Controller transmits the process data to the host computer through RS232 communication protocol; the data in the host computer is transmitted to the data processing server through the local area network; the process data is displayed to the management personnel in wastewater treatment plant through the Web server by the way of the Browser/Server, and the results of water permeability prediction and the membrane fouling early warning are displayed by the way of Client/Server; (2) data pretreatment of the running process: taking the membrane pool operation data as the research object, the characteristic analysis model is established by partial least squares method to extract five principal component variables, which are water flow, water pressure, aeration, ORP in anaerobic zone and nitrate in aerobic zone; these five principal component variables as input variables of the membrane fouling intelligent prediction model, and water permeability as the output variable of the membrane fouling intelligent prediction model; (3) intelligent prediction of membrane fouling: establish soft-computing model to achieve water permeability prediction based on recurrent fuzzy neural network, the structure of recurrent fuzzy neural network contains four layers: input layer, membership function layer, normalized layer and output layer, the network is 5MM1, M is an integer and 2<M<30; connecting weights between input layer and membership function layer are assigned 1, the output of recurrent fuzzy neural network is y(t); the prediction method of water permeability based on recurrent fuzzy neural network is: $\begin{array}{cc}y\left(t\right)=f\left(x\left(t\right)\right)=\underset{j=1}{\overset{M}{.Math.}}.Math.w_j\left(t\right).Math.\underset{i=1}{\overset{5}{.Math.}}.Math.\exp \left[-\frac{{\left[{}_{\mathrm{ij}}\left(t\right).Math.x_i\left(t\right)+{}_{\mathrm{ij}}\left(t\right).Math.O_{\mathrm{ij}}^2\left(t-1\right)-m_{\mathrm{ij}}\left(t\right)\right]}^2}{{\left({}_{\mathrm{ij}}\left(t\right)\right)}^2}\right],& \left(1\right)\end{array}$ where x(t)=[x.sub.1(t), x.sub.2(t), x.sub.3(t), x.sub.4(t), x.sub.5(t)] the output vector at time t, x.sub.1(t) is the value of water flow, x.sub.2(t) is the value of water pressure, x.sub.3(t) is the value of aeration, x.sub.4(t) is the value of ORP in anoxic zone, and x.sub.5(t) is the value of nitrate in aerobic zone, f is the corresponding relation between y(t) and x(t), w.sub.j(t) is the jth weight between normalized layer and output layer, .sub.ij(t)=1 is the weight between the ith neuron in input layer and the jth neuron in membership function layer, m.sub.ij(t) is the ith element of the center values of the jth neuron in the membership function layer and .sub.ij(t) is the ith element of width values of the jth neuron in the membership function layer, .sub.ij(t) is the feedback weight in the membership function layer, O.sup.2.sub.ij(t1) is the feedback value of the membership function layer, where
O.sub.ij.sup.2(t1)=exp{[.sub.ij(t1)x.sub.i(t1)+.sub.ij(t1)O.sub.ij.sup.2(t2)m.sub.ij(t1)].sup.2/(.sub.ij(t1)).sup.2}, (2) where .sub.ij(t1)=1 is the weight between the ith neuron in input layer and the jth neuron in membership function layer, m.sub.ij(t1) is the ith element of the center values of the jth neuron in the membership function layer and .sub.ij(t1) is the ith element of width values of the jth neuron in the membership function layer, .sub.ij(t1) is the feedback weight in the membership function layer at time, O.sup.2.sub.ij(t2) is the feedback value of the membership function layer; the error of recurrent fuzzy neural network is: $\begin{array}{cc}E\left(t\right)=\frac{1}{N}.Math.\underset{i=1}{\overset{N}{.Math.}}.Math.{\left(y_d\left(t\right)-y\left(t\right)\right)}^2,& \left(3\right)\end{array}$ where N is the number of samples, y.sub.d(t) is the output of recurrent fuzzy neural network at time t, y(t) is the actual output at time t, the model is trained as: 1) give a recurrent fuzzy neural network, the initial number of neurons in membership function layer and normalized layer are M, M>2 is a positive integer; the input of recurrent fuzzy neural network is x(1), x(2), . . . , x(t), . . . , x(N), correspondingly, the output is y.sub.d(1), y.sub.d(2), . . . , y.sub.d(t), . . . , y.sub.d(N), the number of training samples is N, expected error value is set to E.sub.d, E.sub.d(0, 0.01), the assignment interval of each variable in the initial center values m.sub.j(1) is [2, 2], m.sub.j(1)=(m.sub.1j(1), m.sub.2j(1), . . . , m.sub.ij(1)), m.sub.ij(1) is the initial value of the ith element of the center values of the jth neuron in the membership function layer, the assignment interval of each variable in the initial width values .sub.j(1) is [0,1], (1)=(.sub.1j(1), .sub.2j(1), . . . , .sub.ij(1)), .sub.ij(1) is the initial value of the ith element of width values of the jth neuron in the membership function layer, .sub.ij(t1) is the feedback weight in the membership function layer at time t1, the assignment interval of the initial feedback connection weight .sub.ij(1) is [0, 1], j=1, 2, . . . , M; the assignment interval of each variable in the initial weights w(1) is [1, 1], w(1)=(w.sub.1(1), w.sub.2(1), . . . , w.sub.j(1)), w.sub.j(1) is the connection weight between the jth neuron of normalized layer and the output layer at the initial time; 2) set the learning step s=1; 3) calculate the output y(t) of recurrent fuzzy neural network according to Eq. (1), exploiting gradient descent algorithm: $\begin{array}{cc}{m}_{\mathrm{ij}}\left(t+1\right)=m_{\mathrm{ij}}\left(t\right)-{}_m.Math.\frac{1}{{}_{\mathrm{ij}}^2\left(t\right)}.Math.\left(y_d\left(t\right)-y\left(t\right)\right).Math.w_j\left(t\right).Math.O_{\mathrm{ij}}\left(t\right)\left[O_{\mathrm{ij}}\left(t\right)-m_{\mathrm{ij}}\left(t\right)\right],& \left(4\right)\\ {}_{\mathrm{ij}}\left(t+1\right)={}_{\mathrm{ij}}\left(t\right)-{}_{}.Math.\frac{1}{{}_{\mathrm{ij}}^3\left(t\right)}.Math.\left(y_d\left(t\right)-y\left(t\right)\right).Math.w_j\left(t\right).Math.O_{\mathrm{ij}}\left(t\right).Math.{.Math.O_{\mathrm{ij}}\left(t\right)-m_{\mathrm{ij}}\left(t\right).Math.}^2,& \left(5\right)\\ .Math.{}_{\mathrm{ij}}\left(t+1\right)={}_{\mathrm{ij}}\left(t\right)-{}_{}\left(y_d\left(t\right)-y\left(t\right)\right).Math.w_j\left(t\right).Math.O_{\mathrm{ij}}\left(t\right).Math.y\left(t-1\right),& \left(6\right)\\ .Math.w_j\left(t+1\right)=w_j\left(t\right)-{}_w\left(y_d\left(t\right)-y\left(t\right)\right).Math.O_{\mathrm{ij}}\left(t\right),& \left(7\right)\end{array}$ where .sub.m is the learning rate of the center m.sub.ij, .sub.m(0, 0.01], .sub. is the learning rate of the width .sub.j, .sub.(0, 0.01], .sub. is the learning rate of the feedback connection weight .sub.ij, .sub.(0, 0.02], .sub.w is the learning rate of the connection weight w.sub.j, .sub.w (0, 0.01], m.sub.ij(t+1) is the ith element of the center values of the jth neuron in the membership function layer at time t+1 and .sub.ij(t+1) is the ith element of width values of the jth neuron in the membership function layer at time t+1, .sub.ij(t+1) is the feedback weight in the membership function layer at time t+1, w.sub.j(t+1) is the connection weight between the jth neuron of normalized layer and the output layer at time t+1; 4) calculate the performance of recurrent fuzzy neural network according to Eq. (3), if E(t)E.sub.d, go to step 3); if E(t)<E.sub.d, stop the training process; (4) intelligent early warning of membrane fouling: establish a comprehensive evaluation model of membrane fouling level based on the predicted values of water permeability combining with other process variables, which is specifically as follows: 1) determine the warning evaluation index of membrane fouling, set U(t)=={u.sub.1(t), u.sub.2(t), u.sub.3(t), y(t)} as the evaluation indicator vector, u.sub.1(t), u.sub.2(t), and u.sub.3(t) represent the values of water flow, water pressure and aeration, y(t) is the predicted values of permeability; 2) establish the membership functions and fuzzy comprehensive assessment matrix, membership functions reflect the relationships between the quality measurements and the defined risk levels, the membership of the evaluation factor is obtained by bringing the measured values into the membership function, the membership degree matrix R(t) can be represented as $\begin{array}{cc}R\left(t\right)={\left(r_{\mathrm{ij}}\left(t\right)\right)}_{44}=\left(\begin{array}{cccc}{r}_{11}\left(t\right)& r_{12}\left(t\right)& r_{13}\left(t\right)& r_{14}\left(t\right)\\ r_{21}\left(t\right)& r_{22}\left(t\right)& r_{23}\left(t\right)& r_{24}\left(t\right)\\ r_{31}\left(t\right)& r_{32}\left(t\right)& r_{33}\left(t\right)& r_{34}\left(t\right)\\ r_{41}\left(t\right)& r_{42}\left(t\right)& r_{43}\left(t\right)& r_{44}\left(t\right)\end{array}\right),& \left(8\right)\end{array}$ where r.sub.ij(t) (i=1, 2, . . . , 4; j=1, 2, . . . , 4) indicates the membership degree of the ith index and the corresponding jth risk rank; the membership degrees of water flow in different risk rank are $\begin{array}{cc}{r}_{11}\left(t\right)=\{\begin{array}{c}1,u_1\left(t\right)200\\ \left(300-u_1\left(t\right)\right)/100,200<u_1\left(t\right)300\\ 0,u_1\left(t\right)>300\end{array},& \left(9\right)\\ r_{12}\left(t\right)=\{\begin{array}{c}0,u_1\left(t\right)200,u_1\left(t\right)>460\\ \left(u_1\left(t\right)-200\right)/100,200<u_1\left(t\right)300\\ \left(460-u_1\left(t\right)\right)/160,300<u_1\left(t\right)460\end{array},& \left(10\right)\\ r_{13}\left(t\right)=\{\begin{array}{c}0,u_1\left(t\right)300,u_1\left(t\right)>1000\\ \left(u_1\left(t\right)-300\right)/160,300<u_1\left(t\right)460\\ \left(1000-u_1\left(t\right)\right)/540,460<u_1\left(t\right)1000\end{array},& \left(11\right)\\ r_{14}\left(t\right)=\{\begin{array}{c}0,u_1\left(t\right)<460\\ \left(u_1\left(t\right)-460\right)/540,460u_1\left(t\right)1000\\ 1,u_1\left(t\right)>1000\end{array},& \left(12\right)\end{array}$ the membership degrees of water pressure in different risk rank are $\begin{array}{cc}{r}_{21}\left(t\right)=\{\begin{array}{c}1,u_2\left(t\right)5\\ \left(10-u_2\left(t\right)\right)/5,5<u_2\left(t\right)10\\ 0,u_2\left(t\right)>10\end{array},& \left(13\right)\\ r_{22}\left(t\right)=\{\begin{array}{c}0,u_2\left(t\right)5,u_2\left(t\right)>15\\ \left(u_2\left(t\right)-5\right)/5,5<u_2\left(t\right)10\\ \left(15-u_2\left(t\right)\right)/5,10<u_2\left(t\right)15\end{array},& \left(14\right)\\ r_{23}\left(t\right)=\{\begin{array}{c}0,u_2\left(t\right)10,u_2\left(t\right)>20\\ \left(u_2\left(t\right)-10\right)/5,10<u_2\left(t\right)15\\ \left(20-u_2\left(t\right)\right)/5,15<u_2\left(t\right)20\end{array},& \left(15\right)\\ r_{24}\left(t\right)=\{\begin{array}{c}0,u_2\left(t\right)<15\\ \left(u_2\left(t\right)-15\right)/5,15u_2\left(t\right)20\\ 1.Math..Math.u_2\left(t\right)>20\end{array},& \left(16\right)\end{array}$ the membership degrees of aeration in different risk rank are $\begin{array}{cc}{r}_{31}\left(t\right)=\{\begin{array}{c}1,u_3\left(t\right)15\\ \left(20-u_3\left(t\right)\right)/5,15<u_3\left(t\right)20\\ 0,u_3\left(t\right)>20\end{array},& \left(17\right)\\ r_{32}\left(t\right)=\{\begin{array}{c}0,u_3\left(t\right)15,u_3\left(t\right)>30\\ \left(u_3\left(t\right)-15\right)/5,15<u_3\left(t\right)20\\ \left(30-u_3\left(t\right)\right)/10,20<u_3\left(t\right)30\end{array},& \left(18\right)\\ r_{33}\left(t\right)=\{\begin{array}{c}0,u_3\left(t\right)20,u_3\left(t\right)>50\\ \left(u_3\left(t\right)-20\right)/10,20<u_3\left(t\right)30\\ \left(50-u_3\left(t\right)\right)/20,30<u_3\left(t\right)50\end{array},& \left(19\right)\\ r_{34}\left(t\right)=\{\begin{array}{c}0,u_3\left(t\right)<30\\ \left(u_3\left(t\right)-30\right)/20,30u_3\left(t\right)50\\ 1,u_3\left(t\right)>50\end{array},& \left(20\right)\end{array}$ the membership degrees of water permeability in different risk rank are $\begin{array}{cc}{r}_{41}\left(t\right)=\{\begin{array}{c}1,y\left(t\right)30\\ \left(60-y\left(t\right)\right)/30,30<y\left(t\right)60\\ 0,y\left(t\right)>60\end{array},& \left(21\right)\\ r_{42}\left(t\right)=\{\begin{array}{c}0,y\left(t\right)30,y\left(t\right)>80\\ \left(y\left(t\right)-30\right)/30,30<y\left(t\right)60\\ \left(80-y\left(t\right)\right)/20,60<y\left(t\right)80\end{array},& \left(22\right)\\ r_{43}\left(t\right)=\{\begin{array}{c}0,y\left(t\right)60,y\left(t\right)>200\\ \left(y\left(t\right)-60\right)/20,60<y\left(t\right)80\\ \left(80-y\left(t\right)\right)/120,80<y\left(t\right)200\end{array},& \left(23\right)\\ r_{44}\left(t\right)=\{\begin{array}{c}0,y\left(t\right)<80\\ \left(y\left(t\right)-80\right)/120,80y\left(t\right)200\\ 1,y\left(t\right)>200\end{array},& \left(24\right)\end{array}$ 3) determine the pollution levels, suppose B(t)=[b.sub.1(t), b.sub.2(t), b.sub.3(t), b.sub.4(t)] is the possibility vector of matrix R(t), (t)=[.sub.1(t), .sub.2(t), .sub.3(t), .sub.4(t)] is the weight vector of matrix R(t); the relation between b.sub.j(t) and .sub.j(t) is
b.sub.j(t)=r.sub.1j(t).sub.1(t)+r.sub.2j(t).sub.2(t)+r.sub.3j(t).sub.3(t)+r.sub.4j(t).sub.4(t),(25)
B(t)=R(t)(t). (26) where b.sub.j(t), j=1,2,3,4, reflects the possibility of the jth risk rank, B(t) can reflect the contribution degree of different risk ranks; then, it will be
B(t)=(t)(t), (27)
.sub.max(t)=max (t), (28) where (t) is the ratio coefficient vector between B(t) and (t); the maximum ratio coefficient is given as .sub.max(t), which is the maximum eigenvalue of R(t); according to the matrix theory, the column ordinal of largest eigenvalue could be considered as the corresponding risk rank; membrane pollution intelligent early warning method consists of process data acquisition and pretreatment, membrane pollution intelligent prediction and early warning; the process data are collected by the acquisition instrument installed on the process site; five principal component variables is extracted by the characteristic analysis model based on partial least squares method; the soft-computing model based on recurrent fuzzy neural network is established to achieve water permeability prediction; the level of membrane fouling is evaluated by establishing a comprehensive evaluation model; the results of water permeability prediction and membrane fouling early warning are displayed in the interface of an early warning system to guide the operation of water plant, which can improve the efficiency and economic benefits of MBR wastewater treatment process.

Description

DESCRIPTION OF DRAWINGS

[0024] FIG. 1 is the overall structure of membrane pollution intelligent early warning system.

[0025] FIG. 2 is the data acquisition hardware platform of the MBR membrane fouling intelligent early warning system.

[0026] FIG. 3 is a network deployment diagram of the MBR membrane fouling intelligent early warning system.

[0027] FIG. 4 is structural topology diagram of the recurrent fuzzy neural network.

[0028] FIG. 5 is the predicted result diagram of the permeability, where the black line with the star is the desired output value of the permeability, and the black line with the point is the predicted value of the recurrent fuzzy neural network.

[0029] FIG. 6 is the prediction error diagram of the permeability.

DETAILED DESCRIPTION OF THE INVENTION

[0030] (1) Design of Membrane Pollution Intelligent Early Warning System and Implementation of Software and Hardware Function Integration

[0031] The hardware platform environment built in the actual wastewater treatment plant is shown in FIG. 2. The running process data is acquired by the acquisition instrument installed at the process site; which is transmitted to the PLC through Modbus communication protocol, and PLC transmits the running process data to the host computer through RS232 communication protocol. The data in the host computer is transmitted to the data processing server through the local area network. The running process data is displayed to the management personnel in wastewater treatment plant through the Web server in the Browser/Server mode, and the water permeability prediction and the membrane fouling early warning result are displayed in the Client/Server mode. The main functions of the developed MBR membrane fouling intelligent early warning system include: {circle around (1)}query of membrane operation parameters, {circle around (2)}online prediction of water permeability, and {circle around (3)}early warning of membrane fouling level.

[0032] The invention adopts the component technology in the software industry to package the membrane fouling data preprocessing module, the membrane fouling intelligent prediction module and the membrane fouling intelligent early warning module as functional modules, which enhances the reusability of the model, and compensates for the blank from the intelligent early warning technology of MBR membrane fouling to the human-computer interaction interface in the actual system operation at home and abroad. This invention adopts the .NET platform for software development, facilitates the creation of ActiveX controls, and expands the usable environment of the software. The fieldbus technology is used to establish a full-process system communication network to realize information transmission between modules, Meanwhile, the MBR membrane fouling intelligent early warning system realizes the connection between the central control room and the various data collection points in the field, which constitutes a centralized early warning system. The system is easy to expand, and each part has independent functions, which can add software and hardware modules according to actual predictions and integrate with other systems to achieve stability and reliability of the system and ensure the early warning accuracy of membrane fouling.

[0033] (2) Implementation of Membrane Fouling Intelligent Early Warning Method

[0034] The invention provides a membrane bioreactor-MBR membrane pollution intelligent early warning method, that the characteristic variable of the MBR membrane water permeability is obtained by feature analysis, the soft-computing model of the MBR membrane water permeability is established by recurrent fuzzy neural network to achieve the intelligent detection of MBR membrane permeability, a comprehensive evaluation model of membrane fouling level is established through the prediction values of membrane permeability combining with other process variables that can be collected by the wastewater treatment plant to realize the judgment of membrane fouling level, which improve the intelligent early warning of membrane fouling in wastewater treatment plant to ensure the normal operation of the wastewater treatment process.

[0035] {circle around (1)}The input variables are collected by the online measuring instrument installed at the process site. Five variables are acquired which parameter information and collection position are shown in Table 1.

TABLE-US-00001 TABLE 1 Process variable type Acquisition Parameter name Unit Acquisition position instrument Water Flow m.sup.3/h Head of MBR pool ViSolid700IQ Water pressure kPa End of MBR pool SensoLyt700IQ Aeration m.sup.3/h Gas pump SensoLyt700IQ Anaerobic zone ORP mV End of anaerobic tank SensoLyt700IQ Aerobic zone nitrate mg/l Secondary settling tank NitraLyt700IQ of aerobic pool end

[0036] {circle around (2)}A soft-computing model is established using recurrent fuzzy neural network. The real-time data is collected to train and test the recurrent fuzzy neural network. 80 samples are selected as testing data. The collected data is shown in Table 2.

[0037] {circle around (3)}Comprehensive evaluation of membrane fouling is established by using the predicted values of water permeability and other relevant acquisition variables (water flow, water pressure, and aeration) to obtain the pollution level of the membrane.

TABLE-US-00002 TABLE 2 Test data of soft-computing model ORP in Nitrate in Number Water Water Anaerobic Aerobic of data flow pressure Aeration zone zone sets (m.sup.3/h) (kPa) (m.sup.3/h) (mV) (mg/l) 1 356.46 28.91 5222.8 164.13 5.84 2 352.92 29.6 7815.39 163.77 5.85 3 347.29 27.86 7815.39 218.82 5.8 4 278.54 24.31 4629.63 223.31 5.84 5 334.58 27.08 4629.63 157.19 5.93 6 328.54 28.04 4629.63 162.83 5.94 7 343.75 27.95 4629.63 167.25 5.9 8 343.96 27.6 6828.7 163.48 5.96 9 332.5 26.91 4722.22 163.48 5.93 10 338.96 27.52 5740.74 154.51 5.95 11 334.79 27.52 5740.74 158.56 5.95 12 338.75 27.52 7751.74 158.56 5.99 13 308.75 26.04 7751.74 142.14 5.99 14 319.58 26.3 4560.19 109.01 6.06 15 343.75 28.3 4560.19 101.85 6.05 16 306.04 25.43 7167.25 86.81 6.16 17 296.25 24.74 7847.22 83.84 6.18 18 303.96 25.17 7847.22 95.56 6.24 19 312.71 25.69 4649.88 66.55 6.27 20 310.42 25.61 5196.76 79.28 6.31 21 320.62 26.13 7896.41 59.46 6.57 22 319.58 26.3 7896.41 50.78 6.72 23 311.67 26.04 4568.87 51.36 7.47 24 314.79 26.13 4424.19 53.24 7.57 25 317.5 26.65 5185.19 51.94 7.62 26 319.79 26.3 5185.19 58.88 7.76 27 300 26.56 5185.19 47.45 7.8 28 312.71 25.61 4539.93 52.52 7.95 29 307.5 25.61 5248.84 51 8.03 30 308.54 26.04 4629.63 52.23 8.2 31 327.92 26.91 7089.12 47.6 8.3 32 327.92 26.91 4577.55 53.24 8.32 33 329.58 26.74 4629.63 63.08 8.77 34 343.75 27.78 5280.67 60.11 8.96 35 333.75 27.95 7638.89 55.27 9.25 36 355 29.69 6973.38 73.06 9.33 37 353.12 29.08 6973.38 64.16 9.42 38 354.38 28.73 6973.38 58.09 9.46 39 340.83 27.52 6973.38 62.79 9.49 40 318.96 26.3 6973.38 55.19 9.5 41 336.67 27.86 6973.38 66.26 9.52 42 340 28.39 6973.38 56.5 9.6 43 330.42 27.6 6973.38 51.36 9.76 44 315 25.87 6973.38 49.91 9.81 45 331.04 27.34 6973.38 50.93 9.87 46 305.21 26.3 6973.38 58.02 9.67 47 313.96 26.48 6973.38 57.73 9.85 48 303.33 25.78 6973.38 62.28 9.96 49 336.25 28.39 6973.38 58.59 10.17 50 363.96 30.21 6973.38 54.9 10.13 51 299.17 25.61 6973.38 241.17 9.74 52 327.5 27.43 6973.38 231.55 9.52 53 323.12 27 5410.88 231.77 9.52 54 317.92 26.74 7989 243.42 9.35 55 317.92 26.56 7867.48 254.99 9.29 56 322.92 26.91 7867.48 262.73 9.25 57 323.33 26.74 4256.37 252.82 9.16 58 327.29 27.43 4832.18 257.6 8.68 59 318.96 26.22 7997.69 246.46 8.66 60 342.5 29.25 7798.03 283.2 8.54 61 345.62 29.08 4673.03 280.31 8.6 62 348.12 29.6 4887.15 236.11 8.57 63 331.04 28.65 7960.07 112.27 8.54 64 339.58 28.21 7968.75 68.58 8.55 65 321.88 27.17 7621.53 73.42 8.76 66 334.17 27.17 4641.2 37.25 9 67 315.42 26.48 5031.83 32.41 9.16 68 328.33 27.78 5031.83 241.68 9.06 69 394.79 31.86 8107.64 185.84 9.51 71 410.21 32.81 4620.95 246.6 9.67 72 360.83 29.25 4664.35 121.09 9.84 73 347.92 28.3 4858.22 220.49 9.82 74 359.58 28.82 8029.51 217.95 9.87 75 392.92 31.25 7873.26 238.72 9.85 76 355 28.82 6368.63 237.12 9.89 77 354.38 28.39 4265.05 226.35 9.91 78 355.21 29.69 4858.22 245.52 9.9 79 376.04 30.47 4858.22 246.09 9.91 80 363.33 29.77 7696.76 240.16 9.92