CN109683524B - Processing method for sampling and synchronizing sampling signals without sampling and holding - Google Patents

Processing method for sampling and synchronizing sampling signals without sampling and holding Download PDF

Info

Publication number
CN109683524B
CN109683524B CN201811595587.XA CN201811595587A CN109683524B CN 109683524 B CN109683524 B CN 109683524B CN 201811595587 A CN201811595587 A CN 201811595587A CN 109683524 B CN109683524 B CN 109683524B
Authority
CN
China
Prior art keywords
sampling
channel
data
calculation
delta
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.)
Active
Application number
CN201811595587.XA
Other languages
Chinese (zh)
Other versions
CN109683524A (en
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.)
Xdge Automation Co ltd
Original Assignee
Xdge Automation Co ltd
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 Xdge Automation Co ltd filed Critical Xdge Automation Co ltd
Priority to CN201811595587.XA priority Critical patent/CN109683524B/en
Publication of CN109683524A publication Critical patent/CN109683524A/en
Application granted granted Critical
Publication of CN109683524B publication Critical patent/CN109683524B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/04Programme control other than numerical control, i.e. in sequence controllers or logic controllers
    • G05B19/042Programme control other than numerical control, i.e. in sequence controllers or logic controllers using digital processors
    • G05B19/0423Input/output
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/20Pc systems
    • G05B2219/25Pc structure of the system
    • G05B2219/25257Microcontroller

Abstract

The invention relates to a processing method for carrying out sampling synchronization on sampling signals without sampling and holding, which comprises the following steps: step 1: measuring analog quantity acquisition channels in the relay protection device to obtain a time interval delta t between each sampling channel; step 2: calculating data of each sampling point to perform interpolation calculation, and deducing a sampling value before delta t; directly compensating the original sampling points by using a secondary interpolation algorithm, storing the original sampling points into a new array, and compensating the sampling data of each channel to the same moment; and step 3: and performing subsequent protection algorithm calculation by using the compensated new data. The invention provides a processing method for synchronizing the data of each channel of a non-sampling hold circuit, which is used for aligning the sampling time of each channel, provides possibility for the calculation of a point level, does not need to change the original circuit design, and saves the cost and the development period.

Description

Processing method for sampling and synchronizing sampling signals without sampling and holding
Technical Field
The invention relates to the technical field of relay protection, in particular to a method for carrying out sampling signal synchronization processing on sampling signals without sampling and holding.
Background
An analog quantity acquisition board in the relay protection device does not perform sampling and holding processing on acquisition signals of all channels at the beginning of design, but sequentially samples data of all channels. However, in the subsequent device upgrading process, sampling data of all paths at the same time is needed, and when a hardware sampling circuit is not allowed due to the limitation of cost and development period, interpolation processing needs to be carried out on the original sampling data, so that a group of sampling data at the same time can be simulated.
Therefore, the problems of development cost and period are solved, and the down compatibility of the produced equipment can be realized.
Disclosure of Invention
The invention aims to provide a software processing method for synchronizing signals of all paths on a non-sampling hold circuit. The invention provides a processing method for synchronizing the data of each channel of a non-sampling hold circuit, which is used for aligning the sampling time of each channel, provides possibility for the calculation of a point level, does not need to change the original circuit design, and saves the cost and the development period.
In order to achieve the purpose, the invention adopts the following technical means:
a processing method for carrying out sampling synchronization on sampling signals without sampling and holding comprises the following steps:
step 1: measuring analog quantity acquisition channels in the relay protection device to obtain a time interval delta t between each sampling channel;
step 2: performing interpolation calculation on each sampling point data, and deducing a sampling value before a time interval delta t; directly compensating the original sampling points by using a secondary interpolation algorithm, storing the original sampling points into a new array, and compensating the sampling data of each channel to the same moment;
and step 3: and performing subsequent protection algorithm calculation by using the compensated new data.
As a further improvement of the present invention, in step 2,
the quadratic interpolation formula is:
Figure BDA0001921270420000021
Figure BDA0001921270420000022
wherein:
x (k-2), x (k-1), x (k) are three consecutive sampling points; y (n) is the interpolated sample point; t (k) and t (k-1) correspond to the sampling time of x (k), x (k-1), tyAnd p is the calculated compensation coefficient.
As a further improvement of the present invention, when the sampling point needs to be moved for a fixed time, the following process is applied to equation (2):
Figure BDA0001921270420000023
wherein, tdraftFor the time the sample point needs to be moved, for the second path t in turndraftIs delta t microsecond, the second channel is 2 delta t microsecond, and so on.
Compared with the prior art, the invention has the following technical effects:
according to the invention, the delta t required to be compensated is calculated, a quadratic interpolation algorithm is used, the original sampling point is directly compensated, and the original sampling point is stored in a new array, so that the sampling data of each channel is compensated to the same moment. For a hardware acquisition circuit without sampling and holding, the acquired data can basically meet the requirement of sampling data synchronization after an interpolation algorithm. The method does not need to change the circuit of the hardware analog quantity acquisition board, can directly use the analog quantity acquisition board which is put into production, saves the cost, shortens the research and development period and has good compatibility.
Drawings
FIG. 1 is a flow chart embodying the present invention;
FIG. 2 is a first diagram of the simulation effect of the algorithm of the present invention;
FIG. 3 is a diagram of the simulation effect of the algorithm of the present invention.
Detailed Description
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, algorithms required to be used in the description of the embodiments will be briefly described below, and it is obvious that the algorithms in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other algorithms based on these algorithms without creative efforts.
In order to make the objects, technical solutions and advantages of the present invention more apparent, the following further description is made of specific embodiments of the present invention, which are intended to be illustrative rather than restrictive.
The core idea of the invention is as follows: and calculating delta t required to be compensated according to the sampling time interval between the hardware channels, directly compensating the original sampling points by using a quadratic interpolation algorithm, and storing the original sampling points into a new array so as to compensate the sampling data of each channel to the same moment. And performing subsequent protection calculation by using the compensated data.
On the premise of not changing the original hardware link, measuring the sampling time difference delta t between each channel, and supplementing the time difference delta t of sampling points between each channel introduced by hardware asynchronous sampling through a software algorithm. The flow is shown in FIG. 1.
Step 1: sampling data;
step 2: calculating data of each sampling point to perform interpolation calculation, and deducing a sampling value before delta t;
and step 3: and performing subsequent protection algorithm calculation by using the compensated new data.
The invention relates to a processing method for synchronizing data of each channel of a non-sampling hold circuit, which is used for aligning sampling time of each channel, provides possibility for calculating point level without changing the original circuit design, saves cost and development period, and comprises the following specific steps:
the method comprises the following steps: and measuring analog quantity acquisition channels in the relay protection device to obtain the time interval delta t between the sampling channels.
Step two: in the program, the original sampling points are compensated according to a quadratic interpolation formula (1).
Figure BDA0001921270420000041
Figure BDA0001921270420000042
Wherein:
x (k-2), x (k-1), x (k) are three consecutive sampling points; y (n) is the interpolated sample point; t (k) and t (k-1) correspond to the sampling time of x (k), x (k-1), tyAnd p is the calculated compensation coefficient.
When the sampling point needs to be moved for a fixed time, the following process can be applied to equation (2):
Figure BDA0001921270420000043
wherein t isdraftFor the time the sample point needs to be moved, for the second path t in turndraftIs delta t microsecond, the second channel is 2 delta t microsecond, and so on。
And step three, the new sampling array obtained after calculation is used for subsequent protection calculation, so that protection algorithms such as zero sequence removal, charging current calculation and the like can be carried out on the point level.
As shown in FIG. 2 and FIG. 3, the simulation effect graphs of the algorithm of the present invention show that the sampling points after interpolation calculation are substantially the same as the original sampling points without phase shift processing, and the compensation effect is good.
The practical test results of the invention are as follows:
the same analog quantity acquisition board (without sampling and holding) is used in the test, three phases are added with the current 2A with the same amplitude, the voltage 69V with the same amplitude, and the theoretical phase difference of the current and voltage phase angle is 90 degrees:
a. the fundamental values calculated without interpolation are shown in table 1:
TABLE 1
Electric current IA IB IC
Amplitude (mA) 2003 2001 2003
Angle (°) 285.677 386.202 268.625
Voltage of VA VB VC
Amplitude (mV) 69005 69136 69080
Angle (°) 17.120 17.567 17.744
From the table, it can be seen that the IA differs from the other channel angles in table 2:
TABLE 2
IA And IB And IC And VA And VB And VC
By an angle (°) different from IA 0.525 0.948 91.443 91.89 92.067
b. The fundamental values calculated using the interpolation algorithm are shown in table 3:
TABLE 3
Figure BDA0001921270420000051
Figure BDA0001921270420000061
From the table, it can be seen that the IA differs from the other channel angles by table 4:
TABLE 4
IA And IB And IC And VA And VB And VC
By an angle (°) different from IA 0.180 0.244 89.868 89.901 90.009
It can be seen from tables 1 to 4 that, under the same test environment, the phase angle difference of each channel using the interpolation algorithm is effectively compensated, so that under the condition of not changing the hardware circuit of the acquisition board, the correction of the acquisition time difference of each channel is realized, and the corrected sampling point can be applied to zero sequence elimination, subsequent amplitude calculation and the like of the sampling point layer.
The foregoing is a detailed description of the present invention with reference to specific preferred embodiments, and no attempt is made to limit the invention to the particular embodiments disclosed, or modifications and equivalents thereof, since those skilled in the art may make various alterations and equivalents without departing from the spirit and scope of the invention, which is defined by the claims appended hereto.

Claims (1)

