CN116165596A - Special transformer user metering abnormity monitoring method - Google Patents

Abstract

The invention provides a special transformer user metering abnormality monitoring method, which is based on the transformer external characteristic principle, and can infer whether an abnormality or electricity stealing exists in a metering device by analyzing whether the change relation of the voltage of a low-voltage side end of a transformer along with a load meets the transformer external characteristic principle. Compared with the on-site detection method, the method has higher economical efficiency and practicability, can establish a multi-target optimization model for misalignment detection according to the three-phase line voltage balance degree of the self metering point and the voltage fluctuation degree of each phase line, reduces the misjudgment or missed judgment conditions, and improves the accuracy of output results.

Description

Special transformer user metering abnormity monitoring method

Technical Field

The invention belongs to the technical field of special-purpose transformer user abnormal electricity utilization, and particularly relates to a special-purpose transformer user metering abnormality monitoring method.

Background

An important ring in electricity analysis is the search for abnormal electricity usage. The main body of abnormal electricity is electricity stealing behavior. For the power grid, the power transmission loss is increased, and the stable operation of the power grid is jeopardized. For electricity utilization enterprises, the electricity stealing behavior reduces the production cost of the enterprises, and causes unfair competition for other normal production enterprises. For all electricity utilization users, the electricity stealing behavior damages the safe operation of the transmission and distribution network, and the potential safety hazard caused by the electricity stealing behavior threatens the property safety and threatens the personal safety of residents. Therefore, the abnormal electricity utilization behavior is an urgent problem to be solved for any link of electricity utilization. And the electricity stealing behavior mining is to analyze the collected electricity consumption big data, and the data types collected by the system comprise electric quantity data, voltage, current, power factor and other operation condition data. The abnormal electricity utilization data is quickly searched, the reasons for the occurrence of the abnormal data are searched, the information delivery related departments can conduct on-site investigation in a targeted manner, excessive repeated work is avoided, the workload of on-site workers is reduced, and the operation safety of a power grid is efficiently maintained.

Disclosure of Invention

The invention aims to provide a special transformer user abnormal electricity consumption monitoring method which is based on the transformer external characteristic principle, and whether the metering device is abnormal or stolen can be estimated by analyzing whether the change relation of the voltage of the low-voltage side end of the transformer along with the load meets the transformer external characteristic principle or not.

The technical scheme is as follows:

the special-purpose variant user metering abnormality monitoring method is characterized by comprising the following steps of:

step S1: the method for collecting the operation condition data of the user metering point in the electricity consumption information collecting system at least comprises the following steps: voltage, current, active power, reactive power, voltage current phase and transformer model;

step S2: the equivalent impedance parameters R of all transformers are obtained based on the transformer parameter table by collecting the obtained transformer types _T 、X _T ；

Step S3: impedance angle solving: for the alternating acquisition measurement, directly solving the acquired voltage phase and current phase; for meter measurement, calculating an impedance angle by adopting the following formula, wherein in the formula (1), P is active power, Q is reactive power and phi is the impedance angle;

step S4: multiplying the acquired low-voltage side current measured value i by the comprehensive multiplying power k of the electric energy meter ₀ Obtaining a user true practical electric current value I:

I _a ＝k ₀ *i _a (2)

I _b ＝k ₀ *i _b (3)

I _c ＝k ₀ *i _c (4)

step S5: deducing the collected low-voltage side voltage and current to the high-voltage side based on the low-voltage side and high-voltage side vector diagram, wherein the specific deduction process is as follows:

in U _a0 、U _b0 、U _c0 U is the voltage of the primary side of the transformer after conversion _a 、U _b 、U _c For the low-voltage side three-phase voltage collected by the collection system, I _a 、I _b 、I _c Is a low-voltage side three-phase current value;

step S6: solving a three-phase voltage included angle of the high-voltage side after the zero sequence impedance is considered, and the formula is as follows:

ψ _a ＝δ _a +β _a -β _b (11)

ψ _b ＝δ _b +β _b -β _c (12)

ψ _c ＝δ _c +β _c -β _a (13)

