CN113779751A - Low-pressure HPLC (high performance liquid chromatography) platform area topology identification method and system - Google Patents

Low-pressure HPLC (high performance liquid chromatography) platform area topology identification method and system Download PDF

Info

Publication number
CN113779751A
CN113779751A CN202110853927.XA CN202110853927A CN113779751A CN 113779751 A CN113779751 A CN 113779751A CN 202110853927 A CN202110853927 A CN 202110853927A CN 113779751 A CN113779751 A CN 113779751A
Authority
CN
China
Prior art keywords
length
module
branch
frequency
calculating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202110853927.XA
Other languages
Chinese (zh)
Inventor
杨挺
冯少康
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN202110853927.XA priority Critical patent/CN113779751A/en
Publication of CN113779751A publication Critical patent/CN113779751A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J13/00Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network
    • H02J13/00006Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network characterised by information or instructions transport means between the monitoring, controlling or managing units and monitored, controlled or operated power network element or electrical equipment
    • H02J13/00007Circuit arrangements for providing remote indication of network conditions, e.g. an instantaneous record of the open or closed condition of each circuitbreaker in the network; Circuit arrangements for providing remote control of switching means in a power distribution network, e.g. switching in and out of current consumers by using a pulse code signal carried by the network characterised by information or instructions transport means between the monitoring, controlling or managing units and monitored, controlled or operated power network element or electrical equipment using the power network as support for the transmission
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/04Power grid distribution networks
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S40/00Systems for electrical power generation, transmission, distribution or end-user application management characterised by the use of communication or information technologies, or communication or information technology specific aspects supporting them
    • Y04S40/12Systems for electrical power generation, transmission, distribution or end-user application management characterised by the use of communication or information technologies, or communication or information technology specific aspects supporting them characterised by data transport means between the monitoring, controlling or managing units and monitored, controlled or operated electrical equipment
    • Y04S40/121Systems for electrical power generation, transmission, distribution or end-user application management characterised by the use of communication or information technologies, or communication or information technology specific aspects supporting them characterised by data transport means between the monitoring, controlling or managing units and monitored, controlled or operated electrical equipment using the power network as support for the transmission

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

The invention relates to a method and a system for identifying topology of a low-pressure HPLC (high performance liquid chromatography) platform area, which comprises the following steps: step 1, numbering all electric meters according to the distribution area files; step 2, sampling to obtain a frequency response sequence; step 3, performing wavelet denoising on the frequency response sequence; step 4, determining the number of the transmitting and receiving intermediate branches and the branch length by utilizing the characteristic of the trapped wave frequency point; step 5, calculating the length of the bus between the transmitting and receiving nodes by using the average channel attenuation characteristic; step 6, calculating the distance between the carrier wave electric meters in the middle, and adding and deleting middle branch nodes; and 7, identifying the topology of the transformer area. The invention can realize automatic measurement of the line length and accurate and rapid identification of the topological structure of the transformer area without installing an identification device on the intermediate node.

Description