1. A processing method for sampling and synchronizing sampling signals without sampling and holding is characterized by comprising the following steps:
step 1: measuring analog quantity acquisition channels in the relay protection device to obtain a time interval delta t between each sampling channel;
step 2: performing interpolation calculation on each sampling point data, and deducing a sampling value before a time interval delta t; directly compensating the original sampling points by using a secondary interpolation algorithm, storing the original sampling points into a new array, and compensating the sampling data of each channel to the same moment;
and step 3: performing zero sequence removal and charging current calculation on a point level by using the compensated new data;
in step 2, the quadratic interpolation formula is:
Figure FDA0002915400140000011
Figure FDA0002915400140000012
wherein:
x (k-2), x (k-1), x (k) are three consecutive sampling points; y (n) is the interpolated sample point; t (k) and t (k-1) correspond to the sampling time of x (k), x (k-1), tyThe time of the interpolation target point is, and p is a compensation coefficient obtained by calculation;
when the sampling point needs to be moved for a fixed time, the following process is applied to equation (2):
Figure FDA0002915400140000013
wherein, tdraftFor the time the sample point needs to be moved, for the second path t in turndraftIs delta t microsecond, the second channel is 2 delta t microsecond, and so on.
CN201811595587.XA 2018-12-25 2018-12-25 Processing method for sampling and synchronizing sampling signals without sampling and holding Active CN109683524B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811595587.XA CN109683524B (en) 2018-12-25 2018-12-25 Processing method for sampling and synchronizing sampling signals without sampling and holding

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811595587.XA CN109683524B (en) 2018-12-25 2018-12-25 Processing method for sampling and synchronizing sampling signals without sampling and holding

Publications (2)

Publication Number Publication Date
CN109683524A CN109683524A (en) 2019-04-26
CN109683524B true CN109683524B (en) 2021-05-11

Family

ID=66189441

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811595587.XA Active CN109683524B (en) 2018-12-25 2018-12-25 Processing method for sampling and synchronizing sampling signals without sampling and holding

Country Status (1)

Country Link
CN (1) CN109683524B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110750752B (en) * 2019-09-10 2023-12-05 许昌许继软件技术有限公司 Interpolation method and device for analog data
CN114609997B (en) * 2022-03-07 2023-06-27 南京国电南自维美德自动化有限公司 Method for realizing analog input change rate protection of process control station
CN115269593B (en) * 2022-07-20 2023-04-18 清云智通(北京)科技有限公司 Industrial data acquisition and storage method, device and system

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101762742A (en) * 2009-12-16 2010-06-30 南京弘毅电气自动化有限公司 Sampling value correcting method based on fixed-frequency sampling sequence
CN103605023B (en) * 2013-11-20 2016-09-21 国家电网公司 A kind of combining unit time response measuring method and measurement apparatus
CN103969500B (en) * 2014-04-11 2017-09-12 深圳市中电电力技术股份有限公司 Intelligent grid analog acquisition device, method and system
CN103986129B (en) * 2014-05-30 2018-04-20 广西电网公司 Bus bar protecting method based on voltage lock phase networking sampling
US10201059B1 (en) * 2017-08-01 2019-02-05 Kleverness Incorporated Method for analyzing operating parameters for lighting technologies

Also Published As

Publication number Publication date
CN109683524A (en) 2019-04-26

Similar Documents

Publication Publication Date Title
CN109683524B (en) Processing method for sampling and synchronizing sampling signals without sampling and holding
CN103067006B (en) A kind of real-time correction method for time-interleaved A/D conversion system time error
CN102253358A (en) Method for correcting electric energy meter by using correction system, and correction system
CN104375006B (en) A kind of Fast synchronization phase amount correction method
CN105259410A (en) Under-sampling waveform frequency estimation method and device under strong noise interference
CN110837003B (en) Double-window full-phase DFT (discrete Fourier transform) synchronous phasor measurement method and system based on triangular window
CN109581018A (en) Wide-band sampling oscillograph time base error synchroballistic method and device
CN111625769A (en) PMU-SCADA data time synchronization and fusion method based on Lagrange interpolation and cubic exponential smoothing
CN111896802B (en) Frequency self-adaptive sampling method
CN115508619A (en) Inverter impedance measuring method and device, terminal device and storage medium
CN112087232A (en) Method and device for calibrating direct current offset and gain error of TIADC system
CN103605904B (en) Self compensation power system amplitude arithmetic based on error estimation
CN203057112U (en) Real-time correction device according to time errors of time-interleaved digital-to-analog conversion system
US8866517B2 (en) Power system data acquisition systems
CN109490625B (en) Harmonic signal analysis method based on sliding window and semi-definite programming
CN107064610A (en) Being averaged based on group delay
CN105846822A (en) Method and system for measuring mismatch error of time-interleaved analog-to-digital conversion system
CN108982966A (en) Humorous phase angle analysis method based on linear correction algorithm
Zeng et al. Real-time self-adaptive calibration method for high speed acquisition system
Wenxiang et al. Phase measurement about adaptive notch filter based on LMS algorithm
CN115561517B (en) Method, device and equipment for determining split-phase electric quantity of three-phase electric quantity
CN107968396B (en) Thevenin equivalent parameter calculation method based on local phasor measurement information
CN113065304B (en) Loop delay estimation method based on sliding window
CN102045078B (en) FPGA (Field Programmable Gate Array) based software receiver system and implementation method
CN110244116B (en) DC instantaneous power metering circuit and quasi-synchronous calculation method thereof

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
GR01 Patent grant
GR01 Patent grant