wherein beta is _a 、β _b 、β _c For phase shift, i.e. U _i0 I=a, b, c and U _i Included angles between i=a, b, c, δ _a 、δ _b 、δ _c To collect the three-phase voltage included angle, ψ, of the system _a 、Ψ _b 、Ψ _c To consider the three-phase voltage included angle of the high-voltage side after zero sequence impedance, U _a0 、U _b0 、U _c0 An included angle is formed between the two;

step S7: bian Gong from triangle to high voltage side line voltage U _LA 、U _LB 、U _LC ：

U _LA ² ＝U _a0 ² +U _b0 ² -2U _a0 U _b0 cosψ _a (14)

U _LB ² ＝U _b0 ² +U _c0 ² -2U _b0 U _c0 cosψ _b (15)

U _LC ² ＝U _c0 ² +U _a0 ² -2U _c0 U _a0 cosψ _c (16)；

Step S8: deducing a high-voltage side voltage value from the low-voltage side measurement value by the steps S1 to S7, wherein if the measurement is abnormal or electricity is stolen, the deduced high-voltage side voltage value is deviated; suppose there is a metering anomaly or a theft of electricity, and the true current is in a linear relationship with the measured current, namely:

I _a0 ＝k _a I _a (17)

I _b0 ＝k _b I _b (18)

I _c0 ＝k _c I _c (19)

wherein I is _a0 、I _b0 、I _c0 Is the true three-phase current value, I _a 、I _b 、I _c For measuring three-phase current of measuring device, k _a 、k _b 、k _c Is a correction coefficient;

step S9: constructing a result deviation evaluation equation G (k), solving a group of correction coefficients k, minimizing G (k), if the obtained k=1, measuring the value without deviation, if k >1, measuring the value smaller than the actual value, and measuring abnormal or electricity stealing exists, and constructing the following two evaluation equations:

the first is a three-phase balancing method, which does not depend on other metering points, and assumes that the three phases of the external high-voltage line are balanced, the amplitudes are also nearly equal, and the evaluation equation is as follows:

G(k)＝(U _LA -U _LB ) ² +(U _LB -U _LC ) ² +(U _LC -U _LA ) ² (20)

u in _LA 、U _LB 、U _LC The three-phase line voltage is the metering point; when the sum of squares of the three-phase line voltages is minimum, the correction coefficient k is obtained, and because the data acquired by the acquisition system are not all ideal data, a threshold value is required to be set for the correction coefficient, and when the obtained correction coefficient exceeds the set threshold value, the measurement value of the metering device is considered to have deviation;

the second is the voltage fluctuation minimization method, which does not depend on other metering points, and the evaluation equation is as follows, assuming that the voltage of the external high-voltage line is relatively stable:

G _i (k)＝(U _LA,i -U _LA,i-1 ) ² +(U _LB,i -U _LB,i-1 ) ² +(U _LC,i -U _LC,i-1 ) ² (21)

u in _LA,i -U _LA,i-1 、U _LB,i -U _LB,i-1 、U _LC,i -U _LC,i-1 A, B, C three-phase line voltage differences between i time and i-1 time; when the sum of the voltages of the same phase line at different moments is minimum, a correction coefficient is obtained, and when the obtained correction coefficient exceeds a set threshold value, the measured value of the metering device is considered to have deviation;

step S10: establishing a multi-objective optimization model taking three-phase balance and voltage fluctuation into consideration simultaneously:

min{G(k)+G _i (k)} (22)

the misalignment gauge is positioned by thresholding the established multi-objective optimization model.

Further, in step S10, a dual threshold is set to determine the result:

1. setting a value which corresponds to the minimum multi-objective optimization model and is larger than rho as a first re-threshold, namely determining that the metering device is abnormal when the value is larger than rho, wherein the judgment rule is as follows:

min{G(k)+G _i (k)}＞ρ (23)

2. calculating 96 measuring points in one day respectively, calculating the number of measuring points with the number larger than rho corresponding to the correction coefficient k not being 1 or the minimum multi-objective optimization model, taking the number ratio of the abnormal measuring points as a second heavy threshold, and judging the abnormal measuring points according to the judgment rule:

the detection means adopted in the prior art mainly carry out field detection, and the method has low efficiency, consumes manpower and material resources, has limited detection range and has less capture quantity. Compared with the prior art, the method and the device have the advantages that an online detection method is adopted in the preferred scheme, and compared with an on-site detection method, the method and the device are higher in economical efficiency and practicability, a multi-target optimization model for misalignment detection can be built according to the three-phase line voltage balance degree of the self metering point and the voltage fluctuation degree of each phase line, the situations of misjudgment or missed judgment are reduced, and the accuracy of an output result is improved.

Drawings

The invention is described in further detail below with reference to the attached drawings and detailed description:

FIG. 1 is a graph of the normal user low side voltage current in an embodiment of the present invention.

Fig. 2 is a graph of the normal user high side voltage in an embodiment of the present invention.

FIG. 3 is a graph of voltage current on the low side of an abnormal user in accordance with an embodiment of the present invention.

FIG. 4 is a graph of an abnormal user high side voltage profile in accordance with an embodiment of the present invention.

Fig. 5 is a flow chart of a method of an embodiment of the present invention.

Fig. 6 is an equivalent circuit diagram of an embodiment tau type of the present invention.

Fig. 7 is a graph of low-voltage versus high-voltage side vectors (for example, phase a) for an embodiment of the invention.

Detailed Description

In order to make the features and advantages of the present patent more comprehensible, embodiments accompanied with figures are described in detail below:

as shown in fig. 5, the special-change user abnormal electricity consumption monitoring method provided in this embodiment mainly includes the following steps:

step one: collecting operation condition data such as voltage, current, active power, reactive power, voltage current phase, transformer model and the like of a user metering point in an electricity consumption information collecting system;

step two: the transformer equivalent circuit is tau-shaped as shown in figure 6, where R _m 、X _m For exciting impedance, R ₁ 、X _1σ For high voltage side leakage resistance, R ₂ ’、X _2σ ' is the drain impedance of the low side transition to the high side. In general, the exciting current is much smaller than the normal load current of a user, the influence on the line voltage drop is small, and R can be reduced _m 、X _m Neglecting. R is R ₁ 、X ₁ 、R ₂ ’、X _2σ ' R for whole _T 、X _T And (3) representing. The equivalent impedance parameter R of each transformer is obtained by referring to a common transformer parameter table (see annex A) through the known transformer model _T 、X _T ；

Step three: impedance angle solving: for the alternating acquisition measurement, the acquired voltage phase and current phase can be directly solved; for meter measurement, since the meter cannot collect phase information, impedance angle calculation can be performed by adopting the following formula, wherein in the formula (1), P is active power, Q is reactive power, and phi is an impedance angle;

step four: multiplying the acquired low-voltage side current measured value i by the comprehensive multiplying power k of the electric energy meter ₀ Obtaining a user true practical electric current value I:

I _a ＝k ₀ *i _a (2)

I _b ＝k ₀ *i _b (3)

I _c ＝k ₀ *i _c (4)

step five: the collected low-side voltage and current are deduced to the high-side, and as can be seen from the low-side and high-side lateral load diagrams (for example, phase a) in fig. 7, the following is a specific deduction process:

in U _a0 、U _b0 、U _c0 U is the voltage of the primary side of the transformer after conversion _a 、U _b 、U _c For the low-voltage side three-phase voltage collected by the collection system, I _a 、I _b 、I _c Is a low-voltage side three-phase current value.

Step six: solving a three-phase voltage included angle of the high-voltage side after the zero sequence impedance is considered, and the formula is as follows:

ψ _a ＝δ _a +β _a -β _b (11)

ψ _b ＝δ _b +β _b -β _c (12)

ψ _c ＝δ _c +β _c -β _a (13)