Low-pressure HPLC (high performance liquid chromatography) platform area topology identification method and system
Technical Field
The invention belongs to the technical field of topology identification and reconstruction, and relates to a topology identification method and a system, in particular to a low-pressure HPLC (high performance liquid chromatography) platform area topology identification method and a system.
Background
The low voltage stations are at the end of the power system and take on the task of supplying power directly to the users. In recent years, with the development of low-voltage power distribution networks, the informatization and intelligentization levels of transformer areas are greatly improved. However, in the process of constructing the low-voltage distribution network, the distribution lines and the power facilities are often adjusted and changed due to the increase of the electrical loads, and the power grid company cannot adjust the distribution area files in time due to management, technology and other reasons, so that the topology files of the distribution area are inconsistent with the actual topology files. The wrong topology files can influence the development of line loss treatment and electricity larceny prevention troubleshooting work in the transformer area, and great inconvenience is brought to the fine management work of the transformer area. Therefore, the research on the low-voltage transformer area topology identification technology is of great significance.
The traditional station area identification method mainly comprises two types, one type is a communication-based method, the method needs special station area identification equipment such as a station area identification instrument, or integrates a station area topology identification special communication device at intermediate nodes such as a JP cabinet, a branch box, a light power cabinet and the like, the principle of the method is simple, the implementation difficulty is low, but the equipment expenditure is increased, the other type is a large-data-based method, correlation analysis is carried out on data such as power consumption, voltage and the like of a user, although extra equipment does not need to be added, a large amount of user power consumption data is needed, the requirements on data quality and algorithm performance are high, online identification of network topology cannot be achieved, and the practical application difficulty is high.
In view of this, the present invention provides a method and a system for identifying a topology of a low-pressure HPLC platform.
Upon search, no prior art documents that are the same or similar to the present invention have been found.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a method and a system for identifying the topology of a low-pressure HPLC (high performance liquid chromatography) platform area, which can realize automatic measurement of line length and accurate and rapid identification of the topology structure of the platform area under the condition that an identification device is not required to be installed on an intermediate node.
The invention solves the practical problem by adopting the following technical scheme:
a low-pressure HPLC (high performance liquid chromatography) platform area topology identification method comprises the following steps:
step 1, numbering all electric meters according to the distribution area files, drawing in an unidentified electric meter set B1, and initializing a transition electric meter set
Figure BDA0003183392770000022
Identified electricity meter collection
Figure BDA0003183392770000023
Step 2, in a frequency band of 2-30MHz, respectively transmitting and receiving HPLC sweep frequency signals by a carrier electricity meter I in a set B1 and a station area concentrator, sampling to obtain a frequency response sequence H (f), and dividing all electricity meters receiving the carrier signals in the electricity meter I and the set B1 into a set B2; (ii) a
Step 3, performing wavelet denoising on the frequency response sequence;
step 4, determining the number of the transmitting and receiving intermediate branches and the branch length by utilizing the characteristic of the trapped wave frequency point;
step 5, calculating the length of the bus between the transmitting and receiving nodes by using the average channel attenuation characteristic;
step 6, calculating the distance between the carrier wave electric meters in the middle, and adding and deleting middle branch nodes;
step 7, identifying the topology of the transformer area;
the specific method of step 3 is:
performing wavelet denoising on one-dimensional waveform data obtained by sampling, filtering high-frequency noise, wherein the denoising method selects Daubechies wavelet as wavelet basis, a heursure mixed threshold rule and a soft threshold function, and the threshold calculation rule is as follows:
Figure BDA0003183392770000021
if Eta < Crit, then a fixed threshold is chosen, otherwise the smaller of the rigrsure criterion and the sqtwolog criterion is chosen as the present criterion threshold.
Further, the specific steps of step 4 include:
(1) fast Fourier transform extraction of all significant frequency components F in amplitude spectra1~FnJudging the number of branches to be n;
(2) calculating the notch frequency point interval in the signal frequency response curve corresponding to each branch
Figure BDA0003183392770000031
(3) Calculating the length of each branch:
Figure BDA0003183392770000032
in the formula, L is the distributed inductance of the line, and C is the distributed capacitance of the line.
Further, the specific steps of step 5 include:
(1) firstly, according to the Hankel matrix form, the channel frequency response sequence x is constructed into an m multiplied by n signal matrix HxThe matrix structure is as follows:
Figure BDA0003183392770000033
to HxSingular value decomposition is carried out to obtain sub-channel frequency response H1(f)~Hn(f) The singular value decomposition is in the form:
Hx=USVT (4)
(2) windowed DeltafiAnd performing PE operation to extract the peak points of the frequency division rate response curves, wherein the PE operation is as follows:
Figure BDA0003183392770000035
subtracting the structural element Bm from the original input signal Fn to carry out open operation on the signal Fn, thus obtaining a peak point in each window, fitting a peak curve based on a least square method, and approximately replacing the frequency response of a non-branch channel by the curve;
(3) calculating the average channel attenuation statistical characteristic of the non-branch channel, wherein the calculation formula is as follows:
Figure BDA0003183392770000034
wherein f is0For frequency resolution, take f025kHz, frequency range 2-30MHz, so n1=80,n2=1200。
(4) Calculating the length l of the corresponding main line by ACA (l), wherein the ACA (l) is obtained by outdoor actual measurement, and the measurement rule is as follows:
changing the length l of the non-branch line for many times, and respectively calculating ACA under each l to obtain parameters a and b in the formula (7);
ACA(l)=-al+b(dB) (7)
(5) calculating the main line length l corresponding to each frequency division rate response by the ACA (l) of the formula (3)1~lnAnd the bus length l is obtained by taking the average value.
Further, the specific steps of step 6 include:
(1) collecting the carrier electricity meters i and j in the B2 to receive and transmit HPLC sweep frequency signals, and sampling to obtain Hij(f);
(2) Calculating the line length l between i and j according to step 5ij
(3) When l isij-(lbri+lbrj) When the distance is larger than epsilon, the two electric meters are judged to belong to different branch nodes respectively, and the distance between the branch nodes is lij-(lbri+lbrj);
(4) When lij-(lbri+lbrj) When | ≦ epsilon, judging that the branch nodes to which the two electric meters respectively belong are closer, and merging the branch nodes;
(5) when l isij-(lbri+lbrj) When < -epsilon, two electric meters belong to the same branch, and the branch length is (l)bri+lbrj-lij) Adding middle branch nodes;
(6) and (6) circularly executing the step until the line length calculation among all the electric meters in the B2 is completed.
Further, the specific steps of step 7 include:
(1) calculating the distances between all the electric meters in the B2 and the electric meter I, judging the upstream and downstream relations of the nodes according to the distance, generating a cluster topology, dividing all the elements in the B2 into the B3, and enabling all the elements in the B2 to be in a row
Figure BDA0003183392770000041
The topology identification of the cluster where the ammeter I is located is completed;
(2) when in use
Figure BDA0003183392770000042
And (4) returning to the step (2), otherwise, ending the circulation, and finishing the identification of the platform area topology.
A low pressure HPLC zone topology identification system, comprising: the device comprises a signal acquisition module, a signal processing module, a length calculation module and a topology identification module;
the signal acquisition module inserts the broadband carrier communication module of ammeter and concentrator into this system, includes: the electric meter carrier module and the concentrator carrier module are used for acquiring HPLC signal frequency response one-dimensional waveform data of a receiving end; the data storage module is used for storing the integrated data and uploading the integrated data to the signal processing module;
the signal processing module comprises a first statistical analysis module, a system database and a first data management module; the first statistical analysis module is used for performing wavelet denoising on the signal frequency response acquired by the signal acquisition module, filtering high-frequency noise, performing frequency spectrum analysis and extracting peak point and trapped wave frequency point characteristics; the system database and the first data management module are used for storing data and calling the length calculation module;
the length calculation module comprises a second statistical analysis module and a second data management module; the second statistical analysis submodule is used for calculating branch lengths based on the signal frequency selective attenuation characteristics, calculating line lengths based on the average channel attenuation characteristics, and managing the calculated length data by the second data management module and calling the calculated length data by the topology identification module.
The topology identification module consists of a digital signal processor, a power supply unit, a graphic processor, an MCU unit and a GUI, and is used for calling the length calculation module, determining the number of branch nodes and the upstream-downstream relation, generating a distribution room topological graph and outputting the topological graph.
The invention has the advantages and beneficial effects that:
1. the invention fully considers the frequency selective attenuation characteristic and the mathematical relation between the statistical characteristic and the line length of the broadband power line carrier channel, designs a low-voltage HPLC (high performance liquid chromatography) station topology identification method and system based on the broadband power line carrier characteristic, and finally identifies the topology of a station network by the efficient and creative accurate extraction and calculation method of the line branch number and the line length as described in the steps 4 and 5 and the extraction method of the branch node number and the position as described in the step 6. The invention realizes automatic measurement of line length, accurately and rapidly identifies the topological structure of the transformer area without installing an identification device on an intermediate node, and ensures the normal development of the line loss management and electricity larceny prevention troubleshooting work of the transformer area.
2. The invention discloses a topology identification method suitable for a low-voltage HPLC (high performance liquid chromatography) platform area. And in the later stage, the topological structure between the receiving and transmitting nodes is identified by repeatedly adding and deleting intermediate nodes and determining the upstream and downstream relations of the nodes, and finally, the cluster topology is generated, so that the accurate identification of the whole distribution area topology is realized.
Drawings
FIG. 1 is a process flow diagram of the method of the present invention;
FIG. 2(a) is a schematic diagram illustrating a node deletion rule according to the present invention;
FIG. 2(b) is a schematic diagram of node deletion rule of the present invention;
FIG. 2(c) is a schematic diagram of node deletion rule of the present invention (III);
FIG. 3 is a diagram of a topology identification system for a low pressure HPLC station of the present invention;
FIG. 4 is a topological diagram of an experimental test single cluster block area of the present invention;
FIG. 5 is a graph of the signal frequency response of the present invention;
FIG. 6 is an amplitude spectrum of the present invention;
fig. 7 is a table area topology identification experiment test result diagram of the present invention.
Detailed Description
The embodiments of the invention will be described in further detail below with reference to the accompanying drawings:
a method for identifying topology of low-pressure HPLC (high-pressure liquid chromatography) platform area, as shown in figure 1, comprises the following steps:
step 1, numbering all electric meters according to the distribution area files, drawing in an unidentified electric meter set B1, and initializing a transition electric meter set
Figure BDA0003183392770000061
Identified electricity meter collection
Figure BDA0003183392770000062
Step 2, in a frequency band of 2-30MHz, respectively transmitting and receiving HPLC sweep frequency signals by a carrier electricity meter I in a set B1 and a station area concentrator, sampling to obtain a frequency response sequence H (f), and dividing all electricity meters receiving the carrier signals in the electricity meter I and the set B1 into a set B2; (ii) a
Step 3, performing wavelet denoising on the frequency response sequence;
the specific method of the step 3 comprises the following steps:
wavelet denoising is carried out on the one-dimensional waveform data obtained by sampling, high-frequency noise is filtered, and the denoising method selects Daubechies wavelets as wavelet bases, a heursure mixed threshold rule and a soft threshold function. The threshold calculation rule is as follows:
Figure BDA0003183392770000071
if Eta < Crit, then a fixed threshold is chosen, otherwise the smaller of the rigrsure criterion and the sqtwolog criterion is chosen as the present criterion threshold.
Step 4, determining the number of the transmitting and receiving intermediate branches and the branch length by utilizing the characteristic of the trapped wave frequency point;
the specific steps of the step 4 comprise:
(1) fast Fourier transform extraction of all significant frequency components F in amplitude spectra1~FnJudging the number of branches to be n;
(2) calculating the notch frequency point interval in the signal frequency response curve corresponding to each branch
Figure BDA0003183392770000072
(3) Calculating the length of each branch:
Figure BDA0003183392770000073
in the formula, L is the distributed inductance of the line, and C is the distributed capacitance of the line.
Step 5, calculating the length of the bus between the transmitting and receiving nodes by using the average channel attenuation characteristic;
the specific steps of the step 5 comprise:
(1) firstly, according to the Hankel matrix form, the channel frequency response sequence x is constructed into an m multiplied by n signal matrix HxThe matrix structure is as follows:
Figure BDA0003183392770000074
to HxSingular value decomposition is carried out to obtain sub-channel frequency response H1(f)~Hn(f) The singular value decomposition is in the form:
Hx=USVT (4)
(2) windowed DeltafiAnd performing PE operation to extract the peak points of the frequency division rate response curves, wherein the PE operation is as follows:
Figure BDA0003183392770000082
the result of opening operation is carried out on the signal F [ n ] by subtracting the structural element Bm from the original input signal F [ n ], so that the peak point in each window can be obtained, a peak curve is fitted based on the least square method, and the curve is used for approximately replacing the frequency response of the non-branch channel.
(3) Calculating the average channel attenuation statistical characteristic of the non-branch channel, wherein the calculation formula is as follows:
Figure BDA0003183392770000081
wherein f is0For frequency resolution, take f025kHz, frequency range 2-30MHz, so n1=80,n2=1200。
(4) Calculating the length l of the corresponding main line by ACA (l), wherein the ACA (l) is obtained by outdoor actual measurement, and the measurement rule is as follows:
changing the length l of the non-branch line for many times, and respectively calculating ACA under each l to obtain parameters a and b in the formula (7);
ACA(l)=-al+b(dB) (7)
(5) calculating the main line length l corresponding to each frequency division rate response by the ACA (l) of the formula (3)1~lnAnd averaging to obtain the bus length l;
step 6, calculating the distance between the carrier wave electric meters in the middle, and adding and deleting middle branch nodes;
and (5) solving the line length among the electric meters in the set B2 based on the line length solving principle in the step 5.
The specific steps of the step 6 comprise:
(1) collecting the carrier electricity meters i and j in the B2 to receive and transmit HPLC sweep frequency signals, and sampling to obtain Hij(f);
(2) Calculating the line length l between i and j according to step 5ij
(3) When l isij-(lbri+lbrj) When the distance is larger than epsilon, the two electric meters are judged to belong to different branch nodes respectively, and the distance between the branch nodes is lij-(lbri+lbrj) As shown in FIG. 2 (a);
(4) when lij-(lbri+lbrj) When | ≦ ε, it is determined that the branch nodes to which the two ammeters belong are closer, and the branch nodes are merged as shown in FIG. 2 (b);
(5) when l isij-(lbri+lbrj) When < -epsilon, two electric meters belong to the same branch, and the branch length is (l)bri+lbrj-lij) 2, adding an intermediate branch node, as shown in FIG. 2 (c);
(6) circularly executing the step 6 until the line length calculation among all the electric meters in the B2 is completed;
step 7, identifying the topology of the transformer area
The specific steps of the step 7 comprise:
(1) calculating the distances between all the electric meters in the B2 and the electric meter I, judging the upstream and downstream relations of the nodes according to the distance, generating a cluster topology, dividing all the elements in the B2 into the B3, and enabling all the elements in the B2 to be in a row
Figure BDA0003183392770000092
The topology identification of the cluster where the ammeter I is located is completed;
(2) when in use
Figure BDA0003183392770000091
And (4) returning to the step (2), otherwise, ending the circulation, and finishing the identification of the platform area topology.
A low pressure HPLC station topology identification system, as shown in fig. 3, comprising: the device comprises a signal acquisition module, a signal processing module, a length calculation module and a topology identification module;
the signal acquisition module inserts the broadband carrier communication module of ammeter and concentrator into this system, includes: the electric meter carrier module and the concentrator carrier module are used for acquiring HPLC signal frequency response one-dimensional waveform data of a receiving end; the data storage module is used for storing the integrated data and uploading the integrated data to the signal processing module;
the signal processing module comprises a first statistical analysis module, a system database and a first data management module; the first statistical analysis module is used for performing wavelet denoising on the signal frequency response acquired by the signal acquisition module, filtering high-frequency noise, performing frequency spectrum analysis and extracting peak point and trapped wave frequency point characteristics; the system database and the first data management module are used for storing data and calling the length calculation module;
the length calculation module comprises a second statistical analysis module and a second data management module; the second statistical analysis submodule is used for calculating branch lengths based on the signal frequency selective attenuation characteristics, calculating line lengths based on the average channel attenuation characteristics, and managing the calculated length data by the second data management module and calling the calculated length data by the topology identification module.
The topology identification module consists of a digital signal processor, a power supply unit, a graphic processor, an MCU unit and a GUI, and is used for calling the length calculation module, determining the number of branch nodes and the upstream-downstream relation, generating a distribution room topological graph and outputting the topological graph.
In order to verify the effectiveness of the low-voltage HPLC station topology identification method based on the broadband power line carrier channel characteristics, the method is applied to a simple typical station established based on a transmission line theory for analyzing the topology identification effect.
As shown in fig. 4, a typical network topology of a cluster of stations is established based on transmission line theory, lengths of each line segment are marked in the figure, where a and B are a transmitting end and a receiving end, nodes 1 to 5 are branch terminal nodes, and loads Z of the transmitting end and the receiving end are respectivelySAnd ZLAre all 50 omega, and all branch terminal nodes are open circuits. The line model used was a cross-linked polyethylene insulated polyvinyl chloride jacketed (YJV) cable, with the line parameters shown in table 1.
TABLE 1 YJV line parameters
Figure BDA0003183392770000101
The topology identification is carried out according to the steps, the frequency response of the signal received by the node B and the amplitude spectrum thereof are shown in fig. 5 and fig. 6, and obvious frequency components are extracted from the amplitude spectrum, wherein the obvious frequency components are 0.8928MHz, 0.9643MHz, 1.679MHz, 1.964MHz, 2.143MHz and 2.357MHz respectively. The generated topology is shown in fig. 7, and the comparison between the calculation result and the actual length of each segment of the line is shown in table 2.
TABLE 2 comparison of theoretical calculation of line length with actual length results
Figure BDA0003183392770000102
As can be seen from fig. 7 and fig. 4, the method provided by the present invention basically realizes accurate identification of the topology of the distribution room, the topology of the distribution room including the branch condition and the upstream and downstream relationships of the nodes is correctly generated, and there are no nodes identified by errors and few nodes identified by errors. As can be seen from the table, the line length error rate is small overall, up to 6.72%. For nodes a and B, the error rate of the calculation of the branch length is small and is controlled within 2%, while the reason that the error of the line length calculated based on the statistical characteristics of the average channel attenuation is large is that the linear relationship ACA (l) of the ACA with the unbranched line length is obtained by actual measurement and approximation using a curve fitting method, wherein the existing error is the main reason for the large error rate of the line length.
It should be emphasized that the examples described herein are illustrative and not restrictive, and thus the present invention includes, but is not limited to, those examples described in this detailed description, as well as other embodiments that can be derived from the teachings of the present invention by those skilled in the art and that are within the scope of the present invention.

