CN112039558B

CN112039558B - Network synchronous clock-based distribution room identification method and device

Info

Publication number: CN112039558B
Application number: CN202010954399.2A
Authority: CN
Inventors: 庞浩
Original assignee: Henan Intelligent Microelectronics Co ltd
Current assignee: Henan Intelligent Microelectronics Co ltd
Priority date: 2020-09-11
Filing date: 2020-09-11
Publication date: 2021-07-09
Anticipated expiration: 2040-09-11
Also published as: CN112039558A

Abstract

The invention relates to a method and a device for identifying a distribution room based on a network synchronous clock, which comprises the following steps: 1) calculating the phase difference of the network, and calculating the voltage phase difference R of the node relative to the main node of the network_{S_C1}(ii) a 2) Calculating the phase difference of the non-network-access node, and calculating the voltage phase difference R of the node relative to the main node of the non-network-access node_{S_C2}(ii) a 3) Comparison of R_{S_C1}And R_{S_C2}Size, if R_{S_C1}If the node is smaller, the node belongs to a power supply area corresponding to the network; if R is_{S_C2}If the node is smaller, the node belongs to a power supply area corresponding to the network which is not accessed; and a processing algorithm of clock synchronization feedback is adopted when the voltage phase difference of the network master node relative to the network master node is calculated. The invention has the beneficial effects that: based on the beacon messages of the accessed network CCO1 and the un-accessed network CCO2', a synchronous clock scale for phase difference measurement is established, so that the phases of the two networks can be accurately compared.

Description

Network synchronous clock-based distribution room identification method and device

Technical Field

The invention relates to the technical field of power detection, in particular to a method and a device for identifying a distribution room based on a network synchronous clock.

Background

In a power distribution network, due to various problems of manual installation, error recording or untimely recording caused by transformation of a transformer area, and the like, the transformer area house change relationship between a carrier electric energy meter and a terminal corresponding to a power supply transformer is wrong.

In order to solve the problem, the phase difference between the electric energy meter of the power supply network and the terminal voltage is tracked by means of phase difference analysis of voltages at different positions of the power distribution network, so that the station-area user variation relationship between the electric energy meter and the terminal corresponding to the power supply transformer is identified.

For example, the invention of the application of vibration table area identification in the power distribution network (patent application number: 201711343618.8) in the invention patent of China (method, system and device for identifying the table area) proposes a method for obtaining a zero-crossing phase by using voltage curve fitting and using the zero-crossing phase for identifying the table area. In order to obtain an accurate zero-crossing phase measurement value, the method adopts the steps of sampling a voltage signal analog-to-digital conversion circuit ADC, storing data, and then performing curve fitting on the data, so that a zero-crossing phase value is determined and is used for station area identification.

This prior art has some problems in implementation: according to the current state of the carrier communication technology, only when a communication module on an ammeter is accessed to a terminal network, a synchronous clock is established with the accessed network, and the synchronous clock is not established with an adjacent network. For an electric energy meter which is already connected to the network, if the network connected to the electric energy meter is suspected to be wrong (the network is in one-to-one correspondence with the station areas, namely, the relationship between the corresponding station areas is suspected to be wrong), the credibility of a calculation result is poor due to the fact that the two networks do not have synchronous clock scales; in addition, various interference factors can not judge the real area relationship of the electric energy meter at all.

In addition, in the prior art, an analog-to-digital conversion circuit is adopted to obtain a voltage curve, and further data fitting and processing are performed, so that the circuit system is complex, has high requirement on sampling frequency and large calculation amount, and is difficult to realize.

Disclosure of Invention

The application aims to provide a method and a device for identifying a transformer area based on a network synchronous clock, which are used for solving the problem of how to verify the subscriber-to-variable relationship of the transformer area of a network-accessed node.

In order to achieve the above object, the present invention provides a method for identifying a distribution room based on a network synchronous clock, comprising:

1) calculating the phase difference of the network, and calculating the voltage phase difference R of the node relative to the main node of the network_{S_C1}；

2) Calculating the phase difference of the non-network-access node, and calculating the voltage phase difference R of the node relative to the main node of the non-network-access node_{S_C2}；

