Electrodialysis water treatment control system and control method

Abstract

The invention discloses an electrodialysis water treatment control system and control method, and relates to the technical field of electrodialysis water treatment control, and the electrodialysis water treatment control system comprises a controllable data analysis module used for analyzing the relationship between controllable data and concentration variation; obtaining optimal controllable data based on the relationship between the controllable data and the concentration variation; the processing time analysis module is used for obtaining the relationship between the solution concentration and the processing time under the optimal controllable data condition, and obtaining the optimal processing time based on the relationship between the solution concentration and the processing time; and the control module is used for carrying out electrodialysis water treatment by using the optimal controllable data and the optimal treatment time. The system is used for solving the problems that only single controllable data is analyzed and the controllable data is not screened in the existing electrodialysis water treatment technology; and therefore, the optimal control data of electrodialysis water treatment cannot be accurately found out.

Description

Technical Field

[0001] The present invention relates to the technical field of electrodialysis water treatment control, and in particular to an electrodialysis water treatment control system and a control method. Background Art

[0002] Electrodialysis is a physical and chemical process that utilizes the selective permeability of ion exchange membranes to make anions and cations in water migrate directionally and selectively permeate through the membrane under the action of an external DC electric field, thereby achieving the separation of ions from water. The treated water usually has a lower salt content and is suitable for a variety of uses.

[0003] However, there are many influencing factors in the electrodialysis water treatment process, such as current density, inlet flow rate and temperature, and these are all controllable data. Single data analysis and control cannot make the electrodialysis treatment efficiency reach the best, and the acquired control data may be single, so it is necessary to screen the data. For example, in the patent application with publication number CN118579873A, a system and control method based on electrodialysis wastewater treatment is disclosed. This scheme only analyzes the inlet flow rate data, resulting in that the electrodialysis water treatment efficiency cannot reach the highest. The existing electrodialysis water treatment technology only analyzes a single controllable data and does not screen the controllable data, resulting in the inability to accurately find the optimal control data for electrodialysis water treatment. Summary of the invention

[0004] The present invention aims to solve one of the technical problems in the prior art to at least a certain extent, by analyzing the relationship between controllable data and concentration variation; obtaining optimal controllable data based on the relationship between controllable data and concentration variation; obtaining the relationship between solution concentration and processing time under the conditions of optimal controllable data, and obtaining the optimal processing time based on the relationship between solution concentration and processing time, thereby solving the problem in the existing electrodialysis water treatment technology that only a single controllable data is analyzed and the controllable data is not screened, resulting in the inability to accurately find the optimal control data for electrodialysis water treatment.

[0005] To achieve the above-mentioned object, in a first aspect, the present invention provides an electrodialysis water treatment control system, comprising: a data acquisition module, a controllable data analysis module, a processing time analysis module and a control module;

[0006] The data acquisition module is used to acquire the main controllable data in the electrodialysis water treatment process, the main controllable data including current density and liquid inlet flow rate; marking the change in the concentration of the solution in the fresh water chamber within the first time as the concentration change;

[0007] The controllable data analysis module is used to analyze the relationship between the controllable data and the concentration variation; based on the relationship between the controllable data and the concentration variation, the optimal controllable data is obtained;

[0008] The processing time analysis module is used to obtain the relationship between solution concentration and processing time under the optimal controllable data conditions, and obtain the optimal processing time based on the relationship between solution concentration and processing time;

[0009] The control module is used to perform electrodialysis water treatment using optimal controllable data and optimal treatment time.

[0010] Furthermore, the controllable data analysis module is configured with a permutation and combination strategy, and the permutation and combination strategy includes:

[0011] Get the adjustable range of current density, marked as current density range;

[0012] Divide the current density range into N equal current density ranges, marked as current equal division ranges;

[0013] Obtain the maximum value of each current division range and mark it as the current density score;

[0014] Get the adjustable range of the inlet flow rate, marked as the inlet flow rate range;

[0015] Divide the liquid inlet flow rate range into M equal liquid inlet flow rate ranges, marked as flow rate equal division ranges;

[0016] Obtain the maximum value of each flow rate division range and mark it as the liquid inlet flow rate score;

[0017] Take one current density score and one liquid inlet flow rate score and combine them to obtain N*M combinations, which are marked as score combinations.

[0018] Furthermore, the controllable data analysis module is configured with a threshold acquisition strategy, and the threshold acquisition strategy includes:

[0019] Obtain the concentration change of the first quantity under each score combination, marked as the score concentration change;

[0020] The first screening threshold value and the second screening threshold value of the concentration change under each score combination are obtained by using the abnormal data threshold value acquisition method;

[0021] Under each score combination, the score concentration changes that are less than the first screening threshold and greater than the second screening threshold are deleted, and the score combination corresponding to the deleted score concentration changes is deleted at the same time, and the remaining current density score, liquid inlet flow score and score concentration change after deletion are marked as remaining current density score, remaining liquid inlet flow score and remaining score concentration change, respectively.

[0022] Furthermore, the abnormal data threshold acquisition method includes:

[0023] Label the first quantity Ds;

[0024] Arrange the score concentration changes from small to large, and set each score concentration change to correspond to an arrangement number, which is a positive integer;

[0025] Determine whether k1*Ds is an integer. If it is an integer, the score concentration change amount of the arrangement number k1*Ds is taken as the first score. If it is not an integer, the average of the score concentration change amounts corresponding to the arrangement numbers on both sides of k1*Ds is taken as the first score; wherein k1 is the first score ratio;

[0026] Determine whether k2*Ds is an integer. If it is an integer, the concentration change of the score with the arrangement number k2*Ds is taken as the second score. If it is not an integer, the average of the concentration change of the scores corresponding to the arrangement numbers on both sides of k1*Ds is taken as the second score; wherein k2 is the second score ratio;

[0027] The first score ratio and the second score ratio satisfy the relationship: k2 = [(1 / k1)-1]*k1;

[0028] The first screening threshold is calculated as: Dy1 = B1-b*(B2-B1); wherein Dy1 is the first screening threshold, B1 is the first score, B2 is the second score, and b is the screening range constant;

[0029] The second screening threshold is calculated as: Dy2=B2+b*(B2-B1); wherein Dy2 is the second screening threshold, B1 is the first score, B2 is the second score, and b is the screening range constant.

[0030] Furthermore, the controllable data analysis module is configured with a scatter plot drawing strategy, and the scatter plot drawing strategy includes:

[0031] A plane rectangular coordinate system is established with the current density score as the X-axis, the liquid flow score as the Y-axis, and the score concentration change as the Z-axis, and is marked as the concentration change coordinate system;

[0032] The residual current density, the residual inlet flow rate score and the residual score concentration change of each group are taken as data points in the concentration change coordinate system and marked as concentration relationship data points;

[0033] All concentration relationship data points are merged into the concentration change coordinate system to obtain a scatter plot, which is marked as a concentration relationship scatter plot.

[0034] Furthermore, the controllable data analysis module is configured with an optimal controllable data acquisition strategy, and the optimal controllable data acquisition strategy includes:

[0035] Set the concentration relationship fitting function to: Ndb = p1*Dm 2 +p2*Yl 2+p3*Dm+p4*Yl+p5;

[0036] Where Ndb is the residual current density, Dm is the residual inlet flow fraction, Yl is the residual fraction concentration change, p1, p2, p3, p4 and p5 are constants, Dm>0, Yl>0;

[0037] The concentration relationship scatter plot is fitted with a concentration relationship fitting function to obtain a concentration relationship surface;

[0038] In the concentration variation coordinate system, when the maximum value of the Z axis of the concentration relationship surface is obtained, the corresponding values ​​of the X axis and the Y axis are marked as the optimal current density and the optimal liquid flow rate, respectively.

[0039] Furthermore, the processing time analysis module is configured with a data screening strategy, and the data screening strategy includes:

[0040] Under the conditions of the second optimal current density and the optimal liquid inlet flow rate, the solution concentration is obtained at every second time interval starting from the processing time of 0 seconds;

[0041] The abnormal data screening method was used to obtain the screening solution concentration at each treatment time.

[0042] Furthermore, the abnormal data screening method includes:

[0043] Get all solution concentration ranges under a treatment time, marked as concentration range;

[0044] Divide the concentration range into z equal ranges, marked as concentration division ranges;

[0045] Count the solution concentration frequency in each concentration range and mark it as concentration frequency;

[0046] Draw a concentration frequency histogram with solution concentration as the X-axis, concentration frequency as the Y-axis, and concentration division range as the histogram interval;

[0047] In the concentration frequency histogram, judge from the left and right sides to the middle of the histogram in turn. If the concentration frequency is less than or equal to the concentration frequency threshold, delete the concentration frequency in the concentration frequency histogram. If the concentration frequency is greater than the concentration frequency threshold, stop judging.

[0048] Obtain the remaining solution concentrations in the concentration frequency histogram and mark them as screening solution concentrations.

[0049] Furthermore, the processing time analysis module is configured with a time acquisition strategy, and the time acquisition strategy includes:

[0050] A plane rectangular coordinate system is established with the treatment time as the X-axis and the screening solution concentration as the Y-axis, and is marked as the treatment time coordinate system;

[0051] The treatment time and the concentration of the screening solution are used as coordinate points of the treatment time coordinate system and are marked as time coordinate points;

[0052] All time coordinate points are imported into the processing time coordinate system to obtain a time scatter plot;

[0053] Set the time prediction fitting function to:

[0054] Where Ndg is the screening concentration, e is the natural constant, t is the treatment time, and o1 and o2 are constants;

[0055] The time scatter plot is fitted with the time prediction fitting function to obtain the time prediction fitting line and the specific values ​​of o1 and o2;

[0056] The processing time of each interval L is drawn as a tangent on the prediction fitting line;

[0057] In the processing time coordinate system, the absolute value of the tangent slope is obtained from left to right. If the absolute value of the tangent slope is less than the slope threshold, the processing time corresponding to the tangent point at this time is obtained and marked as the optimal processing time.