Claims (7)

1. A low-pressure HPLC (high performance liquid chromatography) platform area topology identification method is characterized by comprising the following steps: the method comprises the following steps:
step 1, numbering all electric meters according to the distribution area files, drawing in an unidentified electric meter set B1, and initializing a transition electric meter set
Figure FDA0003183392760000011
Identified electricity meter collection
Figure FDA0003183392760000012
Step 2, in a frequency band of 2-30MHz, respectively transmitting and receiving HPLC sweep frequency signals by a carrier electricity meter I in a set B1 and a station area concentrator, sampling to obtain a frequency response sequence H (f), and dividing all electricity meters receiving the carrier signals in the electricity meter I and the set B1 into a set B2; (ii) a
Step 3, performing wavelet denoising on the frequency response sequence;
step 4, determining the number of the transmitting and receiving intermediate branches and the branch length by utilizing the characteristic of the trapped wave frequency point;
step 5, calculating the length of the bus between the transmitting and receiving nodes by using the average channel attenuation characteristic;
step 6, calculating the distance between the carrier wave electric meters in the middle, and adding and deleting middle branch nodes;
and 7, identifying the topology of the transformer area.
2. A low pressure HPLC zone topology identification method according to claim 1, characterized by: the specific method of the step 3 comprises the following steps:
performing wavelet denoising on one-dimensional waveform data obtained by sampling, filtering high-frequency noise, wherein the denoising method selects Daubechies wavelet as wavelet basis, a heursure mixed threshold rule and a soft threshold function, and the threshold calculation rule is as follows:
Figure FDA0003183392760000013
if Eta < Crit, then a fixed threshold is chosen, otherwise the smaller of the rigrsure criterion and the sqtwolog criterion is chosen as the present criterion threshold.
3. A low pressure HPLC zone topology identification method according to claim 1, characterized by: the specific steps of the step 4 comprise:
(1) fast Fourier transform extraction of all significant frequency components F in amplitude spectra1~FnJudging the number of branches to be n;
(2) calculating the notch frequency point interval in the signal frequency response curve corresponding to each branch
Figure FDA0003183392760000021
(3) Calculating the length of each branch:
Figure FDA0003183392760000022
in the formula, L is the distributed inductance of the line, and C is the distributed capacitance of the line.
4. A low pressure HPLC zone topology identification method according to claim 1, characterized by: the specific steps of the step 5 comprise:
(1) firstly, according to the Hankel matrix form, the channel frequency response sequence x is constructed into an m multiplied by n signal matrix HxThe matrix structure is as follows:
Figure FDA0003183392760000023
to HxSingular value decomposition is carried out to obtain sub-channel frequency response H1(f)~Hn(f) The singular value decomposition is in the form:
Hx=USVT (4)
(2) windowed DeltafiAnd performing PE operation to extract the peak points of the frequency division rate response curves, wherein the PE operation is as follows:
Figure FDA0003183392760000024
subtracting the structural element Bm from the original input signal Fn to carry out open operation on the signal Fn, thus obtaining a peak point in each window, fitting a peak curve based on a least square method, and approximately replacing the frequency response of a non-branch channel by the curve;
(3) calculating the average channel attenuation statistical characteristic of the non-branch channel, wherein the calculation formula is as follows:
Figure FDA0003183392760000025
wherein f is0For frequency resolution, take f025kHz, frequency range 2-30MHz, so n1=80,n2=1200;
(4) Calculating the length l of the corresponding main line by ACA (l), wherein the ACA (l) is obtained by outdoor actual measurement, and the measurement rule is as follows:
changing the length l of the non-branch line for many times, and respectively calculating ACA under each l to obtain parameters a and b in the formula (7);
ACA(l)=-al+b(dB) (7)
(5) calculating the main line length l corresponding to each frequency division rate response by the ACA (l) of the formula (3)1~lnAnd the bus length l is obtained by taking the average value.
5. A low pressure HPLC zone topology identification method according to claim 1, characterized by: the specific steps of the step 6 comprise:
(1) collecting the carrier electricity meters i and j in the B2 to receive and transmit HPLC sweep frequency signals, and sampling to obtain Hij(f);
(2) Calculating the line length l between i and j according to step 5ij
(3) When l isij-(lbri+lbrj) When the distance is larger than epsilon, the two electric meters are judged to belong to different branch nodes respectively, and the distance between the branch nodes is lij-(lbri+lbrj);
(4) When lij-(lbri+lbrj) When | ≦ epsilon, judging that the branch nodes to which the two electric meters respectively belong are closer, and merging the branch nodes;
(5) when l isij-(lbri+lbrj) When < -epsilon, two electric meters belong to the same branch, and the branch length is (l)bri+lbrj-lij) Adding middle branch nodes;
(6) and (6) circularly executing the step until the line length calculation among all the electric meters in the B2 is completed.
6. A low pressure HPLC zone topology identification method according to claim 1, characterized by: the specific steps of the step 7 comprise:
(1) calculating the distances between all the electric meters in the B2 and the electric meter I, judging the upstream and downstream relations of the nodes according to the distance, generating a cluster topology, dividing all the elements in the B2 into the B3, and enabling all the elements in the B2 to be in a row
Figure FDA0003183392760000031
The topology identification of the cluster where the ammeter I is located is completed;
(2) when in use
Figure FDA0003183392760000032
And (4) returning to the step (2), otherwise, ending the circulation, and finishing the identification of the platform area topology.
7. A low pressure HPLC district topology identification system is characterized in that: the method comprises the following steps: the device comprises a signal acquisition module, a signal processing module, a length calculation module and a topology identification module;
the signal acquisition module inserts the broadband carrier communication module of ammeter and concentrator into this system, includes: the electric meter carrier module and the concentrator carrier module are used for acquiring HPLC signal frequency response one-dimensional waveform data of a receiving end; the data storage module is used for storing the integrated data and uploading the integrated data to the signal processing module;
the signal processing module comprises a first statistical analysis module, a system database and a first data management module; the first statistical analysis module is used for performing wavelet denoising on the signal frequency response acquired by the signal acquisition module, filtering high-frequency noise, performing frequency spectrum analysis and extracting peak point and trapped wave frequency point characteristics; the system database and the first data management module are used for storing data and calling the length calculation module;
the length calculation module comprises a second statistical analysis module and a second data management module; the second statistical analysis submodule is used for calculating branch length based on the signal frequency selective attenuation characteristic, calculating line length based on the average channel attenuation characteristic, and managing the calculated length data by the second data management module and calling the length data by the topology identification module;
the topology identification module consists of a digital signal processor, a power supply unit, a graphic processor, an MCU unit and a GUI, and is used for calling the length calculation module, determining the number of branch nodes and the upstream-downstream relation, generating a distribution room topological graph and outputting the topological graph.
CN202110853927.XA 2021-07-28 2021-07-28 Low-pressure HPLC (high performance liquid chromatography) platform area topology identification method and system Pending CN113779751A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110853927.XA CN113779751A (en) 2021-07-28 2021-07-28 Low-pressure HPLC (high performance liquid chromatography) platform area topology identification method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110853927.XA CN113779751A (en) 2021-07-28 2021-07-28 Low-pressure HPLC (high performance liquid chromatography) platform area topology identification method and system

