Clearing method for inter-provincial auxiliary service

Abstract

A clearing method for inter-provincial auxiliary service belongs to the technical field of service clearing, and comprises the following steps: receiving an inter-provincial auxiliary service value of an inter-provincial device, the inter-provincial auxiliary service value comprising a plurality of subgroups; determining an inter-provincial auxiliary service value break variable of each subgroup; the coordination quantity of each sampling point in the inter-province auxiliary service value is calculated; calculating a Z score of each sampling point according to a Z-score method, multiplying the coordination amount by the Z score to obtain a last Z score of each sampling point, and performing detection according to the last Z score; by calculating the coordination quantity of each sampling point, the correlation degree between the sampling point and the other sampling point is determined, and the wrong determination condition which is often formed when the sampling point with small standard deviation is detected by the Z-fractional method is prevented, so that the correctness and the credibility of the stability verification of the inter-province auxiliary service value are improved.

Description

Technical Field

[0001] The present invention belongs to the technical field of service clearing, and in particular relates to a clearing method for inter-provincial auxiliary services. Background Art

[0002] During the development of electricity market, two types of ancillary service clearing methods have been experienced: sequential method and co-optimization method.

[0003] In the clearing method for inter-provincial ancillary services, the existing technical solution with patent publication number "CN111509700A" is currently generally used to achieve it, which includes constructing an inter-provincial electricity spot market clearing model, where the system load balance constraint of the electricity spot market clearing model is based on the output of inter-provincial coal-fired power units in a set period of time, the output of gas-fired power units in a set period of time, the output of non-market units in a set period of time, the power of inter-provincial interconnection lines in a set period of time, the total load of inter-provincial interconnection lines in a set period of time, the output of inter-provincial coal-fired power units in a set period of time, the output of gas-fired power units in a set period of time, the output of non-market units in a set period of time or the power of inter-provincial interconnection lines in a set period of time is the inter-provincial ancillary service value.

[0004] On the other hand, the auxiliary service also needs to perform stability checks on the inter-provincial auxiliary service values, that is, first find out the abnormal values ​​in the inter-provincial auxiliary service values, specifically obtain the abnormal values ​​through the Z-score method, and compare them with the previous abnormal values ​​to achieve the stability check of the inter-provincial auxiliary service values. The key to the Z-score method is: by giving each information point (that is, a single inter-provincial auxiliary service value) a Z score based on the standard deviation, if the absolute value of the Z score of an information point exceeds three, it means that the information point is an abnormal point.

[0005] However, the current Z-score method is key to obtaining the standard deviation of information points, without involving the correlation between information points. In the range of small standard deviation, even if an information point has some correlation with other information points (just like, it often belongs to a hidden subcluster), it is often misidentified as an abnormal value by the Z-score method due to its small standard deviation. In addition, during the process of obtaining inter-provincial auxiliary service values, the inter-provincial auxiliary service values ​​obtained often have abnormal values ​​due to factors such as scenes. Such abnormal values ​​caused by such factors as scenes will disrupt the precise calculation of the correlation between information points and the standard deviation, thereby weakening the abnormal point detection function of the Z-score method. Summary of the invention

[0006] In order to solve the defects in the prior art, the present invention proposes a clearing device and method for inter-provincial auxiliary services, which collects the inter-provincial auxiliary service value of the inter-provincial device, and the inter-provincial auxiliary service value contains multiple subgroups; identifies the inter-provincial auxiliary service value mutation amount of each subgroup; calculates the coordination amount of each sampling point in the inter-provincial auxiliary service value; calculates the Z score of each sampling point according to the Z-score method, multiplies the coordination amount by the Z score to obtain the final Z score of each sampling point, and performs detection based on the final Z score; by calculating the coordination amount of each sampling point, the correlation between the sampling point and other sampling points is identified, thereby preventing the Z-score method from often forming a wrong determination when performing detection on sampling points with small standard deviations, thereby improving the accuracy and reliability of the stability check of the inter-provincial auxiliary service value.

[0007] The present invention utilizes the following technical solutions.

[0008] A method for clearing inter-provincial auxiliary services, comprising:

[0009] Construct an inter-provincial electricity spot market clearing model, where the system load balance constraint of the electricity spot market clearing model is based on the output of inter-provincial coal-fired power units in a set period of time, the output of gas-fired power units in a set period of time, the output of non-market units in a set period of time, the power of inter-provincial interconnection lines in a set period of time, and the total load of inter-provincial interconnection lines in a set period of time. The output of inter-provincial coal-fired power units in a set period of time, the output of gas-fired power units in a set period of time, the output of non-market units in a set period of time or the power of inter-provincial interconnection lines in a set period of time is the inter-provincial auxiliary service value;