[0058] In a second aspect, the present invention provides an electrodialysis water treatment control method, comprising the following steps: obtaining main controllable data in the electrodialysis water treatment process, the main controllable data including current density and liquid inlet flow rate; marking the change in solution concentration in the fresh water chamber within a first time as a concentration change;

[0059] Analyze the relationship between controllable data and concentration variation; obtain optimal controllable data based on the relationship between controllable data and concentration variation;

[0060] Obtain the relationship between solution concentration and treatment time under optimal controllable data conditions, and obtain the optimal treatment time based on the relationship between solution concentration and treatment time;

[0061] Use the best controllable data and the best processing time for electrodialysis water treatment.

[0062] Beneficial effects of the present invention: The present invention analyzes the relationship between controllable data and concentration variation, obtains optimal controllable data based on the relationship between controllable data and concentration variation, obtains the relationship between solution concentration and treatment time under the condition of optimal controllable data, and obtains the optimal treatment time based on the relationship between solution concentration and treatment time. The advantages are that the current density and the inlet flow rate are analyzed in combination to obtain the settings of the current density and the inlet flow rate when the electrodialysis water treatment efficiency reaches a high level, and obtain the minimum time when the electrodialysis water treatment reaches the peak of the fresh water concentration, and obtain the minimum time when the electrodialysis water treatment is completed, thereby improving the accuracy and efficiency of finding the optimal control data for the electrodialysis water treatment.

[0063] The present invention obtains a concentration relationship surface by fitting a concentration relationship fitting function to a concentration relationship scatter plot. The advantage is that the accuracy of finding the optimal control data for electrodialysis water treatment can be improved by intuitively finding the optimal current density and the optimal liquid inlet flow rate. BRIEF DESCRIPTION OF THE DRAWINGS

[0064] Figure 1 is a functional block diagram of the system of the present invention;

[0065] Figure 2 is a schematic diagram of a concentration relationship scatter plot of the present invention;

[0066] Figure 3 is a schematic diagram of a concentration frequency histogram of the present invention;

[0067] Figure 4 is a schematic diagram of a time scatter plot of the present invention;

[0068] Figure 5 A flow chart of the steps of the present invention. DETAILED DESCRIPTION

[0069] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0070] Example 1, please refer to Figure 1 As shown, an electrodialysis water treatment control system includes: a data acquisition module, a controllable data analysis module, a processing time analysis module and a control module;

[0071] The data acquisition module is used to obtain the main controllable data in the electrodialysis water treatment process, and the main controllable data include current density and liquid inlet flow rate; the change in the solution concentration of the desalination chamber within the first time is marked as the concentration change; the first time setting here needs to be relatively small, specifically set to 5 to 20 minutes, and a time length can be arbitrarily selected within the range of 5 to 20 minutes as the first time. If the first time is too long, the change in the solution concentration is basically the same, and the electrodialysis water efficiency cannot be distinguished. For example, the first time is 10 minutes; this example is applied to the desalination of concentrated brine, and the solution concentration is expressed as 10% NaCl, which is simplified to 10% below;

[0072] The controllable data analysis module is used to analyze the relationship between the controllable data and the concentration variation; based on the relationship between the controllable data and the concentration variation, the optimal controllable data is obtained;

[0073] The controllable data analysis module is configured with permutation and combination strategies, which include:

[0074] Get the adjustable range of current density, marked as current density range;

[0075] Divide the current density range into N equal current density ranges, marked as current equal division ranges;

[0076] Obtain the maximum value of each current division range and mark it as the current density score;

[0077] Get the adjustable range of the inlet flow rate, marked as the inlet flow rate range;

[0078] Divide the liquid inlet flow rate range into M equal liquid inlet flow rate ranges, marked as flow rate equal division ranges;

[0079] Obtain the maximum value of each flow rate equal division range and mark it as the liquid inlet flow rate score;

[0080] Take one current density score and one liquid inlet flow rate score and combine them to obtain N*M combinations, which are marked as score combinations; by dividing the data, the amount of calculation is reduced and the data is made more accurate;

[0081] In practical applications, the current density range is: 0-50mA / cm 2 , divide the current density range into 5 current equal division ranges, such as 0-10mA / cm 2 With 10-20mA / cm 2 The maximum value of each current division range is obtained, that is, the current density score is: 10mA / cm 2 With 20mA / cm 2 etc., obtain the inlet flow rate range of 0-30L / h; divide the inlet flow rate range into 3 flow rate equalization ranges of 0-10L / h, 10-20L / h and 20-30L / h; obtain the maximum value of each flow rate equalization range of 10L / h, 120L / h and 30L / h, combine the current density score and the inlet flow rate score by one each, and obtain 15 score combinations, for example, a score combination current density score is 10mA / cm 2 The liquid flow rate is 10L / h.

[0082] The controllable data analysis module is configured with a threshold acquisition strategy, which includes:

[0083] Obtain the concentration change of the first quantity under each score combination, marked as the score concentration change;

[0084] The first screening threshold value and the second screening threshold value of the concentration change under each score combination are obtained by using the abnormal data threshold value acquisition method;

[0085] Under each score combination, the score concentration changes that are less than the first screening threshold and greater than the second screening threshold are deleted, and the score combination corresponding to the deleted score concentration changes is deleted at the same time, and the remaining current density score, liquid inlet flow score and score concentration change after deletion are marked as the remaining current density score, the remaining liquid inlet flow score and the remaining score concentration change respectively; delete the data that is too large or too small, which may be abnormal data;

[0086] In practical applications, for example, in a score combination, the current density score is 10mA / cm 2 When the liquid flow rate score is 10L / h, if the first number is 20, the score concentration changes are arranged from small to large as 1.21%, 1.41%, ..., 1.86%; the first screening threshold is 1.31%; the second screening threshold is 1.81%, then 1.21% and 1.86% are deleted;

[0087] The abnormal data threshold acquisition method includes:

[0088] Label the first quantity Ds;

[0089] Arrange the score concentration changes from small to large, and set each score concentration change to correspond to an arrangement number, which is a positive integer;

[0090] Determine whether k1*Ds is an integer. If it is an integer, take the score concentration change of the arrangement number k1*Ds as the first score. If it is not an integer, find the average of the score concentration changes corresponding to the arrangement numbers on both sides of k1*Ds as the first score; k1 is the first score ratio; the first score ratio ranges from 0 to 0.5, and the first score ratio is optimally set to 0.25;

[0091] Determine whether k2*Ds is an integer. If it is an integer, the score concentration change amount of the arrangement number k2*Ds is used as the second score. If it is not an integer, the average of the score concentration change amounts corresponding to the arrangement numbers on both sides of k1*Ds is taken as the second score; k2 is the second score ratio; the first score ratio range is between 0.5 and 1;

[0092] The first score ratio and the second score ratio satisfy the relationship: k2 = [(1 / k1)-1]*k1;

[0093] The first screening threshold is calculated as: Dy1 = B1-b*(B2-B1); wherein Dy1 is the first screening threshold, B1 is the first score, B2 is the second score, and b is the screening range constant;

[0094] The second screening threshold is calculated as: Dy2 = B2 + b*(B2-B1); wherein Dy2 is the second screening threshold, B1 is the first score, B2 is the second score, and b is the screening range constant; the setting of b determines the size of the first screening threshold and the second screening threshold; the setting of b is generally between 1 and 2;

[0095] In practical applications, for example, in a score combination, the current density score is 10mA / cm 2 When the liquid flow rate score is 10L / h, if the first number is 20, the score concentration changes are arranged from small to large as 1.21%, 1.41%, ..., 1.86%; set k1=0.25, k2=[(1 / k1)-1]*k1=0.75, judge k1*Ds=5 as an integer, and take the score concentration change with arrangement number 5 as 1.51% first score, judge k2*Ds=15 as an integer, and take the score concentration change with arrangement number 15 as 1.61% first score; b is set to 2, and the first screening threshold is calculated as: Dy1=B1-b*(B2-B1)=1.31%; the second screening threshold is calculated as: Dy2=B2+b*(B2-B1)=1.81%.

[0096] The controllable data analysis module is configured with a scatter plot drawing strategy, which includes:

[0097] A plane rectangular coordinate system is established with the current density score as the X-axis, the liquid flow score as the Y-axis, and the score concentration change as the Z-axis, and is marked as the concentration change coordinate system;

[0098] The residual current density, the residual inlet flow rate score and the residual score concentration change of each group are taken as data points in the concentration change coordinate system and marked as concentration relationship data points;

[0099] All concentration relationship data points are merged into the concentration change coordinate system to obtain a scatter plot, which is marked as a concentration relationship scatter plot.

[0100] The controllable data analysis module is configured with an optimal controllable data acquisition strategy, which includes:

[0101] Set the concentration relationship fitting function to: Ndb = p1*Dm 2 +p2*Yl 2 +p3*Dm+p4*Yl+p5;

[0102] Where Ndb is the residual current density, Dm is the residual inlet flow fraction, Yl is the residual fraction concentration change, p1, p2, p3, p4 and p5 are constants, Dm>0, Yl>0;

[0103] The concentration relationship scatter plot is fitted with a concentration relationship fitting function to obtain a concentration relationship surface;

[0104] In the concentration variation coordinate system, when the maximum value of the Z axis of the concentration relationship surface is obtained, the corresponding values ​​of the X axis and the Y axis are marked as the optimal current density and the optimal liquid flow rate, respectively.

[0105] In practical application, please participate Figure 2 As shown, observing the concentration relationship curve, the maximum value of the Z axis of the surface is 6%, and the corresponding values ​​of the X axis and the Y axis are 30mA / cm 2 With 20L / h, the optimal current density is 30mA / cm 2 The optimal liquid inlet flow rate is 20L / h.

[0106] The processing time analysis module is used to obtain the relationship between solution concentration and processing time under the optimal controllable data conditions, and obtain the optimal processing time based on the relationship between solution concentration and processing time;

[0107] The processing time analysis module is configured with data filtering strategies, which include:

[0108] Under the conditions of the second optimal current density and the optimal liquid inlet flow rate, the solution concentration is obtained at every second time interval starting from the processing time of 0 seconds;

[0109] The abnormal data screening method was used to obtain the screening solution concentration at each treatment time.

[0110] Abnormal data screening methods include:

[0111] Get all solution concentration ranges under a treatment time, marked as concentration range;

[0112] Divide the concentration range into z equal ranges, marked as concentration division ranges;

[0113] Count the solution concentration frequency in each concentration range and mark it as concentration frequency;

[0114] Draw a concentration frequency histogram with solution concentration as the X-axis, concentration frequency as the Y-axis, and concentration division range as the histogram interval;

[0115] In the concentration frequency histogram, judgment is made from the left and right sides to the middle of the histogram. If the concentration frequency is less than or equal to the concentration frequency threshold, the concentration frequency is deleted from the concentration frequency histogram. If the concentration frequency is greater than the concentration frequency threshold, the judgment is stopped; the concentration frequency less than or equal to the concentration frequency threshold is marked as an abnormally small frequency; the concentration frequency threshold means that the frequency distribution is small, such as the concentration frequency threshold is set to 1;

[0116] Obtain the remaining solution concentrations in the concentration frequency histogram and mark them as screening solution concentrations;

[0117] In practical application, please participate Figure 3 As shown, the concentration range of 8.1-8.6 is obtained when the processing time is 30s; the concentration range is divided into 5 identical concentration division ranges, the concentration frequency is counted, and a concentration frequency histogram is drawn. In the concentration frequency histogram, judgments are made from the left side to the middle of the histogram. If the concentration frequency 2 is greater than the concentration frequency threshold 1, the judgment is stopped; in the concentration frequency histogram, judgments are made from the right side to the middle of the histogram. If the concentration frequency 1 is equal to the concentration frequency threshold 1, the concentration frequency is deleted from the concentration frequency histogram, and the concentration frequency 0 is continued to be equal to the concentration frequency threshold 1. The concentration frequency is deleted from the concentration frequency histogram, and the concentration frequency 2 is continued to be greater than the concentration frequency threshold 1, and the judgment is stopped.

[0118] The processing time analysis module is configured with a time acquisition strategy, which includes:

[0119] A plane rectangular coordinate system is established with the treatment time as the X-axis and the screening solution concentration as the Y-axis, and is marked as the treatment time coordinate system;

[0120] The treatment time and the concentration of the screening solution are used as coordinate points of the treatment time coordinate system and are marked as time coordinate points;

[0121] All time coordinate points are imported into the processing time coordinate system to obtain a time scatter plot;

[0122] Set the time prediction fitting function to:

[0123] Where Ndg is the screening concentration, e is the natural constant, t is the treatment time, and o1 and o2 are constants;

[0124] The time scatter plot is fitted with the time prediction fitting function to obtain the time prediction fitting line and the specific values ​​of o1 and o2;

[0125] The processing time of each interval L is drawn as a tangent on the prediction fitting line;

[0126] Obtain the absolute value of the tangent slope from left to right in the processing time coordinate system. If the absolute value of the tangent slope is less than the slope threshold, obtain the processing time corresponding to the tangent point at this time and mark it as the optimal processing time. The smaller the absolute value of the tangent slope, the slower the conversion. The slope threshold is 0.1.

[0127] In practical application, please participate Figure 4As shown, the prediction fitting line is obtained by fitting. Starting from the start time, every 10 minutes, the absolute values ​​of the tangent slopes are obtained from left to right in the prediction fitting line processing time coordinate system, namely |-0.81|, |-0.71|, ..., |-0.06|. If the absolute value of the tangent slope is 0.06, which is less than the slope threshold of 0.1, the processing time corresponding to the tangent point is 100 minutes.

[0128] The control module is used to perform electrodialysis water treatment using optimal controllable data and optimal treatment time;

[0129] In practical applications, the optimal current density is 30 mA / cm2 when the solution concentration is 14%. 2 The optimal liquid inlet flow rate is 20L / h and the set time is 100min.

[0130] Example 2, please refer to Figure 5 As shown, a method for controlling electrodialysis water treatment comprises the following steps:

[0131] Step S1, obtaining the main controllable data in the electrodialysis water treatment process, the main controllable data including current density and liquid inlet flow rate; marking the change in the concentration of the solution in the fresh water chamber within the first time as the concentration change;

[0132] Step S2, analyzing the relationship between the controllable data and the concentration variation; obtaining the optimal controllable data based on the relationship between the controllable data and the concentration variation; Step S2 includes the following sub-steps:

[0133] Step S201, obtaining an adjustable range of current density, marked as a current density range;

[0134] Step S202, dividing the current density range into N equal current density ranges, marked as current equal division ranges;

[0135] Step S203, obtaining the maximum value of each current equally divided range, and marking it as a current density score;

[0136] Step S204, obtaining an adjustable range of the liquid inlet flow rate, marked as the liquid inlet flow rate range;

[0137] Step S205, dividing the liquid inlet flow rate range into M equal liquid inlet flow rate ranges, marked as flow rate equal division ranges;

[0138] Step S206, obtaining the maximum value of each flow rate equal division range, and marking it as the liquid inlet flow rate score;

[0139] Step S207, combining one current density score and one liquid inlet flow rate score to obtain N*M combinations, which are marked as score combinations;

[0140] Step S208, obtaining the concentration change of the first quantity under each score combination, marked as the score concentration change;

[0141] Step S209, using the abnormal data threshold acquisition method to obtain the first screening threshold and the second screening threshold of the concentration change under each score combination; Step S209 includes the following sub-steps:

[0142] Step S20901, marking the first quantity as Ds;

[0143] Step S20902, arranging the score concentration changes from small to large, setting each score concentration change to correspond to an arrangement number, which is a positive integer;

[0144] Step S20903, determine whether k1*Ds is an integer. If it is an integer, the score concentration change amount of the arrangement number k1*Ds is used as the first score. If it is not an integer, the average of the score concentration change amounts corresponding to the arrangement numbers on both sides of k1*Ds is calculated as the first score; wherein k1 is the first score ratio;

[0145] Step S20904, determine whether k2*Ds is an integer. If it is an integer, the score concentration change amount of the arrangement number k2*Ds is used as the second score. If it is not an integer, the average value of the score concentration change amount corresponding to the arrangement number on both sides of k1*Ds is obtained as the second score; wherein k2 is the second score ratio; the first score ratio and the second score ratio satisfy the relationship: k2 = [(1 / k1)-1]*k1;

[0146] Step S20905, calculating and obtaining the first screening threshold value: Dy1 = B1-b*(B2-B1); wherein Dy1 is the first screening threshold value, B1 is the first score, B2 is the second score, and b is the screening range constant;

[0147] Step S20906, calculate and obtain the second screening threshold value: Dy2=B2+b*(B2-B1); wherein Dy2 is the second screening threshold value, B1 is the first score, B2 is the second score, and b is the screening range constant.

[0148] Step S210, establishing a plane rectangular coordinate system with the current density score as the X-axis, the liquid flow score as the Y-axis, and the score concentration change as the Z-axis, marked as the concentration change coordinate system;

[0149] Step S211, taking the residual current density, the residual inlet flow rate score and the residual score concentration change of each group as data points in the concentration change coordinate system, and marking them as concentration relationship data points;

[0150] Step S212, integrating all concentration relationship data points into the concentration variation coordinate system to obtain a scatter plot, which is marked as a concentration relationship scatter plot;

[0151] Step S213, setting the concentration relationship fitting function to: Ndb = p1*Dm 2 +p2*Yl 2 +p3*Dm+p4*Yl+p5; where Ndb is the residual current density, Dm is the residual inlet flow fraction, Yl is the residual fraction concentration change, p1, p2, p3, p4 and p5 are constants, Dm>0, Yl>0;

[0152] Step S213, fitting the concentration relationship scatter diagram with a concentration relationship fitting function to obtain a concentration relationship surface;

[0153] Step S214, in the concentration variation coordinate system, the corresponding values ​​of the X-axis and the Y-axis when the maximum value of the Z-axis of the concentration relationship surface is obtained, and are marked as the optimal current density and the optimal liquid inlet flow rate respectively.

[0154] Step S3, obtaining the relationship between the solution concentration and the processing time under the optimal controllable data conditions, and obtaining the optimal processing time based on the relationship between the solution concentration and the processing time; Step S3 includes the following sub-steps:

[0155] Step S301, under the conditions of the second optimal current density and the optimal liquid inlet flow rate, the solution concentration is obtained at every second time interval starting from the processing time of 0 seconds;

[0156] Step S302, using an abnormal data screening method to obtain the concentration of the screening solution at each processing time; Step S302 includes the following sub-steps:

[0157] Step S30201, obtaining all solution concentration ranges under a processing time, marked as concentration ranges;

[0158] Step S30202, dividing the concentration range into z identical ranges, marked as concentration division ranges;

[0159] Step S30203, counting the solution concentration frequency of each concentration division range, and marking it as concentration frequency;

[0160] Step S30204, drawing a concentration frequency histogram with solution concentration as the X-axis, concentration frequency as the Y-axis, and concentration division range as the histogram interval;

[0161] Step S30205, in the concentration frequency histogram, judging is performed from the left and right sides to the middle of the histogram in sequence. If the concentration frequency is less than or equal to the concentration frequency threshold, the concentration frequency is deleted from the concentration frequency histogram. If the concentration frequency is greater than the concentration frequency threshold, the judging is stopped.

[0162] Step S30205, obtaining the remaining solution concentrations in the concentration frequency histogram and marking them as screening solution concentrations.

