Comprehensive Harmonic Responsibility Calculation Based on Different Weighting Methods

Dividing the responsibility of harmonics between customer and utility sides is necessary to improve the objectivity, rationality, and scientificalness of power measurement and power quality evaluation. On the basis of expert experience and customer needs and considering the background harmonic fluctuation, this study proposes two reasonable conditions to evaluate harmonic responsibility. At the same time, the reference index of harmonic responsibility calculation and the problem of different frequency harmonics are considered in the comprehensive calculation. On the basis of the IEEE Std. 1459–2010 power theory, the index set of harmonic responsibility is established, and several common subjective and objective weighting methods are used to weigh each index. On the basis of optimization theory, an optimization model is established by constructing a Lagrange function in finding the condition extremum to unify the subjective and objective information. Finally, the calculation method is verified by a Norton equivalent model, and the harmonic responsibility of each index is calculated by the harmonic current vector method; thus, calculating the actual impedance of customer and utility sides is unnecessary, making the calculation simple and effective. Results of the comparative analysis show that the comprehensive evaluation method of harmonic responsibility with the combinatorial weighting method can not only meet the different needs of different groups for harmonic indicators but also reflect the background harmonic fluctuations sensitively. In this manner, the subjective and objective information are unified, and the shortcomings of a single weighting method are overcome. Hence, the evaluation results are accurate, practical, and reasonable.