Publications (1)

Publication Number Publication Date
CN113779751A true CN113779751A (en) 2021-12-10

Family

ID=78836208

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110853927.XA Pending CN113779751A (en) 2021-07-28 2021-07-28 Low-pressure HPLC (high performance liquid chromatography) platform area topology identification method and system

Country Status (1)

Country Link
CN (1) CN113779751A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114330527A (en) * 2021-12-23 2022-04-12 国网北京市电力公司 Building ammeter distribution identification method, system, device and storage medium
CN117148010A (en) * 2023-09-01 2023-12-01 北京佳悦灏源科技有限公司 HPLC signal characteristic and voltage correlation-based platform region identification method
CN117175572A (en) * 2023-09-08 2023-12-05 北京佳悦灏源科技有限公司 Topology identification method based on HPLC carrier signal combination

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108832626A (en) * 2018-07-06 2018-11-16 北京中电普华信息技术有限公司 A kind of low-voltage power distribution station area topology identification method and system
CN109390934A (en) * 2018-10-12 2019-02-26 江苏佳源科技有限公司 A kind of intelligent low-voltage platform area power circuit automated topology system
CN111462470A (en) * 2020-05-19 2020-07-28 广州穗华能源科技有限公司 Low-voltage distribution area topology identification method based on power line carrier attenuation characteristics
CN111611543A (en) * 2020-04-29 2020-09-01 北京智芯微电子科技有限公司 Low-voltage distribution area user network topology identification method and system
CN112436866A (en) * 2020-11-10 2021-03-02 四川能信科技股份有限公司 Method for calculating circuit topology based on power carrier radiation signal intensity distribution characteristics
WO2021043317A1 (en) * 2019-09-06 2021-03-11 北京市腾河智慧能源科技有限公司 Transformer area identification method and method for constructing transformer area line topography

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108832626A (en) * 2018-07-06 2018-11-16 北京中电普华信息技术有限公司 A kind of low-voltage power distribution station area topology identification method and system
CN109390934A (en) * 2018-10-12 2019-02-26 江苏佳源科技有限公司 A kind of intelligent low-voltage platform area power circuit automated topology system
WO2021043317A1 (en) * 2019-09-06 2021-03-11 北京市腾河智慧能源科技有限公司 Transformer area identification method and method for constructing transformer area line topography
CN111611543A (en) * 2020-04-29 2020-09-01 北京智芯微电子科技有限公司 Low-voltage distribution area user network topology identification method and system
CN111462470A (en) * 2020-05-19 2020-07-28 广州穗华能源科技有限公司 Low-voltage distribution area topology identification method based on power line carrier attenuation characteristics
CN112436866A (en) * 2020-11-10 2021-03-02 四川能信科技股份有限公司 Method for calculating circuit topology based on power carrier radiation signal intensity distribution characteristics

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨志淳;沈煜;杨帆;乐健;宿磊;雷杨;: "基于数据关联分析的低压配电网拓扑识别方法", 电测与仪表, no. 18 *
陆欣;冷安辉;毛珊珊;洪海敏;王春;: "基于高频同步采集与边缘计算的低压配电网电气拓扑发现技术", 电信科学, no. 08 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114330527A (en) * 2021-12-23 2022-04-12 国网北京市电力公司 Building ammeter distribution identification method, system, device and storage medium
CN114330527B (en) * 2021-12-23 2023-10-13 国网北京市电力公司 Building ammeter distribution identification method, system, device and storage medium
CN117148010A (en) * 2023-09-01 2023-12-01 北京佳悦灏源科技有限公司 HPLC signal characteristic and voltage correlation-based platform region identification method
CN117148010B (en) * 2023-09-01 2024-04-05 北京佳悦灏源科技有限公司 HPLC signal characteristic and voltage correlation-based platform region identification method
CN117175572A (en) * 2023-09-08 2023-12-05 北京佳悦灏源科技有限公司 Topology identification method based on HPLC carrier signal combination
CN117175572B (en) * 2023-09-08 2024-04-02 北京佳悦灏源科技有限公司 Topology identification method based on HPLC carrier signal combination