[0163] Step S303, establishing a plane rectangular coordinate system with the processing time as the X-axis and the concentration of the screening solution as the Y-axis, marked as the processing time coordinate system;

[0164] Step S304, taking the processing time and the concentration of the screening solution as coordinate points of the processing time coordinate system, and marking them as time coordinate points;

[0165] Step S305, importing all time coordinate points into the processing time coordinate system to obtain a time scatter plot;

[0166] Step S306, setting the time prediction fitting function to: Where Ndg is the screening concentration, e is the natural constant, t is the treatment time, and o1 and o2 are constants;

[0167] Step S307, fitting the time scatter plot with a time prediction fitting function to obtain a time prediction fitting line and specific values ​​of o1 and o2;

[0168] Step S308, draw a tangent line on the prediction fitting line for each processing time interval L;

[0169] Step S309, obtain the absolute value of the tangent slope from left to right in the processing time coordinate system. If the absolute value of the tangent slope is less than the slope threshold, obtain the processing time corresponding to the tangent point at this time and mark it as the optimal processing time.

[0170] Step S4, using the optimal controllable data and the optimal processing time to perform electrodialysis water treatment.

[0171] It should be understood by those skilled in the art that the embodiments of the present invention can be provided as methods, systems or computer program products. Therefore, the present invention can take the form of a complete hardware embodiment, a complete software embodiment or an embodiment combining software and hardware. Moreover, the present invention can take the form of a computer program product implemented on one or more computer-usable storage media containing computer-usable program codes. Among them, the storage medium can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (Static Random Access Memory, referred to as SRAM), electrically erasable programmable read-only memory (Electrically Erasable Programmable Read-Only Memory, referred to as EEPROM), erasable programmable read-only memory (Erasable Programmable Read Only Memory, referred to as EPROM), programmable read-only memory (Programmable Red-Only Memory, referred to as PROM), read-only memory (Read-Only Memory, referred to as ROM), magnetic memory, flash memory, disk or optical disk. These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0172] In the embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. The device embodiments described above are merely schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some communication interfaces, and the indirect coupling or communication connection of devices or units can be electrical, mechanical or other forms.

Claims

1. An electrodialysis water treatment control system, characterized in that: include: Data acquisition module, controllable data analysis module, processing time analysis module and control module; The data acquisition module is used to acquire the main controllable data in the electrodialysis water treatment process, the main controllable data including current density and liquid inlet flow rate; marking the change in the concentration of the solution in the fresh water chamber within the first time as the concentration change; The controllable data analysis module is used to analyze the relationship between the controllable data and the concentration variation; Obtaining optimal controllable data based on the relationship between controllable data and concentration variation; The processing time analysis module is used to obtain the relationship between solution concentration and processing time under the optimal controllable data conditions, and obtain the optimal processing time based on the relationship between solution concentration and processing time; The control module is used to perform electrodialysis water treatment using optimal controllable data and optimal treatment time.

2. An electrodialysis water treatment control system according to claim 1, characterized in that: The controllable data analysis module is configured with a permutation and combination strategy, and the permutation and combination strategy includes: Get the adjustable range of current density, marked as current density range; Divide the current density range into N equal current density ranges, marked as current equal division ranges; Obtain the maximum value of each current division range and mark it as the current density score; Get the adjustable range of the inlet flow rate, marked as the inlet flow rate range; Divide the liquid inlet flow rate range into M equal liquid inlet flow rate ranges, marked as flow rate equal division ranges; Obtain the maximum value of each flow rate division range and mark it as the liquid inlet flow rate score; Take one current density score and one liquid inlet flow rate score and combine them to obtain N*M combinations, which are marked as score combinations.

3. An electrodialysis water treatment control system according to claim 2, characterized in that: The controllable data analysis module is configured with a threshold acquisition strategy, and the threshold acquisition strategy includes: Obtain the concentration change of the first quantity under each score combination, marked as the score concentration change; The first screening threshold value and the second screening threshold value of the concentration change under each score combination are obtained by using the abnormal data threshold value acquisition method; Under each score combination, the score concentration changes that are less than the first screening threshold and greater than the second screening threshold are deleted, and the score combination corresponding to the deleted score concentration changes is deleted at the same time, and the remaining current density score, liquid inlet flow score and score concentration change after deletion are marked as remaining current density score, remaining liquid inlet flow score and remaining score concentration change, respectively.

4. An electrodialysis water treatment control system according to claim 3, characterized in that: The abnormal data threshold acquisition method includes: Label the first quantity Ds; Arrange the score concentration changes from small to large, and set each score concentration change to correspond to an arrangement number, which is a positive integer; Determine whether k1*Ds is an integer. If it is an integer, the score concentration change amount of the arrangement number k1*Ds is taken as the first score. If it is not an integer, the average of the score concentration change amounts corresponding to the arrangement numbers on both sides of k1*Ds is taken as the first score; wherein k1 is the first score ratio; Determine whether k2*Ds is an integer. If it is an integer, the concentration change of the score with the arrangement number k2*Ds is taken as the second score. If it is not an integer, the average of the concentration change of the scores corresponding to the arrangement numbers on both sides of k1*Ds is taken as the second score; wherein k2 is the second score ratio; The first score ratio and the second score ratio satisfy the relationship: k2 = [(1 / k1)-1]*k1; The first screening threshold is calculated as: Dy1 = B1-b*(B2-B1); wherein Dy1 is the first screening threshold, B1 is the first score, B2 is the second score, and b is the screening range constant; The second screening threshold is calculated as: Dy2=B2+b*(B2-B1); wherein Dy2 is the second screening threshold, B1 is the first score, B2 is the second score, and b is the screening range constant.