Introduction
With the rapid development of power electronics and nonlinear loads, harmonics in power systems have gradually become serious, thereby attracting increasing scientific research attention. First, the harmonics have been analyzed in detail [1]. With the development of the subject, researchers have proposed several positioning methods of the main harmonic source [2][3][4]; however, these methods only locate the main harmonic source, but they cannot calculate the harmonic responsibility of each harmonic source quantitatively. Reference [5] first proposed the calculation of harmonic responsibility in 2002; the model of calculation of harmonic responsibility is shown in Figure 1. With the development of this topic, many researchers begin to focus on modeling approach of utility side and customer side [6][7][8], According to calculation model of harmonic responsibility, which is based on the Norton equivalent model shown in Figure 1, the calculation of parameters of the utility side and customer side is obviously very important, thus proposing fluctuation method, binary linear regression, and so on. Compared with harmonic current vector method, these methods show some practical deficiencies.
The traditional calculation method of harmonic responsibility mainly focuses on the contributions of the harmonic current and the harmonic voltage at the point of common coupling (PCC). Therefore, Reference [16] proposed the index of harmonic current and harmonic voltage based on superposition theory. Some researchers have explored the comprehensive evaluation of harmonic responsibility and obtained some results in the evaluation method. In Reference [17], the PCC point current was divided into distorted and non-distorted currents, which were used as the basis for the division of harmonic current at the utility and customer sides, but only applicable to the case where the customer side was pure resistance load. In Reference [18], the load current was divided into two parts, namely, nonlinear and linear currents, which were used as the basis for the division of the harmonic responsibility of the customer side. The quantitative evaluation index was provided; however, this method was unsuitable for the condition of utility side distortion. In Reference [19], the principles of harmonic current index and harmonic voltage index were analyzed in detail, and the two principles were relatively different; thus, the calculation results of harmonic responsibility were also different. Reference [20] indicated that the fluctuation of harmonic voltage and current is large in the case of background harmonic fluctuation. Moreover, the calculation results of existing methods may have large deviation, and they cannot reflect the variation characteristics of harmonic voltage and current nor consider the effect of background harmonic fluctuation. In addition, the traditional index based on the superposition theory was used to calculate the harmonic responsibility under the specific harmonic frequency, which lacks further research on the combination of different frequencies. Therefore, the relevant staff often is often confused when selecting the harmonic responsibility calculation index due to the lack of unity [19]. Currently, Emanuel proposed the latest power theory standard IEEE Std. 1459-2010, which defines different power physical quantities under different conditions [21] and provides a standard for relevant workers to measure physical quantities under harmonic conditions. Similar to the calculations of harmonic responsibility, the comprehensive evaluation of power quality also needs to combine different indicators [22][23][24][25], which provide enlightenment for comprehensive harmonic responsibility calculation.
For the responsibility evaluation, determining the weight that corresponds to the evaluation index is important. The change of the weight value has a direct effect on the rationality of the comprehensive evaluation and thus on the comprehensive evaluation of harmonic responsibility. In view of the aforementioned problems, this study proposes the idea of the comprehensive harmonic responsibility calculation, which can evaluate the harmonic responsibility of each of the harmonic sources more comprehensively and accurately compared to the previously published work. On the basis of the actual situation, this study presents two reasonable judgment conditions for the According to calculation model of harmonic responsibility, which is based on the Norton equivalent model shown in Figure 1, the calculation of parameters of the utility side and customer side is obviously very important, thus proposing fluctuation method, binary linear regression, and so on. Compared with harmonic current vector method, these methods show some practical deficiencies.
The traditional calculation method of harmonic responsibility mainly focuses on the contributions of the harmonic current and the harmonic voltage at the point of common coupling (PCC). Therefore, Reference [16] proposed the index of harmonic current and harmonic voltage based on superposition theory. Some researchers have explored the comprehensive evaluation of harmonic responsibility and obtained some results in the evaluation method. In Reference [17], the PCC point current was divided into distorted and non-distorted currents, which were used as the basis for the division of harmonic current at the utility and customer sides, but only applicable to the case where the customer side was pure resistance load. In Reference [18], the load current was divided into two parts, namely, nonlinear and linear currents, which were used as the basis for the division of the harmonic responsibility of the customer side. The quantitative evaluation index was provided; however, this method was unsuitable for the condition of utility side distortion. In Reference [19], the principles of harmonic current index and harmonic voltage index were analyzed in detail, and the two principles were relatively different; thus, the calculation results of harmonic responsibility were also different. Reference [20] indicated that the fluctuation of harmonic voltage and current is large in the case of background harmonic fluctuation. Moreover, the calculation results of existing methods may have large deviation, and they cannot reflect the variation characteristics of harmonic voltage and current nor consider the effect of background harmonic fluctuation. In addition, the traditional index based on the superposition theory was used to calculate the harmonic responsibility under the specific harmonic frequency, which lacks further research on the combination of different frequencies. Therefore, the relevant staff often is often confused when selecting the harmonic responsibility calculation index due to the lack of unity [19]. Currently, Emanuel proposed the latest power theory standard IEEE Std. 1459-2010, which defines different power physical quantities under different conditions [21] and provides a standard for relevant workers to measure physical quantities under harmonic conditions. Similar to the calculations of harmonic responsibility, the comprehensive evaluation of power quality also needs to combine different indicators [22][23][24][25], which provide enlightenment for comprehensive harmonic responsibility calculation.
For the responsibility evaluation, determining the weight that corresponds to the evaluation index is important. The change of the weight value has a direct effect on the rationality of the comprehensive evaluation and thus on the comprehensive evaluation of harmonic responsibility. In view of the aforementioned problems, this study proposes the idea of the comprehensive harmonic responsibility calculation, which can evaluate the harmonic responsibility of each of the harmonic sources more comprehensively and accurately compared to the previously published work. On the basis of the actual situation, this study presents two reasonable judgment conditions for the comprehensive calculation method of the harmonic responsibility based on different weighting methods, establishes the indicator set of harmonic responsibility based on IEEE Std. 1459-2010 power theory, and selects several commonly used subjective and objective weighting methods. At the same time, on the basis of optimization theory, the subjective and objective weighting methods are combined as the combinatorial weighting method (CWM) for assigning the evaluation indicators.
The remainder of this paper is organized as follows. Section 2 presents the process of comprehensive harmonic responsibility calculation. Section 3 introduces and analyzes different methods of empowerment. Section 4 establishes an example of harmonic responsibility calculation using Matlab/Simpower system. The simulation results verify the practicability of the different weighting methods in the comprehensive harmonic responsibility calculation. Section 5 presents the results of the CWM on the basis of the proposed two reasonable judgment conditions.

Comprehensive Harmonic Responsibility Calculation
In practice, different harmonic responsibility indicators have different effects on power systems and users. For example, voltage distortion has a considerable effect on power systems; thus, the customer side should take greater responsibility for distorted voltage. The injected harmonic current also has a more serious effect on the customer side. Thus, the utility side should take greater responsibility for the distorted current, and the opposite applies to the customer side. In practice, the fifth harmonic problem is likely to be caused by the system, whereas the third harmonic problem is likely to be caused by the user. Thus, the customer side should bear less responsibility for the fifth harmonic and greater responsibility for the third harmonic. By contrast, under the background harmonic fluctuation condition, the fluctuation of harmonic voltage and current on the utility and customer sides are relatively large, and the harmonic parameters in the actual power grid change with time. Therefore, evaluation results should reflect the changes of data sensitively. In sum up, the following conditions should be fully considered for the reasonable division of harmonic liability: Condition 1: The evaluation results should meet the different needs of customer side and utility side for harmonic indicators.
Condition 2: The evaluation results should reflect the change of background harmonic sensitively. Therefore, the comprehensive calculation steps of harmonic responsibility are as follows. First, the index set of harmonic responsibility calculation is established based on IEEE Std. 1459-2010. In IEEE Std. 1459-2010, the apparent power is decomposed in detail, and the equation can be expressed as follows: where U H and I H are the RMS values of harmonic voltage and harmonic current, respectively; U 1 and I 1 are the RMS values of fundamental voltage and fundamental current, respectively; S 1 is The physical quantities corresponding to each index can be expressed as follows: where DVh is the voltage distortion power under hth harmonic; DIh is the current distortion power under hth harmonic; SHh is the harmonic apparent power under hth harmonic, respectively. Because the content of SH in SN is the smallest in practice, only the product of harmonic voltage and current at the same frequency is considered.
Second, the calculation model of harmonic contribution is established. On the basis of the model of harmonic responsibility in Figure 1, the calculation model of the harmonic contribution of the utility and customer sides based on superposition theory on the PCC point can be established in Figure 3.

Utility side
Customer side On the basis of the index set established in Figure 2 and the calculation model of the harmonic contribution of the utility and customer sides, the harmonic responsibility based on each index in the index set of utility/customer side is calculated, respectively. The physical quantities corresponding to each index can be expressed as follows: where D Vh is the voltage distortion power under hth harmonic; D Ih is the current distortion power under hth harmonic; S Hh is the harmonic apparent power under hth harmonic, respectively. Because the content of S H in S N is the smallest in practice, only the product of harmonic voltage and current at the same frequency is considered. Second, the calculation model of harmonic contribution is established. On the basis of the model of harmonic responsibility in Figure 1, the calculation model of the harmonic contribution of the utility and customer sides based on superposition theory on the PCC point can be established in Figure 3. Equation (1) indicates that IEEE Std. 1450-2010 separates the fundamental and harmonic parts. Nonfundamental apparent power SN represents all physical quantities that are relevant to harmonics.
According to the decomposition of apparent power in IEEE Std. 1459-2010, the index set of harmonic responsibility calculation is established in Figure 2. The physical quantities corresponding to each index can be expressed as follows: where DVh is the voltage distortion power under hth harmonic; DIh is the current distortion power under hth harmonic; SHh is the harmonic apparent power under hth harmonic, respectively. Because the content of SH in SN is the smallest in practice, only the product of harmonic voltage and current at the same frequency is considered.
Second, the calculation model of harmonic contribution is established. On the basis of the model of harmonic responsibility in Figure 1, the calculation model of the harmonic contribution of the utility and customer sides based on superposition theory on the PCC point can be established in Figure 3. On the basis of the index set established in Figure 2 and the calculation model of the harmonic contribution of the utility and customer sides, the harmonic responsibility based on each index in the index set of utility/customer side is calculated, respectively.  On the basis of the index set established in Figure 2 and the calculation model of the harmonic contribution of the utility and customer sides, the harmonic responsibility based on each index in the index set of utility/customer side is calculated, respectively.
Finally, the comprehensive harmonic responsibility is calculated based on different weighting methods. The weight W i of each index in the index set is calculated by different weighting methods, which will be introduced in detail in Section 3. Then, the comprehensive harmonic responsibility of the utility/customer side is calculated by combining the calculation results of the harmonic responsibility and the weight of each index, respectively, as shown as follows: where R i−u and R i−c are the harmonic responsibilities based on each index of the utility and customer sides, respectively; W i−u and W i−c are the weights of each index of the utility and customer sides, respectively; R * u and R * c are the initial values of the comprehensive harmonic responsibility of the utility and customer sides, respectively; R u and R c are the comprehensive harmonic responsibility of the utility and customer sides, respectively; and n is the number of index i.