[0010] The clearing methods for inter-provincial auxiliary services also include:

[0011] Step 1: collecting the inter-provincial auxiliary service value of the inter-provincial device, where the inter-provincial auxiliary service value contains multiple subgroups;

[0012] Step 2: Identify the inter-provincial ancillary service value mutation of each subgroup;

[0013] Step 3: Calculate the coordination amount of each sampling point within the inter-provincial auxiliary service value;

[0014] Step 4: Calculate the Z score of each sampling point according to the Z-score method, multiply the coordination amount by the Z score to obtain the final Z score of each sampling point, and perform detection according to the final Z score.

[0015] Furthermore, in step 1, the inter-provincial device is an inter-provincial coal-fired power unit, a gas-fired power unit or an inter-provincial interconnection line.

[0016] Furthermore, in step 1, in order to sample the inter-provincial auxiliary service value of the inter-provincial device in real time, a power transmitter can be installed on the cable where the inter-provincial device is located, so that when the inter-provincial device is an inter-provincial coal-fired power unit, a gas-fired power unit or an inter-provincial interconnection line, samples can be taken to obtain the output of the inter-provincial coal-fired power unit in a set period of time, the output of the gas-fired power unit in a set period of time, the output of the non-marketized unit in a set period of time or the power of the inter-provincial interconnection line in a set period of time.

[0017] Further, in step 1, the sampling speed of the power transmitter is 500 times / s, so as to sample the inter-provincial auxiliary service value of the inter-provincial device, and then transfer the sampled inter-provincial auxiliary service value to the information table, and then perform service value sorting on the sampled inter-provincial auxiliary service value.

[0018] Furthermore, in step 1, the method for sorting the service values ​​is the Wiener filtering algorithm.

[0019] Furthermore, in step 1, the sampled inter-provincial auxiliary service values ​​are cut into a plurality of subgroups according to the time interval of 300 seconds formed between the sampling moments as a cutting unit.

[0020] Furthermore, in step 2, the mutation amount of the inter-provincial auxiliary service value of each subgroup is obtained through the Laplace transform method, that is, the Laplace transform method is used to perform time-frequency conversion on the sampled inter-provincial auxiliary service values, so as to separate the inter-provincial auxiliary service values ​​into phasors of different frequencies, so as to analyze the phasors outside the industrial frequency in the inter-provincial auxiliary service values, and then obtain the amplitude of each phasor outside the industrial frequency and the amplitude of the phasor at the industrial frequency, so as to respectively regard them as the actual amount of each phasor outside the industrial frequency and the actual amount of the phasor at the industrial frequency.

[0021] Furthermore, in step 2, the inter-provincial auxiliary service value mutation amount of each subgroup is determined based on the actual amount of each phase quantity outside the power frequency and the actual amount of the phase quantity at the power frequency;

[0022] The calculation equation of the inter-provincial auxiliary service value mutation is:

[0023]

[0024] In the equation, h γ represents the inter-provincial auxiliary service value mutation of the γth subgroup, V i Represents the actual quantity of the phasor outside the power frequency with power i, where i = 2, 3, 4..., V1 represents the actual quantity of the phasor of the power frequency.

[0025] Furthermore, in step 2, the method for obtaining V1 may also be: calculating the average of the phase quantities of the power frequency of all sampling points in a subgroup, and the average is the actual quantity of the phase quantity of the power frequency.

[0026] Furthermore, in step 3, the operation equation of the coordination amount is:

[0027]

[0028] In the equation, η γ，j represents the coordination amount of the jth sampling point in the γth subgroup; θ γ represents the amount obtained by subtracting the inter-provincial auxiliary service value mutation of the γth subgroup from the inter-provincial auxiliary service value mutation of the other subgroups; ν (j,l) represents the value of the j-th sampling point in the γ-th subgroup minus the value of the l-th sampling point, p represents the total number of additional sampling points in the γ-th subgroup other than the j-th sampling point, and e represents the Euler number.

[0029] Furthermore, in step 4, when performing detection based on the final Z score of each sampling point, a stable critical value is pre-defined for each sampling point in the same subgroup. When detecting whether a sampling point is abnormal, the final Z score of the sampling point is compared with the stable critical value corresponding to the sampling point. If the final Z score of the sampling point is greater than its stable critical value, then the sampling point is an abnormal point, otherwise it is a reasonable point. Then, all the abnormal points in the subgroup are totaled. If the total number of abnormal points in a subgroup exceeds half of the total number of sampling points, the maintainer is informed to perform maintenance.

[0030] A clearing device for inter-provincial auxiliary services, comprising:

[0031] A subgroup module, which is used to collect the inter-province auxiliary service value of the inter-province device, and the inter-province auxiliary service value contains multiple subgroups;

[0032] A recognition module, which is used to recognize the mutation amount of the inter-provincial auxiliary service value of each subgroup;

[0033] A calculation module, which is used to calculate the coordination amount of each sampling point in the inter-provincial auxiliary service value;

[0034] The detection module is used to calculate the Z score of each sampling point according to the Z-score method, multiply the coordination amount by the Z score to obtain the final Z score of each sampling point, and perform detection according to the final Z score.

[0035] The beneficial effects of the present invention are that, compared with the prior art, the technical effects of the present invention include:

[0036] By calculating the coordination value of each sampling point, the correlation between the sampling point and other sampling points is determined, avoiding the misdetermination often caused by the Z-score method when detecting sampling points with small standard deviations, thereby improving the accuracy and reliability of the stability check of inter-provincial auxiliary service values. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] Figure 1 is a partial flow chart of the clearing method for inter-provincial auxiliary services described in the present invention;

[0038] Figure 2 It is a partial structural diagram of the clearing device for inter-provincial auxiliary services described in the present invention. DETAILED DESCRIPTION

[0039] In order to make the purpose, technical solution and advantages of the present invention clearer, the following will be combined with the drawings in the embodiments of the present invention to clearly and completely express the technical solution of the present invention. The embodiments expressed in this application are only some embodiments of the present invention, not all embodiments. According to the spirit of the present invention, other embodiments obtained by ordinary technicians in this field without making creative work are all within the protection scope of the present invention.

[0040] like Figure 1 As shown, a method for clearing inter-provincial auxiliary services according to the present invention includes:

[0041] Construct an inter-provincial electricity spot market clearing model, where the system load balance constraint of the electricity spot market clearing model is obtained based on the output of inter-provincial coal-fired power units in a set period of time, the output of gas-fired power units in a set period of time, the output of non-market units in a set period of time, the power of inter-provincial interconnection lines in a set period of time, and the total load of inter-provincial interconnection lines in a set period of time. The output of inter-provincial coal-fired power units in a set period of time, the output of gas-fired power units in a set period of time, the output of non-market units in a set period of time, or the power of inter-provincial interconnection lines in a set period of time is the inter-provincial auxiliary service value; the system load balance constraint of the electricity spot market clearing model is achieved based on the output of inter-provincial coal-fired power units in a set period of time, the output of gas-fired power units in a set period of time, the output of non-market units in a set period of time, the power of inter-provincial interconnection lines in a set period of time, and the total load of inter-provincial interconnection lines in a set period of time, such as the corresponding method in the prior art solution with patent publication number "CN111509700A".

[0042] The clearing methods for inter-provincial auxiliary services also include:

[0043] Step 1: collecting the inter-provincial auxiliary service value of the inter-provincial device, where the inter-provincial auxiliary service value contains multiple subgroups;

[0044] In a preferred but non-limiting embodiment of the present invention, in step 1, the inter-provincial device is an inter-provincial coal-fired power unit, a gas-fired power unit or an inter-provincial connecting line.

[0045] In a preferred but non-limiting embodiment of the present invention, in step 1, the inter-provincial auxiliary service value of the inter-provincial device is to be sampled in real time, and a power transmitter can be installed on the cable where the inter-provincial device is located, so that when the inter-provincial device is an inter-provincial coal-fired power unit, a gas-fired power unit or an inter-provincial interconnection line, samples can be taken to obtain the output of the inter-provincial coal-fired power unit in a set time period, the output of the gas-fired power unit in a set time period, the output of the non-marketized unit in a set time period or the power of the inter-provincial interconnection line in a set time period.