5. An electrodialysis water treatment control system according to claim 4, characterized in that: The controllable data analysis module is configured with a scatter plot drawing strategy, and the scatter plot drawing strategy includes: A plane rectangular coordinate system is established with the current density score as the X-axis, the liquid flow score as the Y-axis, and the score concentration change as the Z-axis, and is marked as the concentration change coordinate system; The residual current density, the residual inlet flow rate score and the residual score concentration change of each group are taken as data points in the concentration change coordinate system and marked as concentration relationship data points; All concentration relationship data points are merged into the concentration change coordinate system to obtain a scatter plot, which is marked as a concentration relationship scatter plot.

6. An electrodialysis water treatment control system according to claim 5, characterized in that: The controllable data analysis module is configured with an optimal controllable data acquisition strategy, and the optimal controllable data acquisition strategy includes: Set the concentration relationship fitting function to: Ndb = p1*Dm 2 +p2*Yl 2 +p3*Dm+p4*Yl+p5; Where Ndb is the residual current density, Dm is the residual inlet flow fraction, Yl is the residual fraction concentration change, p1, p2, p3, p4 and p5 are constants, Dm>0, Yl>0; The concentration relationship scatter plot is fitted with a concentration relationship fitting function to obtain a concentration relationship surface; In the concentration variation coordinate system, when the maximum value of the Z axis of the concentration relationship surface is obtained, the corresponding values ​​of the X axis and the Y axis are marked as the optimal current density and the optimal liquid flow rate, respectively.

7. An electrodialysis water treatment control system according to claim 6, characterized in that: The processing time analysis module is configured with a data screening strategy, which includes: Under the conditions of the second optimal current density and the optimal liquid inlet flow rate, the solution concentration is obtained at every second time interval starting from the processing time of 0 seconds; The abnormal data screening method was used to obtain the screening solution concentration at each treatment time.

8. An electrodialysis water treatment control system according to claim 7, characterized in that: Abnormal data screening methods include: Get all solution concentration ranges under a treatment time, marked as concentration range; Divide the concentration range into z equal ranges, marked as concentration division ranges; Count the solution concentration frequency in each concentration range and mark it as concentration frequency; Draw a concentration frequency histogram with solution concentration as the X-axis, concentration frequency as the Y-axis, and concentration division range as the histogram interval; In the concentration frequency histogram, judge from the left and right sides to the middle of the histogram in turn. If the concentration frequency is less than or equal to the concentration frequency threshold, delete the concentration frequency in the concentration frequency histogram. If the concentration frequency is greater than the concentration frequency threshold, stop judging. Obtain the remaining solution concentrations in the concentration frequency histogram and mark them as screening solution concentrations.

9. An electrodialysis water treatment control system according to claim 8, characterized in that: The processing time analysis module is configured with a time acquisition strategy, and the time acquisition strategy includes: A plane rectangular coordinate system is established with the treatment time as the X-axis and the screening solution concentration as the Y-axis, and is marked as the treatment time coordinate system; The treatment time and the concentration of the screening solution are used as coordinate points of the treatment time coordinate system and are marked as time coordinate points; All time coordinate points are imported into the processing time coordinate system to obtain a time scatter plot; Set the time prediction fitting function to: Where Ndg is the screening concentration, e is the natural constant, t is the treatment time, and o1 and o2 are constants; The time scatter plot is fitted with the time prediction fitting function to obtain the time prediction fitting line and the specific values ​​of o1 and o2; The processing time of each interval L is drawn as a tangent on the prediction fitting line; In the processing time coordinate system, the absolute value of the tangent slope is obtained from left to right. If the absolute value of the tangent slope is less than the slope threshold, the processing time corresponding to the tangent point at this time is obtained and marked as the optimal processing time.

10. An electrodialysis water treatment control method, applicable to an electrodialysis water treatment control system according to any one of claims 1 to 9, characterized in that: The steps include: Obtain the main controllable data in the electrodialysis water treatment process, including current density and liquid inlet flow rate; mark the change in solution concentration in the fresh water chamber within the first time as the concentration change; Analyze the relationship between controllable data and concentration changes; Obtaining optimal controllable data based on the relationship between controllable data and concentration variation; Obtain the relationship between solution concentration and treatment time under optimal controllable data conditions, and obtain the optimal treatment time based on the relationship between solution concentration and treatment time; Use the best controllable data and the best processing time for electrodialysis water treatment.

Images