3) Comparison of R_{S_C1}And R_{S_C2}Size, if R_{S_C1}If the node is smaller, the node belongs to a power supply area corresponding to the network; if R is_{S_C2}If the node is smaller, the node belongs to a power supply area corresponding to the network which is not accessed;

wherein, the master node which does not enter the network is set as CCO2', and the calculation of the phase difference of the local node relative to the network which does not enter the network comprises:

according to the beacon message of CCO2', the local synchronous clock count value A is locked_SGet CCO2' Beacon-locked local clock count value A_S2(k) (ii) a Analysis ofThe beacon message of the CCO2 'obtains the beacon clock count value A of the CCO2' transmitted in the message_C2'(k)；

The CCO2' beacon locks the local clock count value a_S2(k) Synchronizing count value B with non-network_S2(k) Making a difference according to A_C2' (k) processing to obtain synchronization error d of non-network₂(k) For said non-network-entry synchronization error d₂(k) Performing proportional integral control to obtain clock error e_S(k) (ii) a Clock error e_S(k) Accumulating to obtain the non-network-access synchronous count value B_S2(k) (ii) a Analyzing and processing the CCO2 'beacon message to obtain a CCO2' beacon clock count value A transmitted in the message_C2' (k); through proportional integral control, the non-network-accessed synchronous count value B is enabled to be_S2(k) And CCO2' beacon clock count value a_C2' (k) synchronization;

computing the CCO2' synchronization zero crossing time stamp y_C2(i)；

Acquiring zero-crossing sequence y of the node_S(i)；

Two sequences y_C2(i) And y_S(i) Taking a plurality of continuous points to do difference, summation, averaging and conversion to obtain the voltage phase difference value R of the node relative to the master node CCO2' not connected into the network_{S_C2}。

Furthermore, the voltage phase difference value R of the local node relative to the main node CCO2' not entering the network_{S_C2}Comprises the following steps:

wherein, T₀Is a power frequency period, Q is a positive integer, f_SThe clock frequency is the working clock frequency of the node.

Further, the synchronization error d of the non-network-access₂(k) Comprises the following steps:

starting from k-2, calculating the non-entry network according to the following formulaSynchronous count value B_S2(k)：

B_S2(k)＝B_S2(k-1)+(A_C2’(k)-A_C2’(k-1))·(1+e_S2(k-1))；

Wherein: b is_S2(0)＝A_S2(0)，B_S2(1)＝A_S2(1)。

Further, CCO2' synchronizes zero crossing time stamps y_C2(i) Comprises the following steps:

y_C2(i)＝B_S2(k)+(y_C2’(i)-A_C2’(k))·(1+e_S(k))+C_C2，

where the i-numbered sequence is different from the k-numbered sequence, B_S2(k)、A_C2'(k)、e_S(k) Get and y_C2' (i) temporally closest data, i ═ l, l +1, …, l + Q-1, C_C2For the delay compensation value of the zero-crossing time scale of the non-accessed network, when the main node CCO2' and the sub-node STA of the non-accessed network are connected to the same standard voltage source, the delay compensation value is set by setting C_C2The voltage phase difference between the voltage measured by the node and the CCO2' is counted by the counter value R_{S_C2}Is 0.

Furthermore, the voltage phase difference value R of the local node relative to the main network access node CCO1_{S_C1}Comprises the following steps:

wherein, T₀Is a power frequency period, Q is a positive integer, f_SFor the operating clock frequency of the local node, y_C1(i) The timestamp count value is zero crossings by CCO 1.

Further, the zero-crossing sequence y of the node_S(i) Calculated according to the following formula:

wherein i is m.N_P，j＝1,2,…,N_P；N_PIs the window length; c_C1For local zero-crossing smoothing metersThe delay compensation value of the numerical value is obtained by setting the delay compensation value C when the main node CCO1 and the sub-node STA are connected with the same standard voltage source_C1So that the voltage phase difference value R measured by STA with CCO1_{S_C1}Is 0;

and: t is_S(m)＝M_S(m)-M_S(m-1)；M_s(m) is N_PLocking the local synchronous clock count value A for zero crossings for window lengths_S1(i) And (6) averaging.