[0046] In a preferred but non-limiting embodiment of the present invention, in step 1, the sampling speed of the power transmitter is 500 times / s, so as to perform sampling on the inter-provincial auxiliary service value of the inter-provincial device to prevent insufficient sampling and insufficient trend of the inter-provincial auxiliary service value, and then the sampled inter-provincial auxiliary service value is transferred to the information table, and then the sampled inter-provincial auxiliary service value is sorted. By performing service value sorting on the sampled inter-provincial auxiliary service value, abnormal values ​​in the inter-provincial auxiliary service value that are often caused by sampling errors, clutter disturbances, etc. are cleared, thereby preventing abnormal values ​​from often disturbing the calculation value of the Z-score method, thereby improving the accuracy of the Z-score method.

[0047] In a preferred but non-limiting embodiment of the present invention, in step 1, the method for sorting the service values ​​is the Wiener filtering algorithm.

[0048] Abnormal values ​​can be removed or interpolated to complete the value, such as replacing the abnormal value with the average of a pair of reasonable values ​​adjacent to the abnormal value.

[0049] Performing service value sorting can ensure the accuracy and credibility of inter-provincial auxiliary service values, thereby improving the accuracy and credibility of subsequent inter-provincial auxiliary service value mutation calculations and preventing noise from disrupting the Z-score method during detection.

[0050] In a preferred but non-limiting embodiment of the present invention, in step 1, the sampled inter-provincial auxiliary service values ​​are cut into a plurality of subgroups according to the time interval of 300 seconds formed between the sampling moments as a cutting unit.

[0051] Step 2: Identify the inter-provincial ancillary service value mutation of each subgroup;

[0052] In a preferred but non-limiting embodiment of the present invention, in step 2, the mutation amount of the inter-provincial auxiliary service value of each subgroup is obtained via the Laplace transform method, that is, the Laplace transform method is used to perform time-frequency conversion on the sampled inter-provincial auxiliary service value, so as to separate the inter-provincial auxiliary service value into phasors of different frequencies, so as to analyze the phasors outside the industrial frequency in the inter-provincial auxiliary service value, and then obtain the amplitude of each phasor outside the industrial frequency and the amplitude of the phasor of the industrial frequency, so as to respectively regard them as the actual amount of each phasor outside the industrial frequency and the actual amount of the phasor of the industrial frequency.

[0053] In a preferred but non-limiting embodiment of the present invention, in step 2, the inter-provincial auxiliary service value mutation amount of each subgroup is determined based on the actual amount of each phase quantity outside the power frequency and the actual amount of the phase quantity at the power frequency;

[0054] The calculation equation of the inter-provincial auxiliary service value mutation is:

[0055]

[0056] In the equation, h γ represents the inter-provincial auxiliary service value mutation of the γth subgroup, V i Represents the actual quantity of the phasor outside the power frequency with power i, where i = 2, 3, 4..., V1 represents the actual quantity of the phasor of the power frequency.

[0057] In a preferred but non-limiting embodiment of the present invention, in step 2, the method for obtaining V1 may also be: calculating the average of the power frequency phasor of all sampling points in a subgroup, and the average is the actual value of the power frequency phasor.

[0058] By calculating the mutation amount of the inter-provincial auxiliary service value of a subgroup, the error amount of the inter-provincial auxiliary service value in the subgroup can be grasped. That is to say, the larger the mutation amount of the inter-provincial auxiliary service value of a subgroup, the higher the error of the inter-provincial auxiliary service value of the subgroup from the theoretical inter-provincial auxiliary service value of the subgroup, and the higher the difference between the sampling points, that is, the lower the correlation.

[0059] In addition, the larger the mutation amount of the inter-provincial auxiliary service value, the larger the phase quantity of the inter-provincial device where the inter-provincial auxiliary service value is sourced, which means that there are often hidden inter-provincial device or working defects. Therefore, by calculating the mutation amount of the inter-provincial auxiliary service value of the inter-provincial device from time to time, the defects can be found in time and the corresponding methods can be used to perform advance maintenance or overhaul, thereby reducing the maintenance costs caused by the failure of the inter-provincial device.

[0060] Step 3: Calculate the coordination amount of each sampling point in the inter-provincial auxiliary service value; the sampling point is the inter-provincial auxiliary service value sampled at the corresponding sampling time.

[0061] In this application, the coordination amount is used to represent the degree of association between a sampling point and other sampling points in its subgroup.

[0062] In a preferred but non-limiting embodiment of the present invention, in step 3, the operation equation of the coordination amount is:

[0063]