Different Weighting Methods of Comprehensive Harmonic Responsibility Calculation
Equations (8)- (11) indicate that the calculation of the weight of each index W i is an important step and has a deep influence on the result of the comprehensive harmonic responsibility calculation of the utility and customer sides. Therefore, several common weighting methods are introduced as follows.

Analytic Hierarchy Process (AHP))
AHP was proposed by American mathematician T.L. Saaty. It is used to calculate the subjective weight of each index in the index set [26][27][28]. This method can avoid consistency checking using the calculation process of a judgment matrix to reduce the complexity of calculation. However, this approach is difficult to be convinced due to its small quantitative data and several qualitative components. The calculation steps are as follows.
First, the judgment matrix is calculated. At the same layer of indexes, indexes are placed in order; two adjacent indexes are compared with each other. The element b ij of the judgment matrix B represents the importance of index i to index j. Table 1 lists the values of element b ij . Table 1. Value of b i j of judgment matrix B.

Value of b ij
Verbal Judgments of Preferences Judgment matrix B can be expressed as: On the basis of the calculation of judgment matrix B, the subjective weight W i−AHP of index i can be expressed as: where b ij is the element at row i and column j of B, and n is the number of indexes.

Expert Consultation Method (ECM)
ECM is a type of market approach simulation that uses experts' knowledge, experience, and analytical judgment to verify a reference index. ECM is used to calculate the subjective weight of each index in the index set. This method is simple and practicable and is convenient in application. However, it is greatly influenced by subjective factors that may affect the accuracy of the conclusion [29]. Four steps are involved in ECM.
First, the original weight of each index in the index set is calculated. On the basis of personal experience and understanding, some experts are invited to evaluate the weight of each index; these weights are called original weights. The sum of the original weights given by the same expert should be equal to 1: where o li represents the original weight of index i given by expert l. Second, the average weight of index i is calculated: where o i represents the average weight of index i, and d represents the total number of experts. Third, the offset o * li of weight of index i is calculated: Fourth, the modified weight W i−ECM of index i is calculated by o * li and u i , as shown as follows:

Entropy Weighting Method (EWM)
EWM was proposed by Clausius. It is used to calculate the objective weight of each index in the index set. The entropy value is obtained by calculation; the larger it is, the higher the disorder of the event will be, and thus, the smaller the contribution to the final evaluation result will be [30,31]. By contrast, the higher the order of the event is, the greater the contribution to the final evaluation result will be. EWM complies with objective rules and has strict mathematical significance. However, it is not necessarily in line with actual needs because the subjective intentions of policy makers are often neglected. The calculation steps are as follows. First, an evaluation of each index in the index set can be used to obtain fuzzy mapping. On the basis of the definition of fuzzy transformation, fuzzy mapping can determine a fuzzy relation F called fuzzy evaluation matrix.
where f nm represents the fuzzy relation between index n and rank m.
On the basis of the calculation of judgment matrix F, H i , expressed as Equation (18), is the entropy of index expressed as Equation (19), which is the sum of the entropy of each index.
where m is the number of ranks of each index. The objective weight W i−EWM of index i can be expressed as:

Variation Coefficient Method (VCM)
VCM is an objective weighting method, and it is similar to the entropy weighting method. If the variation of index i is large, then different ranks of index i can be distinguished well, and the weight should then be large [32,33]. This method can eliminate the influence of dimension and measurement scale, but it lacks attention to the practical significance of indexes. Four steps are used in the application of VCM.
On the basis of fuzzy evaluation matrix F, which is expressed as Equation (14), the first step is to calculate the average value of index i: The second step is to calculate the standard deviation of index i: The third step is to calculate the value of variation of index i: The fourth step is to calculate the weight W i−VCM of index i:

Combinational Weighting Method (CWM)
The basic idea of CWM is to combine the weights that are calculated by the single weighting method mentioned above, considering subjective preferences of policy makers and the objectivity of index data, thereby achieving the unity of subjective and objective information. Common combination methods include addition combinations and multiplication combinations. The addition combination method is extremely inflexible to distinguish the good and bad information from the subjective and objective weight information. The multiplication combination method will cause a "multiplication effect", which is only applicable to the case of having many indicators and uniform weight distributions. Therefore, on the basis of the idea that the deviation between the weights of the CWM and the single weighting method is as small as possible, the following model is established: where ω i is the weight of index i after combination, R is the evaluation matrix established on the harmonic responsibility of each index, n is the number of index i, 0 < α < 1 is the coefficient of subjective preference, and 0 < 1 − α < 1 is the coefficient of objective preference. By setting the value of α, the analyst's preference for subjective and objective weights can be adjusted flexibly.
The Lagrange function is established to solve the model as follows: On the basis of the necessary conditions for the existence of extremes, the following is set: which can be simplified as: where: By solving the matrix equation (Equation (29)), the following results can be obtained:

Example and Simulation Analysis
As an example, the simulation model is established in Figure 4. The utility side contains a certain amount of background harmonics, whereas the customer side is an uncontrollable rectifier circuit. The parameters of the simulation model are set in Table 2. where: By solving the matrix equation (Equation (29)), the following results can be obtained:

Example and Simulation Analysis
As an example, the simulation model is established in Figure 4. The utility side contains a certain amount of background harmonics, whereas the customer side is an uncontrollable rectifier circuit. The parameters of the simulation model are set in Table 2.

Parameter
Comprehensively considering the third and fifth harmonics, the index set is established based on Figure 2. Then, on the basis of the measurement data on the PCC point, the parameters of the Norton equivalent model of the utility and customer sides are obtained by the improved impedance Comprehensively considering the third and fifth harmonics, the index set is established based on Figure 2. Then, on the basis of the measurement data on the PCC point, the parameters of the Norton equivalent model of the utility and customer sides are obtained by the improved impedance method [15], and the reference impedance of the utility side is selected as the actual impedance in the simulation. . . .
. I Uh are the harmonic current sources under the hth harmonic of the customer and utility sides, respectively.
The harmonic responsibility of different indexes in the index set can be calculated on the basis of the calculation of utility/customer side parameters. Take the 3rd harmonic as an example, and the responsibility based on index 1, 3, 5 can be calculated as follows: (49) Table 3 shows the calculation results. The harmonic responsibilities based on different indexes in the index set should be obtained together to calculate the harmonic responsibility of the utility/customer side. Thus, the calculation of each index weight is important.
In practice, the harmonic parameters change with time under the background harmonic fluctuation; thus, harmonic responsibility calculation should be considered under the background harmonic fluctuation. The cases selected for this study are shown in Table 4 to enable the behavioral analysis of the methods. In this case, the random number model is used to simulate the fluctuation of the customer side, and it can generate a uniformly distributed random number with a variance of 3 and a mean of 0, which is shown in the Appendix A. On the basis of the maximum and minimum corresponding physical quantities of each index in the index set, each index is divided into four levels on average within observation time, which is assumed to be 1 s. On the basis of the above level classification standard, the time statistics of the four power indexes that appear in the corresponding levels can be obtained in Table 5. On the basis of Table 5, the fuzzy evaluation matrix T case1 can be obtained as: AHP method: The index set is established based on Figure 2 for the application of the AHP method. Referring to the experts' experience and users' requirements, the power supply company and the relevant power users are involved in the division of harmonic responsibility simultaneously. Indexes i are in the following order from 1 to 6 as follows: D V3 , D V5 , D I3 , D I5 , S H3 , and S H5 . The following are also set: t 12 = 1.6, t 23 = 1.8, t 34 = 1.8, t 45 = 1.8, and t 56 = 1.2. On the basis of Equation (12), judgment matrix B, which depends on the decision maker, is expressed as follows: Then, based on Equations (51) and (13), the weight W i−c−AHPcase1 of each index of the customer side is shown in Table 6. ECM: For the use of ECM, six experts are invited to evaluate the weight of each index, and indexes are in the following evaluation order: D V3 , D V5 , D I3 , D I5 , S H3 , and S H5 . Then, the original weight matrix O can be assumed as: Then, on the basis of Equations (52) and (14)- (19), the weight W i−c−ECMcase1 of each index of the customer side is shown in Table 6.
EWM: On the basis of the fuzzy evaluation matrix, T case1 is obtained as Equation (50). When using the EWM based on Equations (21)-(23), the weight W i−c−EWMcase1 of each index of the customer side is shown in Table 7. VCM: On the basis of the fuzzy evaluation matrix, T case1 is obtained as Equation (50). On the basis of Equations (24)- (27), the weight W i−c−VCMcase1 of each index of the customer side is shown in Table 7.
CWM: When using the CWM based on Equations (28)-(39) on the basis of the weight calculated above (α = 0.5 indicates that subjective and objective weights are of equal importance), the weight W i−c−CWMcase1 of each index of customer side is shown in Table 8. The weight calculation procedure of each index on the utility side is the same as that on the customer side and will thus not be repeated here. Then, the weight of each index for Case 1 on the utility side based on the subjective, objective, and combinational weighting methods is shown in Tables 6-8, respectively. The comprehensive harmonic responsibility of the utility/customer side for Case 1 can be obtained on the basis of Equations (8)- (11), as shown in Table 9.

Case 2
On the basis of Case 1, assuming that background harmonics under indexes 1 and 5 fluctuate violently, the time for them to be in the fourth stage is 1 s and that in other stages is 0 s during the entire measurement period. Then, the fuzzy evaluation matrix T case2 is given as follows: AHP and ECM: These methods are subjective weighting methods that are qualitatively weighted. They are weighted based on the importance of each index to the power grid and customer. Thus, the calculation results of Case 2 will not change because only the background harmonic situation changes. Taking the customer side as an example, that is EWM and VCM: This approach is based on T case2 , in which the fuzzy evaluation matrix is shown in Equation (53). Then, on the basis of Equations (20)- (27), the weights of each index on the customer side by the entropy weight method (W i−c−EWMcase2 ) and the coefficient of variation method (W i−c−VCMcase2 ) are shown in Table 7.
CWM: By using the CWM based on Equations (28)-(39) on the basis of the weight calculated above, the weight (W i−c−CWMcase2 ) of each index of the customer side is shown in Table 8.
The weight calculation procedure of each index on the utility side is the same as that on the customer side and will thus not be repeated here. Then, the weight of each index of Case 2 on the utility side based on the subjective, objective, and combinational weighting methods is shown in Tables 6-8, respectively. On the basis of Equations (8)- (11), the comprehensive harmonic responsibility of the utility and customer sides for Case 1 can be obtained, as shown in Table 9.

Comparisons of Calculation Results Based on Different Weighting Methods
On the basis of the subjective weighting method, the weight of each index on the customer and utility sides are calculated, as shown in Table 6.
AHP and ECM: Table 6 shows that the weights of the customer and utility sides based on AHP and ECM are different. The distribution trend is consistent with the actual effect of harmonics on the customer and the utility sides. Moreover, given that these subjective weighting methods are based on the experts' comprehensive evaluation of the effect of each index in various aspects from their level of understanding and experiment, the subjective weighting method meets Condition 1. However, in Table 6, the weight of each index corresponding to Case 2 based on these methods has not changed any more than that of Case 1. That is to say, the weights obtained by these methods are not sensitive to the change of harmonics in a power system; thus, it cannot meet Condition 2.
On the basis of the objective weighting method, the weight of each index on the customer and utility sides are calculated, as shown in Table 7.  Table 7 shows that among the objective weights based on EWM and VCM, the weights corresponding to indexes 1 and 5 account for a large proportion in Case 2, which is obviously different from Case 1. This condition is consistent with the presupposition. Thus, the objective weighting method meets Condition 2. However, the weight distribution of the customer side is exactly the same as that of the utility side, which is not consistent with the actual demand of the customer and utility sides for each index. Therefore, the entropy weight method and the coefficient of variation method cannot reflect the different demands of power systems and users for harmonic indexes. Therefore, Condition 1 cannot be satisfied.
On the basis of CWM, the weight of each index on the customer and utility sides are calculated, as shown in Table 8. CWM: The weight in Table 8 shows that the weights distribution based on CWM in Cases 1 and 2 are consistent with the different influences of each index on the customer and utility sides in the actual situation. Thus, the rationality Condition 1 is satisfied. Furthermore, in Case 2, when indexes 1 and 5 change considerably, the corresponding weights change remarkably compared with Case 1. The change trend is consistent with the present situation. Therefore, the CWM can satisfy Conditions 1 and 2.
The analysis of comprehensive harmonic responsibility is obtained by different weighting methods.
On the basis of the weights obtained by different weighting methods, combined with Equations (5)-(8), the harmonic responsibility of the customer/utility side are calculated, as shown in Table 9.
Taking Case 1 as an example, Table 9 shows that in the single weighting method, little difference of harmonic responsibility value is obtained by a similar weighting method. By contrast, the difference obtained by the different weighting methods is large. The calculation results obtained by the CWM are between them. Therefore, selecting the appropriate weighting method is of considerable importance for the accuracy of harmonic responsibility calculation.
The harmonic responsibility of Cases 1 and 2 in Table 9 show that the calculation results based on the subjective weighting method are identical, whereas the harmonic responsibility of Case 2 based on the objective and the CWM is obviously different from that of Case 1. Therefore, neglecting the fluctuation of background harmonic may result in a certain amount of error.
The evaluation results are not sensitive to the changes of index parameters, and the evaluation results are greatly affected by subjective factors. Although objective weighting can reflect this change sensitively, it cannot reflect the different demands of different power users for harmonic indexes. It may give a large weight for some indexes that are considered unimportant by researchers. The CWM proposed in this study combines these subjective and objective weighting methods. It can reflect the subjective will, meet the different needs of different groups for harmonic indexes, and combine with the objective reality. Thus, the evaluation results can reflect the data fluctuation more sensitively, make the subjective and objective information unified, overcome the shortcomings of the single weighting method, and increase the accuracy of the evaluation results. The proposed method is thus more accurate, reasonable, and practical.
In addition to the typical calculation methods proposed in this study, other methods can be the bases of CWMs. Thus, aiming at differences in calculation results based on different CWMs, the next step is to take measurements to combine different calculation results.

Conclusions
This study considers the contradiction between the results of harmonic responsibility caused by different traditional indexes and the harmonics of different frequency harmonics. This study proposes the idea of a comprehensive calculation of harmonic responsibility. The reasonable conditions for the division of the two harmonic responsibilities are proposed on the basis of the different requirements of the utility and customer sides for the division of harmonic responsibilities and the parameter changes caused by the background harmonic fluctuations. On the basis of this condition, the application of different weighting methods in the comprehensive calculation of harmonic responsibilities was analyzed. The main contributions of this study are summarized as follows.

1.
Given the actual influence of harmonics on power supply and consumption and the harmonic characteristics, the reasonable judgment conditions for the division method of two harmonic responsibilities are proposed. It provides a basis for judging whether the calculation method of harmonic responsibility is reasonable.

2.
In comparison with the traditional calculation method of harmonic responsibility, the comprehensive calculation method of harmonic responsibility based on the latest power theory, IEEE Std. 1459-2010, can evaluate the harmonic responsibility. This method combines different frequencies and different indexes under harmonic conditions. The method can evaluate the harmonic responsibility comprehensively and avoid disputes.

3.
On the basis of the proposed reasonable conditions, several commonly used subjective and objective weighting methods are used, and the CWM is established based on optimization theory. The reasonable weight of each index is used for synthesizing the information brought by each index.

4.
Specific simulation cases verify the practicality of the methods proposed in this study, and they also show that calculation results based on the fact that different weighting methods are different.
On the basis of the comparative analysis of the proposed rationality conditions, the CWM is more reliable than the single weighting method. The former not only considers that the weighting results should also change with when the background harmonic changes dynamically, but also enables experts and users to participate in the evaluation process, thereby meeting the needs of users in different situations. The weighting value obtained is more practical and flexible. At the same time, the combination of various weighting methods can overcome the shortcomings of a single weighting method and make the evaluation results reasonable and objective.
The key point of this study is the comprehensive calculation of harmonic responsibility. On this basis, it can be applied to the adjustment scheme or harmonic economic loss assessment considering harmonic responsibility sharing, thereby making it more reasonable and perfect.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
In Section 4.2, after adding the model of a random number that is used to simulate the fluctuation of the customer side in the simulation, measuring waveforms of corresponding physical quantities of indexes within the observation time are shown in Figure A1.

Appendix A
In Section 4.2, after adding the model of a random number that is used to simulate the fluctuation of the customer side in the simulation, measuring waveforms of corresponding physical quantities of indexes within the observation time are shown in Figure A1.