The invention also provides a station area identification based on the network synchronous clock, which comprises a memory, a processor, a communication module for acquiring the network messages and the network messages which are not accessed, a power line voltage signal acquisition module and a computer program stored on the memory; the processor implements the above-described method of station area identification when executing the computer program.

Furthermore, the power line voltage signal acquisition module is a zero-crossing detection and optical coupling isolation circuit.

The invention has the beneficial effects that: based on the beacon messages of the accessed network CCO1 and the un-accessed network CCO2', a synchronous clock scale for phase difference measurement is established, so that the phases of the two networks can be accurately compared. Especially for the non-network, the clock error e between the non-networks is measured by adopting a processing algorithm of clock synchronization feedback_S(k) Establishing local working clock frequency f of sub-node STA_SA base synchronization count.

The invention does not adopt analog-to-digital conversion processing, but realizes the acquisition of the zero-crossing detection signal through the zero-crossing detection and the optical coupling isolation circuit, and filters the noise and the high-frequency interference on the voltage zero-crossing signal by utilizing the filtering function of the circuit. Is easier to realize, and the measuring result is more accurate

The smoothing processing of the local zero-crossing count value, the feedback processing of the network clock without entering and the summation and averaging of the voltage phase difference further filter the interference of noise and jitter on the phase difference measurement result in a software mode.

The invention uses the delay compensation value C of the local zero-crossing smooth counting value_C1And non-network-entry zero crossingTarget delay compensation value C_C2The solution CCO1 or CCO2' adds its own clock count value A to the beacon message_C1(k) Or A_C2'k's processing delay, communication transmission delay, the different modes of sub-node STA recording its own clock, network access synchronization error etc. influence the problem of phase difference measurement, based on checking with the voltage, measure and calibrate, have guaranteed the accurate measurement of sub-node STA voltage phase difference between its network access and not network access.

Drawings

FIG. 1 is a zero crossing detection and opto-isolator circuit according to an embodiment of the present invention;

FIG. 2 is a schematic diagram of voltage phase difference calculation according to an embodiment of the present invention;

FIG. 3 is a diagram of zero crossing locked local synchronous clock count value A in an embodiment of the present invention_S1(i) And a local zero crossing smoothing count value y_S(i) An error test curve;

FIG. 4 is a flowchart illustrating an exemplary non-network-entry clock synchronization feedback process according to an embodiment of the present invention;

FIG. 5 is a diagram illustrating a clock error e of a non-network-access clock synchronization feedback process according to an embodiment of the present invention_S(k) And (6) testing the curve.

Detailed Description

The application scenario of the embodiment is an electricity utilization acquisition system, which comprises a plurality of high-speed carrier communication modules, wherein the carrier communication modules are arranged in terminal equipment (such as an electric energy meter, an acquisition unit and the like). The embodiment relates to 3 modules, including a sub-node STA, two main nodes CCO1 and CCO 2'; the master node CCO1 belongs to the a network and the other master node CCO2' belongs to the B network. The network A and the network B respectively correspond to the power supply area a and the power supply area B.

The child node STA has already joined the a network, so that the network a to which the CCO1 belongs is called "a visited network", and the network B to which the CCO2' belongs is called "a missed network".

Since the sub-node STA has joined the "network entered", the station area to which the sub-node STA belongs is the power supply station area a. However, due to some of the problems noted in the background, the child node STA may not really belong to the power supply station area a. Therefore, the current problem is to identify the subscriber relationship of the child node STA.

The basic idea of the embodiment is as follows: calculating voltage phase difference R between the sub-node STA and the CCO1_{S_C1}(ii) a And calculating the voltage phase difference R between the sub-node STA and the CCO2_{S_C2}. Then compare R_{S_C1}And R_{S_C2}The size of (d); the smaller the phase difference, the closer the relationship between the variables, so the smaller one corresponds to the true relationship between the variables. For example, if R_{S_C1}If the STA is smaller, the sub-node STA belongs to the CCO1 network, namely belongs to the power supply area a; if R is_{S_C2}And smaller, the sub-node STA belongs to the CCO2' network, namely, belongs to the power supply station area b.

In this example, R_{S_C1}And R_{S_C2}Refers to the phase difference obtained from a single measurement.

As another embodiment, considering that an abnormal condition or various interference conditions may occur in a single measurement, generally, the determination is not directly performed with a single measurement result, but needs a period of time, in which multiple measurements are performed, and the comprehensive determination is performed according to the multiple measurement results:

for example, if more than 85% of the measurements are less voltage-different from the CCO1 over a period of time, the child node STA belongs to the CCO1 network.

For another example, if 50 measurements are performed over a period of time, the phase differences obtained from the 50 measurements are summed (with the CCO1 voltage phase difference and with the CCO2' voltage phase difference), and the final determination is made based on which of the summed total phase differences is smaller.

That is, as another embodiment, R_{S_C1}And R_{S_C2}Or may be understood as a comprehensive judgment result of a plurality of measurements over a period of time.

The embodiment adopts a zero-crossing detection and optical coupling isolation circuit to replace an ADC circuit in the prior art. The ZERO-crossing detection and optical coupling isolation input circuit is shown in fig. 1, the circuits AC _ LA and AC _ N are respectively connected with a live wire and a ZERO wire of a 220V power line, the conduction of a triode D20 and a light emitting diode in an optical coupling D24 is triggered at the ZERO-crossing moment of a voltage falling edge, and a high-level ZERO-crossing detection signal is output at the ZERO _ T output end every 20ms according to the falling edge of a power frequency period. C51 charge is discharged, then D20 is turned off, the optocoupler D24 is turned off, and then ZERO _ T goes back low. The circuit can also perform low-frequency filtering on transient noise on the voltage signal through the capacitor C11, the resistor R6 and the capacitor C52.

Compared with an ADC circuit, the zero detection and optical coupling isolation circuit is relatively simple and easy to realize, and the filtering function of the circuit filters noise and high-frequency interference on the voltage zero-crossing signal.

The specific voltage phase difference calculation process is described in detail below. As shown in fig. 2, the method mainly comprises two steps, which are performed simultaneously: firstly, calculating the network phase difference, namely calculating the voltage phase difference R between the output sub-node STA and the CCO1_{S_C1}. Secondly, calculating the voltage phase difference R between the sub-node STA and the CCO2' without network phase difference calculation_{S_C2}。

The method comprises the following steps:

according to the technical specification of interconnection and intercommunication of low-voltage power line broadband carrier communication of national grid company enterprise standard Q/GDW 11612.3-2016, in a low-voltage power line broadband carrier communication network, a main node CCO1 sends a beacon message containing a network reference clock mark, and a sub-node STA in the whole network needs to keep synchronous with the network reference time of the main node CCO 1. The network reference clock is a CCO1 hardware clock f_SCount frequency 25MHz, increment by 1 with 32-bit unsigned number. In a low-voltage power line broadband carrier communication network, a sub-node STA realizes synchronous counting of a network-in local network according to a network-in main node CCO1 of the sub-node STA to obtain a frequency signal f according to a working clock_S32-bit unsigned local synchronous clock count value a at 25MHz synchronized with CCO1 operating clock by synchronization adjustment_S。

The method comprises the steps that a sub-node STA obtains a beacon message and a zero-crossing time mark message of a network access main node CCO1 through a communication network, and voltage phase difference R between the sub-node STA and the network access main node CCO1 is achieved according to the following steps_{S_C1}Measurement:

step 1.1, the sub-node STA realizes synchronous counting processing of the network-access local network according to the unsigned cycle clock count value of the CCO1 contained in the beacon message of the network-access main node CCO1, and obtains the working clock frequency according to the STA asf_SLocal synchronous clock count value A adjusted to be synchronous with CCO1 working clock through synchronization_S。

Step 1.2, the power line voltage signal obtains a zero-crossing detection signal through a zero-crossing detection and optical coupling isolation input circuit, and the local synchronous clock is counted A through triggering of the zero-crossing detection signal_SObtaining the local synchronous clock counting value A of the zero-crossing locking_S1(i) (ii) a I.e. by locking the local synchronous clock count value a_SGeneration of A_S1(i)。

Step 1.3, taking the number N of smooth local zero-crossing count values_PLocking the local synchronous clock count value a for zero crossing at 32_S1(i) Smoothing the local zero-crossing count value when i is m.N_P＝N_P，2N_P，3N_P…, when m is a positive integer, calculate the mean sequence:

j is 0 to N _p1, with N_pThe window length is averaged.

Wherein N is_PTaking the value as a positive integer;

T_S(m)＝M_S(m)-M_S(m-1)；

local zero crossing smoothing count value y_S(i) Calculated according to the following formula:

wherein i is m.N_P，j＝1,2,…,N_P，C_C1For the delay compensation value of the local zero-crossing smooth count value, when the main node CCO1 and the sub-node STA are connected with the same standard voltage source, the set delay compensation value C is used_C1So that the voltage phase difference value R measured by STA with CCO1_{S_C1}Is 0.

Local zero crossing smoothing count value y_S(i) Equivalent to locking the local synchronous clock count value A to the zero crossing_S1(i) Go on heavilyAnd (6) sampling.

According to step 1.3, if the local synchronous clock count value A of the sub-node STA_SCompared with the stable clock count value, the random jitter with the maximum amplitude of 100 exists, and the random error of the zero-crossing detection signal output by the zero-crossing detection and optical coupling isolation input circuit is 5us at most, so that the local synchronous clock count value A is locked by zero crossing_S1(i) Smoothing the count value y with local zero crossings_S(i) The test curve plotted against the deviation from the stable clock timing is shown in fig. 3. As can be seen from FIG. 3, the local zero crossing smoothed count value y_S(i) Comparison with A_S1(i) Noise interference is filtered, and a more stable and accurate zero-crossing counting value is obtained, so that the phase difference can be accurately measured.

Step 1.4, the sub-node STA obtains a partial sequence y of CCO1 zero-crossing time mark count values according to the received CCO1 zero-crossing time mark message_C1(i) I.e. i-l, l +1, …, l + Q-1, corresponds to the zero-crossing sequence y of the child node STA taking the point Q-16_S(i) Then, the voltage phase difference value R of the sub-node STA with respect to the network master node CCO1 is calculated according to the following formula_{S_C1}；

Wherein, T₀Is the power frequency cycle.

The above formula is equivalent to two sequences y_C1(i) And y_S(i) And carrying out difference, summation, averaging and conversion to phase at a plurality of continuous points.

With the above method, a phase difference with extremely high accuracy can be obtained. Since the phase difference calculation is easy to implement because the phase difference is already present in the network, and there are other methods for calculating the phase difference in the prior art, the measurement can be implemented in other methods in the prior art as well as in other embodiments.

Step two:

measuring voltage phase difference value R of quantum node STA relative to network-accessed main node CCO1_{S_C1}Meanwhile, the sub-node STA obtains a beacon message and a zero-crossing time mark message of the master node CCO2' which do not enter the network through the communication network according to the methodThe next step realizes the voltage phase difference R between the sub-node STA and the main node CCO2' which does not enter the network_{S_C2}Measurement:

step 2.1, the communication processing part of the sub-node STA receives the beacon message of the main node CCO2', immediately generates an effective CCO2' beacon receiving mark signal, and the CCO2' beacon receiving mark signal locks the local synchronous clock count value A of the sub-node STA_SGet CCO2' Beacon-locked local clock count value A_S2(k) (ii) a Meanwhile, the beacon message of the master node CCO2 'is analyzed to obtain the beacon clock count value A of CCO2' transmitted in the message_C2'(k)；

Step 2.2, the sub-node STA adopts the processing flow of synchronous feedback of the network clock which is not accessed and is shown in fig. 4:

(1) beacon locking local clock count value A according to CCO2_S2(k) And the non-network-access synchronous count value B_S2(k) Obtaining the synchronization error d of the non-network according to the following formula₂(k)：

Wherein A is_C2' (k) is a CCO2' beacon clock count value transmitted in a message obtained by analyzing and processing a CCO2' beacon message;

(2) with synchronization error d of non-network₂(k) For input, proportional integral control calculation is carried out to obtain a clock error e_S(k) (ii) a Proportional clock error e in proportional integral control_P(k) Calculated according to the following formula:

wherein, K_PIs a proportional control coefficient, set according to the dynamic performance of the non-network-connected clock synchronous feedback system, e_PmaxAnd e_PminThe upper limit threshold and the lower limit threshold of the proportional clock error are respectively set according to the clock requirement of the communication equipment, the calculation of the above formula is started from k to 2, e_P(k) Initial value set to e_P(0)＝e_P(1) 0; integral clock error e in proportional-integral control_I(k) Calculated according to the following formula:

wherein, K_IIs a proportional control coefficient, set according to the dynamic performance of the non-network-connected clock synchronous feedback system, e_ImaxAnd e_IminThe upper limit threshold and the lower limit threshold of the proportional clock error are respectively set according to the clock requirement of the communication equipment, the calculation of the above formula is started from k to 2, e_I(k) Initial value set to e_I(0)＝0，

Finally, the clock error e of the proportional-integral control output_S(k) Calculated according to the following formula:

e_S(k)＝e_P(k)+e_I(k)；

(3) starting from k to 2, the non-network-entry synchronization count value B is calculated according to the following formula_S2(k) And forming feedback:

B_S2(k)＝B_S2(k-1)+(A_C2’(k)-A_C2’(k-1))·(1+e_S2(k-1))

non-network-access synchronous count value B_S2(k) Is set to be: b is_S2(0)＝A_S2(0)，B_S2(1)＝A_S2(1)。

Setting e_Imax＝-e_Imax＝2e-5，ePmax＝-ePmax＝1e-5，K_P＝0.105，K_I0.005625. When the error between the network clock of the network main node CCO1 accessed by the sub-node STA and the network clock of the network main node CCO2' not accessed is 5ppm, the local synchronous clock count value A of the sub-node STA_SCompared with the random jitter with the maximum amplitude of 50 in the stable clock counting value, the beacon sending period is set to be 1.5s, and in order to check the effect of the clock synchronization feedback processing algorithm without network access, the time of the scale is adjustedInitial value of clock error is set as e_P(0)＝e_P(1) 0, the initial value of the integral clock error is set as e_I(0)＝e_I(1) 0. Clock error e of proportional-integral control output_S(k) A test curve tracking the actual clock error by 5ppm is shown in fig. 5. With, clock error e_S(k) Converged, non-network-entered synchronization count value B_S2(k) Will also implement the beacon clock count value A with CCO2_C2' (k) synchronization.

Step 2.3, the child node STA receives the CCO2' zero-crossing time mark y of the CCO2' in the CCO2' zero-crossing time mark message_C2' (i) synchronizing count value B according to non-network_S2(k) CCO2' beacon clock count value a_C2' (k) and clock error e_S(k) And performing zero-crossing time scale adjustment processing according to the following formula to obtain a CCO2' synchronization zero-crossing time scale:

y_C2(i)＝B_S2(k)+(y_C2’(i)-A_C2’(k))·(1+e_S(k))+C_C2，

where the i-numbered sequence is different from the k-numbered sequence, B_S2(k)、A_C2'(k)、e_S(k) Get and y_C2' (i) temporally closest data, i ═ l, l +1, …, l + Q-1, C_C2For the delay compensation value of the zero-crossing time scale of the non-accessed network, when the main node CCO2' and the sub-node STA of the non-accessed network are connected with the same standard voltage source, the delay compensation value C of the zero-crossing time scale of the non-accessed network is set_C2So that the voltage phase difference between STA measured and CCO2' is counted as R_{S_C2}Is 0;

step 2.4, zero-crossing sequence y of sub-node STA_S(i) L, l +1, …, l + Q-1, and calculating the voltage phase difference value R of the sub-node STA relative to the non-network-entry main node CCO2' according to the following formula_{S_C2}：

Wherein, T₀Is the power frequency cycle. The above formula is equivalent to two sequences y_C2(i) And y_S(i) In a plurality of continuous pointsDifference, sum, average, and conversion to phase.

The scheme of the embodiment is applied to the phase difference measurement between the electric energy meter of the low-voltage power supply network and the terminal voltage, the accurate measurement result of 0.002 degree can be obtained, and the measurement precision is high.

After the first step and the second step are completed, two phase differences R are subjected to_{S_C1}、R_{S_C2}And comparing to determine the network and power supply area to which the subnode really belongs.

The above is the introduction of the method of the present invention, a station area identification method based on a network synchronous clock. The invention also relates to an embodiment of the device for realizing the method, which comprises the following steps:

a station area recognition device based on a network synchronous clock comprises a memory, a processor, a communication module for collecting network messages and non-network messages, and a power line voltage signal collection module (in the embodiment, a zero-crossing detection and optical coupling isolation circuit is adopted); a computer program stored on the memory; the processor implements the above-described method of station area identification when executing the computer program. The processor can be a computer, a microprocessor such as an ARM and the like, and a programmable chip such as an FPGA, a DSP and the like.

As will be appreciated by one skilled in the art, 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.

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. 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.

Finally, it should be noted that: although the present invention has been described in detail with reference to the above embodiments, it should be understood by those skilled in the art that: modifications and equivalents may be made to the embodiments of the invention without departing from the spirit and scope of the invention, which is to be covered by the claims.

Claims

1. A method for identifying a distribution area based on a network synchronous clock is characterized by comprising the following steps:

according to the beacon message of CCO2', the local synchronous clock count value A is locked_SGet CCO2' Beacon-locked local clock count value A_S2(k) (ii) a Analyzing the beacon message of the CCO2 'to obtain the beacon clock count value A of the CCO2' transmitted in the message_C2'(k)；

The CCO2' beacon locks the local clock count value a_S2(k) Synchronizing count value B with non-network_S2(k) Making a difference according to A_C2' (k) processing to obtain synchronization without network accessError d₂(k) For said non-network-entry synchronization error d₂(k) Performing proportional integral control to obtain clock error e_S(k) (ii) a Clock error e_S(k) Accumulating to obtain the non-network-access synchronous count value B_S2(k) (ii) a Analyzing and processing the CCO2 'beacon message to obtain a CCO2' beacon clock count value A transmitted in the message_C2' (k); through proportional integral control, the non-network-accessed synchronous count value B is enabled to be_S2(k) And CCO2' beacon clock count value a_C2' (k) synchronization;

computing the CCO2' synchronization zero crossing time stamp y_C2(i)；

Acquiring zero-crossing sequence y of the node_S(i)；

2. The network-synchronized-clock-based station area identification method of claim 1, wherein the voltage phase difference value R of the local node relative to the non-network-connected master node CCO2_{S_C2}Comprises the following steps:

3. The network-synchronized clock-based station area identification method of claim 1, wherein the non-network-synchronized error d₂(k) Comprises the following steps:

starting from k to 2, the non-network-entry synchronization count value B is calculated according to the following formula_S2(k)：

B_S2(k)＝B_S2(k-1)+(A_C2’(k)-A_C2’(k-1))·(1+e_S2(k-1))；

Wherein: b is_S2(0)＝A_S2(0)，B_S2(1)＝A_S2(1)。

4. The network-synchronized clock-based station area identification method of claim 3, wherein the CCO2' synchronizes the zero-crossing time stamp y_C2(i) Comprises the following steps:

y_C2(i)＝B_S2(k)+(y_C2’(i)-A_C2’(k))·(1+e_S(k))+C_C2，

5. The network synchronized clock-based distribution area identification method of claim 1, wherein the voltage phase difference value R of the local node relative to the network access master node CCO1 is_{S_C1}Comprises the following steps:

6. The network synchronous clock-based station area identification method according to claim 1, wherein the zero-crossing sequence y of the node_S(i) AnCalculated according to the following formula:

wherein i is m.N_P，j＝1,2,…,N_P；N_PIs the window length; c_C1For the delay compensation value of the local zero-crossing smooth count value, when the main node CCO1 and the sub-node STA are connected with the same standard voltage source, the set delay compensation value C is used_C1So that the voltage phase difference value R measured by STA with CCO1_{S_C1}Is 0;

7. A station area identification device based on a network synchronous clock is characterized by comprising a memory, a processor, a communication module for acquiring a network message and a network message which are not accessed, a power line voltage signal acquisition module and a computer program stored on the memory; the processor, when executing the computer program, implements the method of any of claims 1-6.

8. The network synchronous clock-based station area identification device according to claim 7, wherein the power line voltage signal acquisition module is a zero-crossing detection and optical coupling isolation circuit.