[0064] In the equation, η γ，j represents the coordination amount of the jth sampling point in the γth subgroup; θ γrepresents the amount obtained by subtracting the inter-provincial auxiliary service value mutation of the γth subgroup from the inter-provincial auxiliary service value mutation of the other subgroups; ν (j,l) represents the value of the j-th sampling point in the γ-th subgroup minus the value of the l-th sampling point, p represents the total number of additional sampling points in the γ-th subgroup other than the j-th sampling point, and e represents the Euler number.

[0065] And θ γ represents the inter-provincial auxiliary service value mutation of the γth subgroup minus the inter-provincial auxiliary service value mutation of the other subgroups, that is, θ γ Can be: θ γ =h γ -h l (γ≠l), where h γ represents the inter-provincial auxiliary service value mutation of the γth subgroup, h l Represents the mutation amount of inter-provincial auxiliary service value of the lth subgroup, and the lth subgroup is a subgroup within the other subgroups above.

[0066] From the above coordination equation, we can understand that The higher the amount of γ，j The lower the amount, the The higher the value, the higher the difference between the jth sampling point in the γth subgroup and the other sampling points in the γth subgroup, that is, the lower the correlation between the jth sampling point and the other sampling points in the γth subgroup, the lower the coordination value of the jth sampling point. In addition, ∑θ γ The higher the amount of γ，j The lower the value of γ The higher the value, the higher the difference between the inter-provincial auxiliary service value mutation amount of the γth subgroup and the inter-provincial auxiliary service value mutation amount of other subgroups, that is, the difference between the sampling points in the γth subgroup is higher than the difference between the sampling points of other subgroups, which means that the correlation between the sampling points in the γth subgroup is lower than the correlation between the sampling points of other subgroups.

[0067] In short, a sampling point The amount or ∑θ γ The higher the amount, the lower the coordination amount of the sampling point.

[0068] Step 4: Calculate the Z score of each sampling point according to the Z-score method, multiply the coordination amount by the Z score to obtain the final Z score of each sampling point, and perform detection according to the final Z score.

[0069] The above Z score is the score obtained based on the Z-score method.

[0070] In a preferred but non-limiting embodiment of the present invention, in step 4, when detection is performed based on the final Z score of each sampling point, a stable critical value is pre-defined for each sampling point in the same subgroup. When detecting whether a sampling point is abnormal, the final Z score of the sampling point is compared with the stable critical value corresponding to the sampling point. If the final Z score of the sampling point is greater than its stable critical value, then the sampling point is an abnormal point, otherwise it is a reasonable point. Then, all the abnormal points in the subgroup are totaled. If the total number of abnormal points in a subgroup exceeds half of the total number of sampling points, it means that the inter-provincial device in the subgroup often has unstable defects, and the maintainer is informed to perform maintenance.

[0071] Since the correlation between sampling points is taken into account when calculating the final Z-score, the Z-score method is prevented from often misdetermining the auxiliary service values ​​between provinces with small standard deviations, thereby improving the accuracy of the information stability check.

[0072] The final Z scores of each sampling point in the same subgroup are accumulated, and then a stable critical value is defined in advance. If the total final Z score of each sampling point in the subgroup is greater than the pre-defined stable critical value, it is determined that the inter-provincial device in the subgroup is unstable, and the maintainer is informed to perform maintenance.

[0073] like Figure 2 As shown, a clearing device for inter-provincial auxiliary services according to the present invention comprises:

[0074] A subgroup module, which is used to collect the inter-province auxiliary service value of the inter-province device, and the inter-province auxiliary service value contains multiple subgroups;

[0075] A recognition module, which is used to recognize the mutation amount of the inter-provincial auxiliary service value of each subgroup;

[0076] A calculation module, which is used to calculate the coordination amount of each sampling point in the inter-provincial auxiliary service value;

[0077] The detection module is used to calculate the Z score of each sampling point according to the Z-score method, multiply the coordination amount by the Z score to obtain the final Z score of each sampling point, and perform detection according to the final Z score.

[0078] The beneficial effects of the present invention are that, compared with the prior art, the technical effects of the present invention include:

[0079] By calculating the coordination value of each sampling point, the correlation between the sampling point and other sampling points is determined, avoiding the misdetermination often caused by the Z-score method when detecting sampling points with small standard deviations, thereby improving the accuracy and reliability of the stability check of inter-provincial auxiliary service values.

[0080] The present disclosure can be a system, method and / or computer program product. The computer program product can include a computer readable backup medium carrying computer readable program instructions for causing a processor to achieve each aspect of the disclosure.

[0081] The computer readable backup medium can be a tangible network capable of retaining and backing up the instructions used by the instruction execution network. The computer readable backup medium can be, but is not limited to, an electrical backup network, a magnetic backup network, an optical backup network, an electromagnetic backup network, a semiconductor backup network, or any suitable combination thereof. Further examples of computer readable backup media (non-exhaustive list) include: portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), static random access memory (SRAM), portable compact disk read-only memory (HD-ROM), digital versatile disk (DXD), memory stick, floppy disk, mechanically encoded network circuits, such as punch cards or protrusions in grooves with instructions backed up thereon, and any suitable combination thereof. Computer readable backup media as used herein is not to be construed as transient signals per se, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through waveguides or other transmission media (e.g., light pulses through power line cables), or electrical signals transmitted through wires.

[0082] The computer readable program instructions expressed herein can be downloaded from the computer readable backup medium to each inference / processing grid line, or downloaded to an external computer or external backup grid line through a wireless network, such as the Internet, a local area network, a wide area network and / or a wireless network. The wireless network can include copper transmission cables, power transmission lines, wireless transmission, routers, firewalls, switches, gateway computers and / or edge service devices. The wireless network adapter card or wireless network interface in each inference / processing grid line receives the computer readable program instructions from the wireless network and forwards the computer readable program instructions for storage in the computer readable backup medium in each inference / processing grid line.

[0083] The computer program instructions for performing the operations of the present disclosure can be assembler instructions, instruction set architecture (ISA) instructions, machine instructions, machine-related instructions, microcode, firmware instructions, conditional setting values, or source code or object code written in any combination of one or more programming languages, including object-oriented programming languages ​​such as Sdaletal A, H++, etc., and conventional procedural programming languages ​​such as "H" language or similar programming languages. The computer-readable program instructions can be executed completely on the client computer, partially on the client computer, as a separate software package, partially on the client computer and partially on the remote computer, or completely on the remote computer or server. In the case of a remote computer, the remote computer can be connected to the client computer through any wireless network, including a local area network (LAb) or a wide area network (WAb), or can be connected to an external computer (just like using an Internet service provider to connect through the Internet). In some embodiments, the electronic circuit is customized by using the state values ​​of computer-readable program instructions, such as a programmable logic circuit, a field programmable gate array (FPGA) or a programmable logic array (PLA), and the electronic circuit can execute computer-readable program instructions to achieve each aspect of the overhead disclosure.

[0084] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the above embodiments, ordinary technicians in the relevant field should understand that the specific implementation methods of the present invention can still be modified or replaced by equivalents, and any modifications or equivalent replacements that do not deviate from the spirit and scope of the present invention should be covered within the protection scope of the claims of the present invention.

Claims

1. A method for clearing inter-provincial auxiliary services, characterized in that: include: Construct an inter-provincial electricity spot market clearing model, where the system load balance constraint of the electricity spot market clearing model is based on the output of inter-provincial coal-fired power units in a set period of time, the output of gas-fired power units in a set period of time, the output of non-market units in a set period of time, the power of inter-provincial interconnection lines in a set period of time, and the total load of inter-provincial interconnection lines in a set period of time. The output of inter-provincial coal-fired power units in a set period of time, the output of gas-fired power units in a set period of time, the output of non-market units in a set period of time or the power of inter-provincial interconnection lines in a set period of time is the inter-provincial auxiliary service value; The clearing methods for inter-provincial auxiliary services also include: Step 1: collecting the inter-provincial auxiliary service value of the inter-provincial device, where the inter-provincial auxiliary service value contains multiple subgroups; Step 2: Identify the inter-provincial ancillary service value mutation of each subgroup; Step 3: Calculate the coordination amount of each sampling point within the inter-provincial auxiliary service value; Step 4: Calculate the Z score of each sampling point according to the Z-score method, multiply the coordination amount by the Z score to obtain the final Z score of each sampling point, and perform detection according to the final Z score.

2. The method for clearing inter-provincial auxiliary services according to claim 1, characterized in that: In step 1, the inter-provincial device is an inter-provincial coal-fired power unit, a gas-fired power unit or an inter-provincial interconnection line.

3. The method for clearing inter-provincial auxiliary services according to claim 2, characterized in that: In step 1, in order to sample the inter-provincial auxiliary service value of the inter-provincial device in real time, a power transmitter can be installed on the cable where the inter-provincial device is located, so that when the inter-provincial device is an inter-provincial coal-fired power unit, a gas-fired power unit or an inter-provincial interconnection line, samples can be taken to obtain the output of the inter-provincial coal-fired power unit in a set period of time, the output of the gas-fired power unit in a set period of time, the output of the non-marketized unit in a set period of time or the power of the inter-provincial interconnection line in a set period of time.

4. The method for clearing inter-provincial auxiliary services according to claim 3, characterized in that: In step 1, the sampling speed of the power transmitter is 500 times / s, so as to perform sampling on the inter-provincial auxiliary service value of the inter-provincial device, and then transfer the sampled inter-provincial auxiliary service value to the information table, and then perform service value sorting on the sampled inter-provincial auxiliary service value; In step 1, the service value sorting method is the Wiener filtering algorithm; In step 1, the sampled inter-provincial auxiliary service values ​​are divided into multiple subgroups according to the 300-s time interval formed by the sampling time as a cutting unit.

5. The method for clearing inter-provincial auxiliary services according to claim 4, characterized in that: In step 2, the mutation amount of the inter-provincial auxiliary service value of each subgroup is obtained through the Laplace transform method, that is, the Laplace transform method is used to perform time-frequency conversion on the sampled inter-provincial auxiliary service values, so as to separate the inter-provincial auxiliary service values ​​into phasors of different frequencies, so as to analyze the phasors outside the industrial frequency in the inter-provincial auxiliary service values, and then obtain the amplitude of each phasor outside the industrial frequency and the amplitude of the phasor of the industrial frequency, so as to respectively regard them as the actual amount of each phasor outside the industrial frequency and the actual amount of the phasor of the industrial frequency.

6. The method for clearing inter-provincial auxiliary services according to claim 5, characterized in that: In step 2, the inter-provincial auxiliary service value mutation amount of each subgroup is determined based on the actual amount of each phase quantity outside the power frequency and the actual amount of the phase quantity at the power frequency; The calculation equation of the inter-provincial auxiliary service value mutation is: In the equation, h γ represents the inter-provincial auxiliary service value mutation of the γth subgroup, V i Represents the actual quantity of the phasor outside the industrial frequency with power i, where i = 2, 3, 4..., V1 represents the actual quantity of the phasor of the industrial frequency.

7. The method for clearing inter-provincial auxiliary services according to claim 6, characterized in that: In step 2, the method for obtaining V1 may also be: calculating the average of the power frequency phasors of all sampling points in a subgroup, where the average is the actual value of the power frequency phasor.

8. The method for clearing inter-provincial auxiliary services according to claim 7, characterized in that: In step 3, the calculation equation of the coordination quantity is: In the equation, η γ，j represents the coordination amount of the jth sampling point in the γth subgroup; represents the amount obtained by subtracting the inter-provincial auxiliary service value mutation of the γth subgroup from the inter-provincial auxiliary service value mutation of the other subgroups; ν (j,l) represents the value of the j-th sampling point in the γ-th subgroup minus the value of the l-th sampling point, p represents the total number of additional sampling points in the γ-th subgroup other than the j-th sampling point, and e represents the Euler number.

9. The method for clearing inter-provincial auxiliary services according to claim 8, characterized in that: In step 4, when performing detection based on the final Z score of each sampling point, a stable critical value is pre-defined for each sampling point in the same subgroup. When detecting whether a sampling point is abnormal, the final Z score of the sampling point is compared with the stable critical value corresponding to the sampling point. If the final Z score of the sampling point is greater than its stable critical value, then the sampling point is an abnormal point, otherwise it is a reasonable point. Then, all the abnormal points in the subgroup are totaled. If the total number of abnormal points in a subgroup exceeds half of the total number of sampling points, the maintainer is informed to perform maintenance.

10. A clearing device for inter-provincial auxiliary services, characterized in that: include: A subgroup module, which is used to collect the inter-province auxiliary service value of the inter-province device, and the inter-province auxiliary service value contains multiple subgroups; A recognition module, which is used to recognize the mutation amount of the inter-provincial auxiliary service value of each subgroup; A calculation module, which is used to calculate the coordination amount of each sampling point in the inter-provincial auxiliary service value; The detection module is used to calculate the Z score of each sampling point according to the Z-score method, multiply the coordination amount by the Z score to obtain the final Z score of each sampling point, and perform detection according to the final Z score.

Images