Disclosure of Invention
In order to solve the defects and shortcomings in the prior art, the method for testing the attenuation time constant of the integration loop of the acquisition unit is provided, the test is carried out in an alternating current-direct current superposition synchronous output mode, the test result is obtained, and meanwhile the influence of discrete test errors on the accuracy of the test result can be avoided.
Specifically, the method for testing the decay time constant of the integration loop of the acquisition unit provided by this embodiment includes:
generating an alternating current and direct current superposition test waveform at an upper computer, and processing data of the alternating current and direct current superposition test waveform to obtain a small-voltage analog signal;
acquiring an optical fiber digital signal output by an integration loop of the acquisition unit, and testing by combining a small-voltage analog signal to obtain an alternating current-direct current transmission ratio;
and calculating the attenuation time constant based on the AC/DC transmission/transformation ratio.
Optionally, an alternating current-direct current superposition test waveform is generated at the upper computer, and data of the alternating current-direct current superposition test waveform are processed to obtain a small-voltage analog signal, including:
setting preset parameters of alternating current amplitude, angular frequency and direct current voltage value at the upper computer, and generating an alternating current and direct current superposition test waveform based on the preset parameters;
and carrying out discretization processing on the data of the alternating current-direct current superposition test waveform, and carrying out D/A conversion to obtain a small-voltage analog signal.
Optionally, the host computer department sets up the preset parameter of alternating current amplitude, angular frequency and direct current voltage value, generates the alternating current-direct current stack test waveform based on preset parameter, includes:
based on a formula, generating a function expression of the AC/DC superposition test waveform,
v(t)=Umsin(ωt)+U0the formula I is shown in the specification,
where v (t) is the output amplitude of the superimposed waveform, UmIs the AC amplitude, omega is the angular frequency, the calculation expression is omega-2 pi f, f is the AC frequency, U0Is a DC value, and U is used for testing to prevent testing error caused by DC overflow after integration0Less than 10 mV.
Optionally, the discretization processing is performed on the data of the alternating current-direct current superposition test waveform, and a small-voltage analog signal is obtained after the D/a conversion, including:
the data of the alternating current-direct current superposition test waveform is subjected to discretization processing as shown in a formula II,
v(n)=v(nTs)=Umsin(ωnTs)+U0in the formula II, the first step is carried out,
where v (n) is the discretization output result, TsThe value range of n is a positive integer for discretized sampling isolation time.
Optionally, acquiring the optical fiber digital signal output by the acquisition unit integrating loop, and testing by combining the small-voltage analog signal to obtain the ac/dc transmission/transformation ratio includes:
collecting analog quantity signals and optical fiber digital quantity signals from a test sample by adopting an FPGA chip;
respectively extracting alternating current and direct current signals of the standard signals and the signals to be detected comprising the collected optical fiber digital quantity signals and the collected small voltage analog signals;
respectively calculating the transmission ratio k of the DC signals according to the extraction results1Transmission ratio k to AC signal2。
Optionally, the method further includes:
after the FPGA chip is adopted to collect the analog quantity signal and the optical fiber digital quantity signal from the test article, the alternating current frequency f is calculated by a fundamental wave frequency measurement method based on an FFT algorithm.
Optionally, the ac/dc signal extraction is performed on the standard signal and the signal to be measured including the collected optical fiber digital quantity signal and the collected small-voltage analog signal, respectively, and includes:
after frequency measurement, Hanning window processing is respectively carried out on the standard signal and the signal to be measured, Fourier transform is carried out on a windowing sequence, and the amplitude and the phase of an alternating current signal are obtained;
and extracting alternating current 8-cycle time window data of the direct current component, and adopting generalized polynomial fitting to calculate the amplitude of the direct current component to eliminate harmonic signals and noise influence.
Optionally, the transmission ratio k of the dc signal is calculated according to the extraction result1Transmission ratio k to AC signal2The method comprises the following steps:
obtaining a transmission calculation formula on a direct current frequency domain and an alternating current frequency domain:
coefficient of direct current transmission1,k1=R2/R1,
Coefficient of transmission of alternating current k2,k2≈1/ωC/R1≈1/ωCR1,
Where C is the capacitance of the integrating output capacitor and R1For the resistance value, R, in series with the integrating output capacitor C in a lossy integrator circuit2The lossy integrator circuit has a resistance value connected in parallel to the integrating output capacitor C.
Optionally, the calculating a decay time constant based on the ac/dc transmission/transformation ratio includes:
transmission ratio k based on direct current signal1Transmission ratio k to AC signal2Calculating the decay time constant by combining the formula III,
Tb=RC≈k1/ωk2≈k1/2πfk2in the formula three, the first step is,
in the formula, TbIs the decay time constant of the integration loop, C is the capacitance of the loop integration capacitor, k1Is a direct current transmission coefficient, k2Is the AC transmission coefficient, and f is the AC frequency.
The beneficial effect that technical scheme that this application provided brought is:
the test is carried out by adopting an alternating current and direct current superposition synchronous output mode, so that the influence of discrete test errors on the accuracy of a test result is avoided. Meanwhile, the test waveform is generated and then discretized, so that discrete errors in the waveform superposition process are avoided.
Example one
At present, electronic transformers mainly use electronic transformers, and the electronic transformers are used for electronic transformers of a relay protection circuit based on the principle of Rogowski coil. The basic principle of Rogowski coil for measuring current is Faraday's law of electromagnetic induction and ampere's loop law, the Rogowski coil is an air-core coil formed by using non-magnetic material as skeleton, and in the air-core coil, it is wound on the non-magnetic skeleton for the second time, and the non-ferromagnetic material can make the linearity of said sensor good, and its unsaturation and hysteresis-free are used, so that the air-core coil possesses excellent frequency response capability.
However, the voltage v (t) output by the rogowski coil for the second time and the primary current ip (t) are in a differential relation, and the electronic transformer based on the rogowski coil principle needs to be restored through an integration link because the output is a signal after differentiation. In order to eliminate the influence of an ideal integrator on infinite integration of a direct current signal, a lossy integrator is generally adopted, that is, a resistor R2 is connected in parallel to an integrating output capacitor C to form a direct current attenuation loop to realize suppression of an output interference signal of the rogowski coil.
The circuit structure of the lossy integrator is shown in fig. 1, and the output voltage signal of the rogowski coil is:
where the proportionality coefficient M is the mutual inductance between the rogowski coil and the current carrying conductor, a constant related only to the structure of the rogowski coil, and i (t) is the measured current.
When the system has a short-circuit fault, the final transient fault current of the electronic transformer is reduced through the differential and integral loops of the rogowski coil as follows:
in the formula I
mIs the magnitude of the periodic component of the fault current,
alpha is fault voltage phase; phi is the angle of the impedance,
T
atime constant, T, of short-circuit loop
aL/R, L and R are short-circuit loop inductance and resistance, respectively, T
bThe integral of the decay time constant of the integral loop is the product T of the integral capacitance and the decay resistance of the loop
b=RC
2. In high voltage networks, the inductance value of the fault loop is typically much larger than the resistance value, which can be approximated as φ ≈ 90 °, so its transient error current is approximately equal to:
from the above formula, the transient state transmission error current of the electronic transformer is mainly determined by the attenuation constant of the system and the attenuation constant of the integral loop, and only when T isb>>TaThe transient error caused by the integration loop can be ignored. So the decay time constant T of the integration loopbIs a main index for determining the transient characteristic of the electronic current transformer based on the Rogowski coil principle. But the integration loop is designed inside the acquisition unit and cannot be obtained through a separate test.
In order to solve the influence of the transient characteristics of the electronic transformer based on the rogowski coil principle, the embodiment of the application provides a method for testing the attenuation time constant of an integration loop of an acquisition unit, as shown in fig. 1, the method comprises the following steps:
11. generating an alternating current and direct current superposition test waveform at an upper computer, and processing data of the alternating current and direct current superposition test waveform to obtain a small-voltage analog signal;
12. acquiring an optical fiber digital signal output by an integration loop of the acquisition unit, and testing by combining a small-voltage analog signal to obtain an alternating current-direct current transmission ratio;
13. and calculating the attenuation time constant based on the AC/DC transmission/transformation ratio.
In implementation, an alternating current-direct current superposition test waveform is generated through an upper computer, test data are discretized and then converted into a small-voltage analog signal through a D/A conversion module, then the test system is constructed by acquiring the analog signal and simultaneously acquiring an optical fiber digital signal output by an acquisition unit, an alternating current-direct current input-output ratio is obtained through testing, a test result is uploaded to the upper computer, and an integral loop of the test result is calculated by the upper computer to obtain an attenuation time constant.
Specifically, the specific content executed in step 11 includes:
111. setting preset parameters of alternating current amplitude, angular frequency and direct current voltage value at the upper computer, and generating an alternating current and direct current superposition test waveform based on the preset parameters;
112. and carrying out discretization processing on the data of the alternating current-direct current superposition test waveform, and carrying out D/A conversion to obtain a small-voltage analog signal.
In implementation, generating an ac-dc superimposed test waveform includes:
based on a formula, generating a function expression of the AC/DC superposition test waveform,
v(t)=Umsin(ωt)+U0the formula I is shown in the specification,
where v (t) is the output amplitude of the superimposed waveform, UmIs the AC amplitude, omega is the angular frequency, the calculation expression is omega-2 pi f, f is the AC frequency, U0Is a DC value, and U is used for testing to prevent testing error caused by DC overflow after integration0Less than 10 mV.
The data with alternating current-direct current stack test waveform carries out the discretization and handles to obtain little voltage analog signal after carrying out D/A conversion, include:
the data of the alternating current-direct current superposition test waveform is subjected to discretization processing as shown in a formula II,
v(n)=v(nTs)=Umsin(ωnTs)+U0in the formula II, the first step is carried out,
where v (n) is the discretization output result, TsFor discretized sample isolation time, the value of nThe range is a positive integer.
In the implementation, the D/A conversion chip adopts AD5683R of ADI company, which is a 16-bit single-channel converter, the relative accuracy of the converter is +/-2 LSB INL, and a 2 ppm/DEG C2.5V reference voltage source is arranged in the converter; the space-saving 2mm multiplied by 2mm 8-pin LFCSP and 10-pin MSOP packaging is adopted, so that more functions can be realized in smaller circuit board space; 2mV total non-adjustment error, without initial calibration or adjustment; 4kV HBM ESD rating, system robustness is achieved.
The step of obtaining the transmission ratio of the alternating current to the direct current provided in the step 12 includes:
121. collecting analog quantity signals and optical fiber digital quantity signals from a test sample by adopting an FPGA chip;
122. respectively extracting alternating current and direct current signals of the standard signals and the signals to be detected comprising the collected optical fiber digital quantity signals and the collected small voltage analog signals;
123. respectively calculating the transmission ratio k of the DC signals according to the extraction results1Transmission ratio k to AC signal2。
In the implementation, step 121 specifically includes: after the FPGA chip is adopted to collect the analog quantity signal and the optical fiber digital quantity signal from the test article, the alternating current frequency f is calculated by a fundamental wave frequency measurement method based on an FFT algorithm.
In addition to alternating current-direct current signal extraction, still include:
1221. respectively carrying out Hanning window processing on the standard signal and the signal to be detected, and carrying out Fourier transform on a windowing sequence to obtain the amplitude and the phase of the alternating current signal;
1222 extracting AC 8 cycle time window data of the DC component, adopting generalized polynomial fitting to calculate the amplitude of the DC component, and eliminating the harmonic signal and noise effect.
Combining the steps, the step of finally obtaining the transmission ratio comprises the following steps:
1223. obtaining a transmission calculation formula on a frequency domain:
coefficient of direct current transmission1,k1=R2/R1,
AC transmission systemNumber k2,k2≈1/ωC/R1≈1/ωCR1,
Where C is the capacitance of the integrating output capacitor and R1For the resistance value, R, in series with the integrating output capacitor C in a lossy integrator circuit2The lossy integrator circuit has a resistance value connected in parallel to the integrating output capacitor C.
Finally, the step 13 of calculating the decay time constant based on the ac/dc transfer ratio includes:
transmission ratio k based on direct current signal1Transmission ratio k to AC signal2Calculating the decay time constant by combining the formula III,
Tb=RC≈k1/ωk2≈k1/2πfk2in the formula three, the first step is,
in the formula, TbIs the decay time constant of the integration loop, C is the capacitance of the loop integration capacitor, k1Is a direct current transmission coefficient, k2Is the AC transmission coefficient, and f is the AC frequency.
In practice, the decay time constant T is obtained in order to obtain an integration loopb=RC2The method comprises the following steps of constructing a function derivation formula, wherein an integral loop is formed by an integral resistor, an integral capacitor, a direct current attenuation resistor and an operational amplifier, and the ideal integral formula is as follows:
in the formula: i (0) is the initial value of i (t)
However, in the amplitude-frequency domain, the characteristic of the integral loop is that the transmission of the fundamental wave cannot be changed due to the DC attenuation, so R2>>1/ω C. Therefore, the following two transformation equations are obtained in the dc and ac frequency domains:
coefficient of direct current transmission1=R2/R1K is the coefficient of transmission of alternating current2≈1/ωC/R1≈1/ωCR1,
Two formulas are divided: k is a radical of1/k2≈(R2/R1)/(1/ωCR1)≈ωCR2The decay time constant can be obtained: t isb=RC2≈k1/ωk2≈k1/2πfk2。
The invention relates to a test scheme for an integral loop decay time constant of an electronic transformer based on the Rogowski coil principle, which aims to improve the transient capability of the electronic transformer and provide guarantee for the safe and stable operation of a power grid. Therefore, the invention has the following technical characteristics:
1. and the alternating current and direct current superposition test is carried out by adopting an alternating current and direct current superposition synchronous output mode, so that the influence of discrete test errors on the accuracy of a test result is avoided.
2. Discretization processing is carried out, and discretization processing is carried out after test waveforms are generated, so that discrete errors in the waveform overlapping process are avoided.
3. Compatibility, because the output protocol of the collector of the electronic transformer manufacturer is basically a proprietary protocol, the self-adaptive analysis of the protocols of multiple manufacturers is supported.
4. And (4) standard source recovery, wherein the standard source is recovered in order to ensure the test precision of the whole test system, so that the precision problem possibly existing in the data source is avoided.
5. Adaptability, the coefficients of the digital quantity and the analog quantity in the acquisition units of various manufacturers are inconsistent, and the input-output ratio is adopted for calculation without concerning the setting of the coefficients in the acquisition units of the manufacturers.
6. The method is easy to realize, the attenuation time constant is a transient index, the testing precision is greatly improved by calculation after the transient index is tested by adopting a steady state method, the accurate measurement of the attenuation time constant can be realized, and the testing method is simple and easy to realize.
The sequence numbers in the above embodiments are merely for description, and do not represent the sequence of the assembly or the use of the components.
The above description is only exemplary of the present application and should not be taken as limiting the present application, as any modification, equivalent replacement, or improvement made within the spirit and principle of the present application should be included in the protection scope of the present application.