wherein beta is _a 、β _b 、β _c For phase shift, i.e. U _i0 (i=a, b, c) and U _i Included angles between (i=a, b, c), δ _a 、δ _b 、δ _c To collect the three-phase voltage included angle, ψ, of the system _a 、Ψ _b 、Ψ _c To consider the three-phase voltage included angle of the high-voltage side after zero sequence impedance, U _a0 、U _b0 、U _c0 And an included angle is formed between the two.

Step seven: bian Gong from triangle to high voltage side line voltage U _LA 、U _LB 、U _LC ：

U _LA ² ＝U _a0 ² +U _b0 ² -2U _a0 U _b0 cosψ _a (14)

U _LB ² ＝U _b0 ² +U _c0 ² -2U _b0 U _c0 cosψ _b (15)

U _LC ² ＝U _c0 ² +U _a0 ² -2U _c0 U _a0 cosψ _c (16)

Step eight: from the first to seventh steps, a high-side voltage value can be deduced from the low-side measurement value, and if the measurement is abnormal or electricity is stolen, the deduced high-side voltage value deviates. Suppose there is a metering anomaly or a theft of electricity, and the true current is in a linear relationship with the measured current, namely:

I _a0 ＝k _a I _a (17)

I _b0 ＝k _b I _b (18)

I _c0 ＝k _c I _c (19)

wherein I is _a0 、I _b0 、I _c0 Is the true three-phase current value, I _a 、I _b 、I _c For measuring three-phase current of measuring device, k _a 、k _b 、k _c Is a correction coefficient.

Step nine: a result deviation evaluation equation G (k) is constructed, a set of correction coefficients k can be solved, G (k) is minimized, if k=1 is found, the measured value has no deviation, if k >1, the measured value is smaller than the actual value, and metering abnormality or electricity theft exists. The following two evaluation equations are now constructed:

the first is a three-phase balancing method, which does not depend on other metering points, and assumes that the three phases of the external high-voltage line are balanced, the amplitudes are also nearly equal, and the evaluation equation is as follows:

G(k)＝(U _LA -U _LB ) ² +(U _LB -U _LC ) ² +(U _LC -U _LA ) ² (20)

u in _LA 、U _LB 、U _LC Is the three-phase line voltage of the metering point. When the sum of squares of the differences between the three phase voltages is minimum, the correction coefficient k is obtained, and because the data acquired by the acquisition system are not all ideal data, a threshold value is required to be determined for the correction coefficient, and when the obtained correction coefficient exceeds the set threshold value, the measurement value of the metering device is considered to have deviation.

The second is the voltage fluctuation minimization method, which does not depend on other metering points, and the evaluation equation is as follows, assuming that the voltage of the external high-voltage line is relatively stable:

G _i (k)＝(U _LA,i -U _LA,i-1 ) ² +(U _LB,i -U _LB,i-1 ) ² +(U _LC,i -U _LC,i-1 ) ² (21)

u in _LA,i -U _LA,i-1 、U _LB,i -U _LB,i-1 、U _LC,i -U _LC,i-1 The line voltages at the three phases i and i-1 of A, B, C are respectively different. When the sum of the voltages of the same phase line at different moments is minimum, a correction coefficient is obtained, and similarly, when the obtained correction coefficient exceeds a set threshold value, the measured value of the metering device is considered to have deviation.

Step ten: according to the two evaluation equations provided in the step nine, a multi-objective optimization model which simultaneously considers three-phase balance and voltage fluctuation is further established:

min{G(k)+G _i (k)} (22)

step eleven: positioning the misalignment metering device by setting a threshold value for the multi-objective optimization model established in the step ten, and judging a result by setting a double threshold value according to the evaluation model provided by the invention:

1. setting a value which corresponds to the minimum multi-objective optimization model and is larger than rho as a first re-threshold, namely determining that the metering device is abnormal when the value is larger than rho, wherein the judgment rule is as follows:

min{G(k)+G _i (k)}＞ρ (23)

2. calculating 96 measuring points in one day respectively, calculating the number of measuring points with the number larger than rho corresponding to the correction coefficient k not being 1 or the minimum multi-objective optimization model, taking the number ratio of the abnormal measuring points as a second heavy threshold, and judging the abnormal measuring points according to the judgment rule:

appendix A transformer equivalent impedance

/>

The following demonstrates the effect of the present solution by means of a specific application example:

(1) And collecting operation condition data such as voltage, current, active power, reactive power, voltage current phase, transformer model and the like of a user metering point in the electricity consumption information collecting system, and multiplying the collected current by the comprehensive multiplying power of the electric energy meter to obtain the real current of the low-voltage side of the user. Taking a normal user as an example, the real voltage-current curve of the low voltage side of the user is shown in fig. 1:

(2) After the impedance parameters and the voltage-current phase angle of the transformer are obtained, the data such as the voltage-current of the low-voltage side are deduced to the high-voltage side line voltage from the fifth step to the seventh step. Under normal conditions, the three phases of the high-voltage side three-phase line voltage are balanced, the balance degree is better than that of the low-voltage side voltage, and the deduced high-voltage side three-phase line voltage is shown in fig. 2:

(3) And (3) after obtaining the high-side three-phase line voltage, performing research and judgment through the multi-objective optimization model established in the steps eight to eleven, and judging the user as a normal user when a proper threshold is established according to experience.

(4) According to comparative analysis, the low-voltage side three-phase voltage and current and the high-voltage side three-phase voltage of an abnormal user are obtained through the same steps, and are shown in fig. 3 and 4, the unbalance degree of the high-voltage side three-phase voltage of the abnormal user can be found to be higher than that of the low-voltage side three-phase voltage, and the abnormal user can be judged to be the abnormal user through the abnormal detection method of the patent.

It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.

The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

The above description is only a preferred embodiment of the present invention, and is not intended to limit the invention in any way, and any person skilled in the art may make modifications or alterations to the disclosed technical content to the equivalent embodiments. However, any simple modification, equivalent variation and variation of the above embodiments according to the technical substance of the present invention still fall within the protection scope of the technical solution of the present invention.

The present patent is not limited to the above-mentioned best mode, any person can obtain other various special-change user metering abnormality monitoring methods under the teaching of the present patent, and all equivalent changes and modifications made according to the claims of the present application shall be covered by the present patent.

Claims (2)

1. The special-purpose variant user metering abnormality monitoring method is characterized by comprising the following steps of:

step S1: the method for collecting the operation condition data of the user metering point in the electricity consumption information collecting system at least comprises the following steps: voltage, current, active power, reactive power, voltage current phase and transformer model;

step S2: the equivalent impedance parameters R of all transformers are obtained based on the transformer parameter table by collecting the obtained transformer types _T 、X _T ；

Step S3: impedance angle solving: for the alternating acquisition measurement, directly solving the acquired voltage phase and current phase; for meter measurement, calculating an impedance angle by adopting the following formula, wherein in the formula (1), P is active power, Q is reactive power and phi is the impedance angle;

step S4: multiplying the acquired low-voltage side current measured value i by the comprehensive multiplying power k of the electric energy meter ₀ Obtaining a user true practical electric current value I:

I _a ＝k ₀ *i _a (2)

I _b ＝k ₀ *i _b (3)

I _c ＝k ₀ *i _c (4)

step S5: deducing the collected low-voltage side voltage and current to the high-voltage side based on the low-voltage side and high-voltage side vector diagram, wherein the specific deduction process is as follows:

in U _a0 、U _b0 、U _c0 U is the voltage of the primary side of the transformer after conversion _a 、U _b 、U _c For the low-voltage side three-phase voltage collected by the collection system, I _a 、I _b 、I _c Is a low-voltage side three-phase current value;

step S6: solving a three-phase voltage included angle of the high-voltage side after the zero sequence impedance is considered, and the formula is as follows:

ψ _a ＝δ _a +β _a -β _b (11)

ψ _b ＝δ _b +β _b -β _c (12)

ψ _c ＝δ _c +β _c -β _a (13)

wherein beta is _a 、β _b 、β _c For phase shift, i.e. U _i0 I=a, b, c and U _i Included angles between i=a, b, c, δ _a 、δ _b 、δ _c To collect the three-phase voltage included angle, ψ, of the system _a 、Ψ _b 、Ψ _c To consider the three-phase voltage included angle of the high-voltage side after zero sequence impedance, U _a0 、U _b0 、U _c0 An included angle is formed between the two;

step S7: bian Gong from triangle to high voltage side line voltage U _LA 、U _LB 、U _LC ：

U _LA ² ＝U _a0 ² +U _b0 ² -2U _a0 U _b0 cosψ _a (14)

U _LB ² ＝U _b0 ² +U _c0 ² -2U _b0 U _c0 co _s ψ _b (15)

U _LC ² ＝U _c0 ² +U _a0 ² -2U _c0 U _a0 cosψ _c (16)；

Step S8: deducing a high-voltage side voltage value from the low-voltage side measurement value by the steps S1 to S7, wherein if the measurement is abnormal or electricity is stolen, the deduced high-voltage side voltage value is deviated; suppose there is a metering anomaly or a theft of electricity, and the true current is in a linear relationship with the measured current, namely:

I _a0 ＝k _a I _a (17)

I _b0 ＝k _b I _b (18)

I _c0 ＝k _c I _c (19)

wherein I is _a0 、I _b0 、I _c0 Is the true three-phase current value, I _a 、I _b 、I _c For measuring three-phase current of measuring device, k _a 、k _b 、k _c Is a correction coefficient;

step S9: constructing a result deviation evaluation equation G (k), solving a group of correction coefficients k, minimizing G (k), if the obtained k=1, measuring the value without deviation, if k >1, measuring the value smaller than the actual value, and measuring abnormal or electricity stealing exists, and constructing the following two evaluation equations:

the first is a three-phase balancing method, which does not depend on other metering points, and assumes that the three phases of the external high-voltage line are balanced, the amplitudes are also nearly equal, and the evaluation equation is as follows:

G(k)＝(U _LA -U _LB ) ² +(U _LB -U _LC ) ² +(U _LC -U _LA ) ² (20)

u in _LA 、U _LB 、U _LC The three-phase line voltage is the metering point; when the sum of squares of the differences between the three phase voltages is minimum, the correction coefficient k is obtained, and since the data acquired by the acquisition system are not all ideal data, a threshold value is set for the correction coefficient, when the obtained correction coefficient exceeds the obtained correction coefficientWhen the threshold value is set, the measured value of the metering device is considered to have deviation;

the second is the voltage fluctuation minimization method, which does not depend on other metering points, and the evaluation equation is as follows, assuming that the voltage of the external high-voltage line is relatively stable:

G _i (k)＝(U _LA,i -U _LA,i-1 ) ² +(U _LB,i -U _LB,i-1 ) ² +(U _LC,i -U _LC,i-1 ) ² (21)

u in _LA,i -U _LA,i-1 、U _LB,i -U _LB,i-1 、U _LC,i -U _LC,i-1 A, B, C three-phase line voltage differences between i time and i-1 time; when the sum of the voltages of the same phase line at different moments is minimum, a correction coefficient is obtained, and when the obtained correction coefficient exceeds a set threshold value, the measured value of the metering device is considered to have deviation;

step S10: establishing a multi-objective optimization model taking three-phase balance and voltage fluctuation into consideration simultaneously:

min{G(k)+G _i (k)} (22)

the misalignment gauge is positioned by thresholding the established multi-objective optimization model.

2. The private transformer user metering anomaly monitoring method of claim 1, wherein: in step S10, a double threshold is set to determine the result:

1. setting a value which corresponds to the minimum multi-objective optimization model and is larger than rho as a first re-threshold, namely determining that the metering device is abnormal when the value is larger than rho, wherein the judgment rule is as follows:

min(G(k)+G _i (k)}＞ρ (23)

2. calculating 96 measuring points in one day respectively, calculating the number of measuring points with the number larger than rho corresponding to the correction coefficient k not being 1 or the minimum multi-objective optimization model, taking the number ratio of the abnormal measuring points as a second heavy threshold, and judging the abnormal measuring points according to the judgment rule:

/>

Images