Similar Documents

Publication Publication Date Title
CN113779751A (en) Low-pressure HPLC (high performance liquid chromatography) platform area topology identification method and system
CN111817289B (en) Low-voltage distribution station physical topology identification method, device, system, terminal and medium
CN111624544B (en) Method and related device for identifying household transformer relation of transformer area
CN110865237A (en) Power distribution station area meter phase relation checking method based on voltage time sequence
CN111239549A (en) Power distribution fault rapid positioning method based on discrete wavelet transform
CN112804318B (en) Ammeter region identification method
CN110579690A (en) Positioning method for small current ground fault section of power distribution network
CN112952828A (en) User station area analysis method for low-voltage distribution network
CN111426905A (en) Power distribution network common bus transformation relation abnormity diagnosis method, device and system
CN113092935A (en) Method, system, equipment and storage medium for identifying small-load line topology
CN112948757A (en) Low-voltage distribution area topology verification method based on improved Pearson correlation coefficient
CN111525560A (en) Transformer area household change identification method and device based on current data of intelligent electric meter
CN114912526A (en) Method and system for identifying distribution area, electronic device and storage medium
CN113866552B (en) Medium voltage distribution network user electricity consumption abnormality diagnosis method based on machine learning
CN113466548B (en) Intelligent ammeter area identification method based on phasor measurement technology
CN112415304B (en) Line identification method based on low-voltage transformer area voltage data curve grouping
CN112182499A (en) Low-voltage distribution network topological structure identification method based on time sequence electric quantity data
CN109361484A (en) A kind of transmission method of electric system time synchronization data
CN113992241A (en) Power frequency communication-based automatic identification and analysis method for distribution room topology
CN110927539B (en) Power distribution network fault section positioning method based on single-ended transient energy spectrum similarity
CN117200241A (en) New forms of energy intelligent voltage regulation system that is incorporated into power networks
CN117421622A (en) Low-voltage transformer area line loss calculation method and system based on clustering algorithm
CN110944336A (en) Time-frequency spectrum sensing method based on limited new information rate
CN112436866B (en) Method for calculating circuit topology based on power carrier radiation signal intensity distribution characteristics
CN115659553A (en) Low-voltage power supply network topology identification method and system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination