Distributed laboratory system punctuality method, apparatus, electronic device, and storage medium

Abstract

Embodiments of the present application disclose a time-keeping method and device for a distributed laboratory system, electronic equipment and a storage medium. Calibration data of a coordinated universal time traceability module and clock difference data and corresponding stabilities of at least two time-keeping clocks in the distributed laboratory system are obtained. A local atomic time is determined according to the clock difference data and the corresponding stabilities. It is determined whether the coordinated universal time traceability module is operating normally. If yes, the local atomic time is calibrated according to the calibration data to obtain a standard time. If not, reference data of a reference clock and prediction data of a virtual coordinated universal time prediction module and corresponding uncertainties thereof are obtained, and the local atomic time is calibrated according to the reference data and the prediction data and the corresponding uncertainties to obtain the standard time. Embodiments of the present application improve the disaster recovery performance and safety of the generation of the standard time of the distributed laboratory system.

Description

TECHNICAL FIELD

[0001] Embodiments of the present application relate to time and frequency metrology technology, and in particular to a distributed laboratory system time keeping method and device, electronic equipment and storage medium. BACKGROUND

[0002] Time and frequency are widely used in scientific research, national economic construction and national security, and at the same time affect various aspects of people's daily life.

[0003] In the prior art, the standard time resources of each laboratory are relatively concentrated, have poor disaster recovery performance and low security. SUMMARY

[0004] The present application provides a distributed laboratory system time keeping method and device, electronic equipment and storage medium to improve the disaster recovery performance and security of the standard time generated by the distributed laboratory system.

[0005] In a first aspect, embodiments of the present application provide a distributed laboratory system time keeping method applied to each laboratory node in a distributed laboratory system, which comprises:

[0006] obtaining calibration data of a coordinated universal time traceability module and clock difference data of at least two time keeping clocks in the distributed laboratory system;

[0007] determining a local atomic time according to each clock difference data and its corresponding stability;

[0008] judging whether the coordinated universal time traceability module is running normally;

[0009] if so, calibrating the local atomic time according to the calibration data to obtain a standard time;

[0010] otherwise, obtaining reference data of a reference clock and prediction data of a virtual coordinated universal time prediction module and their corresponding uncertainties, and calibrating the local atomic time according to the reference data and the prediction data and their corresponding uncertainties to obtain the standard time.

[0011] In a second aspect, embodiments of the present application also provide a distributed laboratory system time keeping device configured in each laboratory node in a distributed laboratory system, which comprises:

[0012] a data acquisition module for obtaining calibration data of a coordinated universal time traceability module and clock difference data of at least two time keeping clocks in the distributed laboratory system;

[0013] a local atomic time generation module for determining a local atomic time according to each clock difference data and its corresponding stability;

[0014] The traceability judgment module is configured to judge whether the coordinated universal time traceability module is in normal operation.

[0015] The standard time correction module is configured to, if yes, calibrate the local atomic time according to the calibration data to obtain the standard time.

[0016] The standard time prediction correction module is configured to, if no, acquire reference data of the reference clock and prediction data of the virtual coordinated universal time prediction module and corresponding uncertainties of the reference data and the prediction data, and calibrate the local atomic time according to the reference data and the prediction data and the corresponding uncertainties to obtain the standard time.

[0017] In a third aspect, an electronic device is also provided, and the electronic device includes:

[0018] one or more processors;

[0019] a storage device configured to store one or more programs;

[0020] When the one or more programs are executed by the one or more processors, the one or more processors implement any one of the time-keeping methods for a distributed laboratory system provided by the embodiments of the present application.

[0021] In a fourth aspect, a storage medium including computer executable instructions is also provided, and the computer executable instructions, when executed by a computer processor, are configured to perform any one of the time-keeping methods for a distributed laboratory system provided by the embodiments of the present application.

[0022] The present application acquires calibration data of a coordinated universal time traceability module and clock difference data of at least two time-keeping clocks in a distributed laboratory system, determines local atomic time according to each clock difference data and corresponding stability, judges whether the coordinated universal time traceability module is in normal operation, calibrates the local atomic time according to the calibration data to obtain standard time if yes, acquires reference data of a reference clock and prediction data of a virtual coordinated universal time prediction module and corresponding uncertainties of the reference data and the prediction data, and calibrates the local atomic time according to the reference data and the prediction data and the corresponding uncertainties to obtain the standard time if no, timely corrects frequency and phase through the reference data and the prediction data of the coordinated universal time to determine the standard time when the coordinated universal time is abnormal, and each laboratory node in the distributed laboratory system has standard time generation correction capability and time data resources, so that the decentralization of standard time generation is realized, and the disaster recovery and security of standard time determination are improved. Therefore, the technical solution of the present application solves the problems of relatively concentrated standard time resources of each laboratory, poor disaster recovery performance, and low security, and achieves the effect of improving the disaster recovery performance and security of standard time generation of the distributed laboratory system. BRIEF DESCRIPTION OF DRAWINGS

[0023] Figure 1 is a flow chart of a distributed laboratory system time keeping method in Embodiment One of the present application;

[0024] Figure 2 is a flow chart of a distributed laboratory system time keeping method in Embodiment Two of the present application;

[0025] Figure 3 is a flow chart of a distributed laboratory system time keeping method in Embodiment Three of the present application;

[0026] Figure 4 is a structural schematic diagram of a distributed laboratory system time keeping device in Embodiment Four of the present application;

[0027] Figure 5 is a structural schematic diagram of an electronic device in Embodiment Five of the present application. DETAILED DESCRIPTION

[0028] In order to make the personnel in the technical field better understand the present application, the technical solutions in the embodiments of the present application will be described clearly and completely below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor should belong to the scope of protection of the present application.

[0029] It should be noted that the terms "first" and "second" and the like in the specification and claims of the present application and the above-mentioned drawings are used to distinguish similar objects, and do not necessarily have to describe a specific order or sequence. It should be understood that the data thus used can be interchanged under appropriate circumstances, so that the embodiments of the present application described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "include" and "have" and any variations thereof are intended to cover non-exclusive inclusion, for example, a process, method, system, product or device that includes a series of steps or units does not have to be limited to those steps or units clearly listed, but can include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0030] Embodiment One

[0031] Figure 1 A flow chart of a distributed laboratory system time keeping method provided for Embodiment One of the present application, the present embodiment can be applicable to the case where each laboratory node in a distributed laboratory system determines a standard time, which method can be executed by a distributed laboratory system time keeping device, which device can be implemented by software and / or hardware and specifically configured in each laboratory node in the distributed laboratory system.

[0032] Referring to Figure 1 The distributed laboratory system time keeping method shown is applied to each laboratory node in a distributed laboratory system, and specifically includes the following steps:

[0033] S110, calibration data of a coordinated universal time tracing module and clock difference data of at least two time keeping clocks in the distributed laboratory system are obtained.

[0034] The distributed laboratory system can be composed of at least two laboratory nodes. Each laboratory node is equipped with an atomic clock group meeting the requirements of the International Bureau, including a reference clock and a time keeping clock. For example, the reference clock can be at least one of a cesium fountain clock and an optical clock. For example, the time keeping clock can be at least one of a rubidium clock, a cesium clock, and a hydrogen clock.

[0035] The coordinated universal time tracing module is used to obtain the deviation of the local standard time from the coordinated universal time (UTC) through the Circular T bulletin (a professional term, a kind of bulletin name) published by the International Bureau downloaded automatically, to ensure that the coordinated universal time tracing uncertainty is controlled within 10 nanoseconds, meeting the requirements of the International Bureau on the laboratory. The calibration data can be the deviation of the local standard time from the coordinated universal time, which can include frequency difference and phase difference, and is used to correct the local standard time.

[0036] The clock difference data can be the frequency difference and phase difference between different laboratory nodes in the distributed laboratory system, and is used to generate the local atomic time. The stability can be an index for measuring the reliability of the clock difference data, and is used to determine the weight of the clock difference data when generating the local atomic time. The laboratory nodes realize high-precision time and frequency comparison through satellite two-way time and frequency transfer (TWSTFT), optical fiber time and frequency transfer, and GNSS time and frequency transfer based on the International GNSS (a professional term, a satellite navigation system name) service organization (IGS) precise ephemeris, to obtain the clock difference data between any two reference clocks and / or time keeping clocks.

[0037] All nodes are interconnected through a redundant secure data network, clock difference data transmission uses the network time security (NTS) protocol defined in IETF RFC 8575 (a professional term, a data transmission protocol), and clock difference data is encapsulated in the BIPM TT-xxx (a professional term, a data model) data model recommended by the International Time Frequency Consultative Committee (CCTF). The system adopts a distributed storage architecture, each node is provided with a primary and backup storage node and real-time synchronization of data copies, and the database has forward error correction and automatic data recovery capability.

[0038] In an optional embodiment, the data network adopts a weight-based multi-copy backup strategy to store the data of the reference clock and / or time keeping clock of each laboratory.

[0039] The data backup quantity of the reference clock and / or the time-keeping clock is determined according to the weight values of the performance indicators of the reference clock and / or the time-keeping clock. The data backup quantity Ni is proportional to the weight value Wi, and satisfies Ni = [k × Wi] + 1, where k is a preset coefficient, and [ ] represents an integral operation. The backup data storage location is dynamically optimized according to the network delay. The weight value is dynamically calculated according to the performance indicators of the reference clock and / or the time-keeping clock of each laboratory, and is updated regularly. The minimum value of the backup quantity Ni is 2, to ensure data redundancy, and the maximum value does not exceed the upper limit Nmax preset by the system; Nmax can be determined by professional technicians according to experience or test, and is not limited in the present application.

[0040] The data of the high-weight laboratory is preferentially stored in a node with low delay and high reliability, to improve the data access efficiency. The backup strategy also includes a data verification mechanism, to ensure the consistency of the backup data and the original data. When the weight value of a laboratory changes, the system automatically adjusts the data backup quantity and the storage space of the laboratory. The backup strategy supports incremental backup, and only synchronizes the data of the reference clock and / or the time-keeping clock that changes, to reduce the network load.

[0041] S120, determining a local atomic time according to the clock difference data and the corresponding stabilities of the clock difference data.

[0042] The local atomic time can be obtained by weighted averaging each atomic clock in the atomic clock group according to the clock difference data and the corresponding stabilities of the clock difference data.

[0043] In an optional embodiment, the fusion weight of each clock difference data is determined according to the corresponding stability of each clock difference data; and the local atomic time is determined according to each clock difference data and the corresponding fusion weight of each clock difference data.

[0044] The fusion weight can be the weight when weighted averaging each clock difference data, and is used to determine the local atomic time. Optionally, after obtaining each clock difference data, the clock difference data collected across nodes is subjected to outlier screening and environmental disturbance compensation preprocessing. The fusion weight of each clock difference data is determined according to the corresponding stability of each clock difference data. The higher the stability, the higher the fusion weight, to ensure the stability of the time to be calibrated. Each clock difference data and the corresponding fusion weight are weighted and averaged to obtain the local atomic time.

[0045] S130, judging whether the time coordination module of the coordinated universal time is normally running.

[0046] When the time coordination module of the coordinated universal time is normally running, the standard time is traced to the coordinated universal time; when the time coordination module of the coordinated universal time is abnormally running, the standard time is traced to the weighted average of the predicted data of the coordinated universal time and the reference data of the reference clock, so it is necessary to judge whether the time coordination module of the coordinated universal time is normally running.

[0047] The abnormal situation of the UTC traceability module can include the following situations: GNSS signal loss for more than 120 minutes, UTC deviation for more than ± 20 ns for 30 minutes, and BIPM connection interruption for 12 hours, etc. When the UTC traceability module meets any abnormal situation, it can be determined that the UTC traceability module is running abnormally.

[0048] S140, if yes, calibrating the local atomic time according to the calibration data to obtain the standard time.

[0049] If yes, that is, the UTC traceability module is running normally, at this time, the local atomic time is calibrated according to the calibration data, and the standard time can be determined. Exemplarily, the calibration data can be frequency and phase, which can be predicted first, and then added to the frequency and phase of the local atomic time to calibrate the local atomic time and obtain the standard time.

[0050] The accuracy of the local atomic time is further improved through further calibration of the calibration data. After obtaining the standard time, the digital fine-tuning phase technology can be used to convert the calculation result into a physical signal output, providing a double-redundant 1PPS signal and a 10MHz frequency output.

[0051] S150, otherwise, obtaining reference data of a reference clock and prediction data of a virtual UTC prediction module and their corresponding uncertainties, and calibrating the local atomic time according to the reference data and the prediction data and their corresponding uncertainties to obtain the standard time.

[0052] Otherwise, that is, the UTC traceability module is running abnormally, at this time, the local atomic time needs to be calibrated according to the reference data and the prediction data and their corresponding uncertainties to determine the standard time. The reference clock is used to correct the time frequency and the time phase. The operation of the reference clock can be intermittent, but the frequency of the reference clock is very stable, so it can be used to correct the time frequency. The reference data can be the frequency and phase of the reference clock, which is used to correct the local time to obtain the standard time.

[0053] The virtual UTC prediction module is used to predict the UTC when the UTC traceability module is running abnormally, to obtain prediction data, which can be used to correct the frequency and phase of the local time with the UTC when the UTC traceability module is running abnormally, to obtain the standard time. Exemplarily, the virtual UTC prediction module predicts the UTC according to the phase difference and frequency difference of the historical UTC to obtain the prediction data.

[0054] According to the benchmark reference data and the prediction data and their corresponding uncertainties, the local atomic time is calibrated to obtain the standard time. The standard time can be determined according to the benchmark reference data and the prediction data, and the local atomic time is calibrated by the fused benchmark data to obtain the standard time. Wherein, the fused reference data can be determined according to the benchmark reference data and the prediction data, and the corresponding weight can be determined according to the uncertainty of the benchmark reference data and the prediction data, and the benchmark reference data and the prediction data are weighted and averaged to obtain the fused benchmark data. At this time, the role of the prediction data is to ensure that the calibration reference does not jump when it transitions from coordinated universal time to the benchmark clock, and to ensure the smooth generation of the standard time. Further, the local atomic time is calibrated by the fused benchmark data to obtain the standard time.

[0055] It should be noted that the related information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data for display, analyzed data, etc.) involved in the present disclosure are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of related data comply with relevant laws, regulations and standards in the relevant region.

[0056] Time and frequency are widely used in scientific research, national economic construction and national security, and at the same time affect all aspects of people's daily life. In scientific research, precise measurement of time and frequency is indispensable support for physics, astronomy, automatic control, geophysics, and even life science research. In the national economic construction, many basic industries, such as transportation, energy and power, finance, communication, network, computer, etc., are highly dependent on accurate measurement of time; in the national security construction, satellite navigation positioning, aerospace, information warfare, precision strike, precise time and frequency measurement plays an irreplaceable role in determining the outcome of the war. The most direct and typical application of time and frequency is satellite positioning system, such as the global positioning system (GPS) of the United States, the GLONASS system of Russia, the GALILEO system of Europe and the BD system of China. It is not an exaggeration to say that the satellite positioning system has profoundly changed the way of human production and life, and greatly promoted the ability of human beings to understand and transform the world.

[0057] Since 1972, Coordinated Universal Time (UTC) has been adopted internationally as a global unified time. UTC is based on about 400 commercial hydrogen atomic clocks and cesium atomic clocks distributed in dozens of laboratories around the world (including NIM of China, LNE-SYRTE of France, NIST of the United States, PTB of Germany, NPL of the United Kingdom, VNIIFTRI of Russia, NICT of Japan, USNO of the United States, etc.), using satellite-based time frequency transfer technology, with Germany PTB as the comparison center, each timekeeping laboratory reports the comparison data to BIPM, BIPM gives different weights to the performance of the timekeeping clock, and the weighted average of all timekeeping clocks is obtained to obtain the free atomic time (EAL), and the international atomic time (TAI) is generated by the reference clock control, and the coordinated universal time UTC is generated after the leap second correction.

[0058] In summary, the comparison center of the coordinated universal time (UTC) is Germany PTB, the management center is France BIPM, and the standard time UTC(k) resources maintained by each country are also concentrated first. The outstanding problems are weak disaster recovery and anti-destroying ability and great security risk.

[0059] The technical scheme of the embodiment, by acquiring the calibration data of the coordinated universal time traceability module and the clock difference data of at least two timekeeping clocks in the distributed laboratory system; according to each clock difference data and its corresponding stability, determine the local atomic time; judge whether the coordinated universal time traceability module is running normally; if yes, calibrate the local atomic time according to the calibration data to obtain the standard time; otherwise, acquire the reference data of the reference clock and the prediction data of the virtual coordinated universal time prediction module and their corresponding uncertainties, and calibrate the local atomic time according to the reference data and the prediction data and their corresponding uncertainties to obtain the standard time. When the coordinated universal time is abnormal, the frequency and phase are corrected in time through the reference data and the prediction data of the coordinated universal time to determine the standard time. At the same time, each laboratory node in the distributed laboratory system has the standard time generation correction ability and time data resources, realizes the decentralization of the standard time generation, and improves the disaster recovery and safety of the standard time determination. Therefore, through the technical scheme of the present application, the problem of relatively concentrated standard time resources of each laboratory, poor disaster recovery performance and low safety is solved, and the effect of improving the disaster recovery performance and safety of the standard time generation of the distributed laboratory system is achieved.

[0060] Embodiment two

[0061] Figure 2 The flowchart of a distributed laboratory system timekeeping method provided for embodiment two of the present application, the technical scheme of the present embodiment is further refined on the basis of the above technical scheme.

[0062] Further, the "otherwise, obtaining reference reference data of the reference clock and prediction data of the virtual coordinated universal time prediction module and corresponding uncertainty thereof, and calibrating the local atomic time according to the reference reference data and the prediction data and corresponding uncertainty thereof to obtain the standard time" is refined as "otherwise, obtaining reference reference data of the reference clock and prediction data of the virtual coordinated universal time prediction module and corresponding uncertainty thereof; determining reference weights of the reference reference data and virtual weights of the prediction data according to uncertainty of the reference reference data and the prediction data; determining fused reference data by weighted averaging the reference reference data and the prediction data based on the reference weights and the virtual weights; calibrating the local atomic time by the fused reference data to obtain the standard time", so as to determine the standard time when the coordinated universal time traceability module is abnormal.

[0063] Referring to Figure 2 A time keeping method of a distributed laboratory system is shown, comprising:

[0064] S210, obtaining calibration data of a coordinated universal time traceability module and clock difference data of at least two time keeping clocks in the distributed laboratory system.

[0065] S220, determining local atomic time according to each clock difference data and corresponding stability thereof.

[0066] S230, judging whether the coordinated universal time traceability module is normally running.

[0067] S240, if yes, calibrating the local atomic time according to the calibration data to obtain the standard time.

[0068] S250, otherwise, obtaining reference reference data of the reference clock and prediction data of the virtual coordinated universal time prediction module and corresponding uncertainty thereof.

[0069] Obtaining reference reference data of the reference clock of each laboratory node and corresponding uncertainty thereof, and prediction data of the virtual coordinated universal time prediction module of the node and corresponding uncertainty thereof.

[0070] In an optional embodiment, when the reference clock is multiple, determining reference weights of each reference clock according to uncertainty of each reference clock in a preset historical time period, and obtaining reference reference data by weighted averaging each frequency data or phase data and corresponding reference weight thereof.

[0071] The reference clock is composed of not less than one international unit of second definition or second secondary definition reproduction device; when the number of reference clocks is greater than or equal to 2, the data weight of each reference clock is determined according to the uncertainty of the reference clock in the preset historical time period, that is, the reference weight. According to the frequency data or phase data of each reference clock and the corresponding reference weight, the weighted average algorithm is used to determine the standard frequency reference data, which is used to participate in the calibration of the local time when the coordinated universal time tracing module is abnormal.

[0072] S260, determining the reference weight of the reference data and the virtual weight of the prediction data according to the uncertainty of the reference data and the prediction data.

[0073] The reference weight is the weight of each reference data, and the virtual weight can be the weight of the prediction data. The uncertainty can be used to measure the accuracy of the reference data and the prediction data, therefore, according to the uncertainty of the reference data and the prediction data, the reference weight of the reference data and the virtual weight of the prediction data are determined to weight average the reference data and the prediction data, and the fusion reference data is determined.

[0074] S270, weighting and averaging the reference data and the prediction data based on the reference weight and the virtual weight to determine the fusion reference data.

[0075] Based on the reference weight and the virtual weight, the reference data and the prediction data are weighted and averaged to obtain the fusion reference data, which can be used to further calibrate the local atomic time and determine the standard time.

[0076] In an optional embodiment, after weighting and averaging the reference data and the prediction data based on the reference weight and the virtual weight to determine the fusion reference data, it further includes: after reaching a first preset time length or the virtual weight is less than a first weight threshold, using a gradual zeroing strategy to reduce the virtual weight to zero.

[0077] The first preset time length can be a pre-set time threshold for reducing the virtual weight, which can be determined by professional technicians according to experience or experiment, and the present application does not make specific limitation. The first weight threshold can be a pre-set minimum value of the virtual weight for reducing the virtual weight. The first weight threshold can be determined by professional technicians according to experience or experiment, and the present application does not make specific limitation.

[0078] The gradual zeroing strategy can be a strategy for reducing the virtual weight in stages to reduce the virtual weight to zero, and the standard time is traced to the reference clock.

[0079] The gradual zeroing strategy can be divided into the following four stages:

[0080] Preparation phase (0-30 minutes) Start the switching pre-check: Verify the reference clock group running status (require at least two reference clocks online), check the environmental parameters (temperature fluctuation <0.5℃) and confirm the storage system redundancy;

[0081] Transition phase (30-180 minutes) Execute linear weight reduction: Reduce the virtual weight by 0.01 per minute, while increasing the reference weight;

[0082] Stable phase (180-300 minutes) Complete the final switch: The virtual weight is reduced to 0, and the reference weight is fixed at 1, and the system outputs the "LOCAL_HOLDOVER" state flag;

[0083] Verification phase (after 300 minutes) Continuous monitoring for 72 hours, checking the time consistency every 15 minutes (require the time stamp deviation <1ns before switching).

[0084] Verification and monitoring system builds a closed-loop evaluation network, sets up system maintenance switch quality indicators, which integrates phase jump amplitude, frequency stability and environmental disturbance factors, and quantifies the switch quality according to 0-100 points. All verification data uses IETF RFC 3161 time stamp and SHA-3 algorithm signature, and is stored in WORM optical disk array.

[0085] Optionally, the reference clock starts the enhanced maintenance mode after the weight is zero. Specifically, the performance tracking unit collects the internal key parameters of the reference clock every 15 minutes; the environmental control unit improves the clock room temperature stability from ±0.5℃ to ±0.2℃; the power system switches to pure battery power mode (eliminates power grid interference).

[0086] Optionally, when the reference clock detects an anomaly, it is handled. Specifically, when the reference clock performance degradation (Allen variance increase >20%) is monitored, the system automatically enters the safety mode: first freeze the current reference weight for 24 hours for observation; if it is confirmed that the degradation is continuous, the reference clock weight is gradually reduced; finally generate an exception analysis report, including fault diagnosis, treatment measures and verification results.

[0087] In an optional embodiment, when the coordinated universal time traceability is restored, the system performs a phased regression process: within 48 hours, gradually reduce the reference weight, while gradually increasing the virtual weight, and switch back to coordinated universal time when the virtual weight is increased to the threshold; and start the backtracking compensation algorithm to avoid phase jump.

[0088] In an optional embodiment, full-dimension monitoring is implemented for timekeeping performance, including: short-term stability evaluation (τ = 1s to 100s) using the overlapping Allan variance method, which can be calculated in real time through FPGA hardware acceleration; long-term stability evaluation (τ = 1 day to 30 days) using the modified Allan variance (MDEV), which can generate performance trend reports every week; frequency accuracy and time deviation relative to coordinated universal time evaluation, which can be performed once a month by comparing UTC(k) with UTC. The evaluation results can be visually displayed, and the visual display includes three key indicator boards: a phase deviation and frequency accuracy trend graph showing the time difference fluctuation of the latest standard time relative to UTC; a frequency stability spectrum graph showing the variance characteristics in the range of 1s-100000s; and a device status matrix reflecting the health status of each atomic clock in real time. All evaluation data are archived in the CCTF-K001 format specified by BIPM, with a preservation period of not less than 10 years.

[0089] In an optional embodiment, timekeeping performance is managed and operated, and the key functions of operation and management include: an abnormality detection module that monitors the Allan variance of each atomic clock and reference clock in real time, automatically triggers an alarm when the threshold of 3σ is exceeded, and automatically marks the abnormal node and performs data rejection; a predictive maintenance module that generates maintenance recommendations 30 days in advance based on device runtime and environmental stress data; a configuration management module that records all parameter change history and supports state rollback at any time point. The operation strictly follows the standard procedures: routine inspection is performed daily, including device status confirmation and environmental parameter review; preventive maintenance is carried out every quarter, including atomic clock magnetic shielding inspection and optical system calibration; annual overhaul includes comprehensive performance testing and key component replacement. All operation records are stored through blockchain technology to ensure data tamper resistance.

[0090] S280, calibrate the local atomic time by fusing the reference data to obtain the standard time.

[0091] The standard time is obtained by fusing the reference data and the local atomic time and performing weighted averaging.

[0092] The technical scheme of the embodiment, by obtaining reference data of a reference clock and prediction data of a virtual coordinated universal time prediction module and their corresponding uncertainties; determining the reference weight of the reference data and the virtual weight of the prediction data according to the uncertainties of the reference data and the prediction data, to improve the reliability of the subsequent standard time; based on the reference weight and the virtual weight, weighted average of the reference data and the prediction data is performed to determine the fused reference data, in the case of abnormality of the coordinated universal time tracing module, the local atomic time is calibrated by fusing the reference data to obtain the standard time, and the accuracy of the standard time is improved.

[0093] Embodiment three

[0094] Figure 3 A flowchart of a distributed laboratory system time keeping method provided for Embodiment Three of the present application, the technical solution of the present embodiment is further refined on the basis of the above technical solution.

[0095] Further, the "generation method of prediction data of the virtual coordinated universal time prediction module" is added, specifically: "obtain historical data of the current time, initial parameter estimates of each predictor and residual sum of squares; determine initial prediction values and time complexity based on the initial parameter estimates of each predictor and the residual sum of squares, and the historical data of the current time; obtain historical data of the next time, and update the initial parameter estimates of each predictor based on the historical data to obtain updated each predictor; determine the prediction value of the next time based on the updated each predictor and the historical data of the next time", to determine the prediction data.

[0096] Referring to Figure 3 A distributed laboratory system time keeping method, comprising:

[0097] S310, obtain calibration data of a coordinated universal time tracing module and clock difference data of at least two time keeping clocks in a distributed laboratory system.

[0098] S320, determine local atomic time based on each clock difference data and its corresponding stability.

[0099] S330, judge whether the coordinated universal time tracing module is running normally.

[0100] S340, if yes, calibrate the local atomic time based on the calibration data to obtain standard time.

[0101] S350, otherwise, obtain historical data of the current time, initial parameter estimates of each predictor and residual sum of squares.

[0102] The historical data can include time difference and phase difference, and after obtaining the historical data, the historical data set can be randomly divided into several subsets. Initialize the intermediate variables of each subset, determine the initial parameter estimates of each predictor based on the intermediate variables. According to the initial parameter estimates, calculate the initial residual sum of squares of each predictor.

[0103] S360, determine the initial prediction value based on the initial parameter estimates of each predictor and the residual sum of squares, and the historical data of the current time.

[0104] Determine the weight of each predictor based on the reciprocal of the residual sum of squares of each predictor; obtain the output of each predictor based on the historical data of the current time, and obtain the initial prediction value of the next time by weighted summing the outputs of each predictor.

[0105] Specifically, the historical data is a first-in-first-out queue with length , at time, is like

[0106]

[0107] is randomly divided into subsets using non-replacement random sampling, and the length of each subset is , and the recommended values of the two satisfy:

[0108]

[0109] The predictor adopts a low-order polynomial model, which is like

[0110]

[0111] , where represents the relative frequency deviation measurement reading of the atomic clock, that is, the data in the historical data set; represents the predictor; represents the frequency drift, represents the frequency, both of which are parameters of the predictor; represents the residual sequence of the predictor.

[0112] S370, obtain the historical data at the next time of the initial prediction value, and update the initial parameter estimation of each predictor according to the historical data to obtain the updated each predictor.

[0113] Obtain the historical data at the next time, update the historical data set, find the single subset containing the oldest data in the last iteration, and replace the oldest data with the latest data; recursively update the corresponding intermediate variables for the single subset involving data update; obtain new predictor parameter estimation based on the updated intermediate variables; re-calculate the residual sum of squares of the single updated predictor from the beginning; recursively update the residual sum of squares of the remaining predictors; and assign weights based on the reciprocal of the residual sum of squares of each predictor.

[0114] In an optional embodiment, the initial parameter estimation of each predictor is obtained by least squares fitting through the inverse covariance matrix and the cumulative vector, or by least squares fitting through the closed formula.

[0115] There are two calculation methods for obtaining the predictor parameter estimation using least squares fitting, corresponding to different initialization and recursive updating methods.

[0116] Method one (using closed formula to calculate least squares fitting):

[0117] The closed formula is like

[0118]

[0119]

[0120] where denote the index of historical data, based on the closed-form formula, four intermediate variables can be set:

[0121]

[0122]

[0123]

[0124]

[0125] These intermediate variables, i.e. the minuend and the subtrahend on the numerator and the denominator in the closed-form formula, are the sample sum or the sample sum of squares of each subset. In the initialization phase, these intermediate variables of each subset are initialized first, which requires traversing each subset, and the time complexity is Then, the intermediate variables are plugged into the closed-form formula, and the initial estimate of the predictor parameters can be obtained. The time complexity can be the time complexity of traversing all historical data. After a series of predictors are obtained, the residual sum of squares and the weight corresponding to each predictor can be calculated by the following formula, which requires traversing the historical data set, and the time complexity is :

[0126]

[0127]

[0128] By weighted summing the outputs of each predictor, the predicted value at the next time can be obtained:

[0129]

[0130] In the next iteration, the oldest data in a subset i.e. is replaced by the latest data i.e. At this time, the intermediate variables corresponding to this subset can be recursively updated with the time complexity :

[0131]

[0132]

[0133]

[0134]

[0135] Based on the updated intermediate variables, the closed-form formula can be applied to obtain the updated parameter estimates of the predictors. For the updated predictors, the residual sum of squares needs to be recalculated from scratch as in the initialization phase, with time complexity . For the non-updated predictors, the residual sum of squares can be updated recursively with time complexity :

[0136]

[0137]

[0138]

[0139] After obtaining all the updated residual sums of squares, the prediction values at the next time instant can be obtained by following the weight calculation and weighted summation steps in the initialization phase, and the iteration continues.

[0140] Method Two (Using the inverse covariance matrix and the cumulative vector to calculate the least squares fit)

[0141] The inverse covariance matrix and the cumulative vector are used as intermediate variables to calculate the least squares fit as follows:

[0142]

[0143] In the initialization phase, the two intermediate variables are initialized by the following steps:

[0144] First, the observation vector and the design matrix of each subset are constructed, which requires traversing each subset with time complexity :

[0145]

[0146]

[0147] where and are the elements of each subset and their corresponding epochs, with the same meaning as in Method One. After obtaining the observation vector and the design matrix, the covariance matrix and the cumulative vector can be calculated by the following formulas:

[0148]

[0149]

[0150] It is particularly important to note that in order to ensure the stability of numerical calculation, the inverse operation in the calculation process is preceded by the regularization of the calculation results of this step, i.e. an arbitrary value of about three orders of magnitude smaller than the diagonal elements is added to the diagonal elements of the matrix to avoid a series of numerical calculation problems caused by possible ill-conditioned matrices. This step of calculation results is first regularized, i.e. an arbitrary value of about three orders of magnitude smaller than the diagonal elements is added to the diagonal elements of the matrix to avoid a series of numerical calculation problems caused by possible ill-conditioned matrices.

[0151] The covariance matrix and the cumulative vector are multiplied to obtain the initial estimate of the predictor parameter, and a series of predictors are obtained. The calculation process and method of the residual sum of squares and the weight are the same, and the time complexity is .

[0152] Similarly, in the next iteration, the oldest data in the subset is replaced by the latest data , i.e. , at this time the intermediate variable corresponding to the subset can be updated recursively with time complexity , this process is explicitly written as the following two parts of lower update and upper update: Lower update (eliminate the influence of old data)

[0153]

[0154]

[0155]

[0156]

[0157] Upper update (eliminate the influence of old data)

[0158]

[0159]

[0160]

[0161] The updated two intermediate variables are multiplied to obtain the updated parameter estimate of the predictor, and the remaining update steps are exactly the same as method one.

[0162] ​​By introducing intermediate / recursion variables, the predictor parameters and the residual sum of squares calculation results of the last iteration are passed to the next iteration, and the next iteration only needs to obtain the new prediction value by the modified update with time complexity of O(1) on this basis. Avoiding repeated traversal of the data set in each iteration, the time complexity of the algorithm is greatly reduced, and the efficiency of the algorithm is improved.

[0163] S380, determining the prediction value of the next moment of the initial prediction value according to the updated each predictor and the historical data of the next moment of the initial prediction value.

[0164] According to the updated each predictor and the historical data of the next moment of the initial prediction value, the prediction value of the next moment is obtained by weighted sum of each predictor output.

[0165] S390, obtaining the reference data of the reference clock and the prediction data of the virtual coordinated universal time prediction module and the corresponding uncertainty, and calibrating the local atomic time according to the reference data and the prediction data and the corresponding uncertainty to obtain the standard time.

[0166] The technical scheme of the embodiment, by obtaining the historical data of the current moment, the initial parameter estimation of each predictor and the residual sum of squares; according to the initial parameter estimation of each predictor and the residual sum of squares, and the historical data of the current moment, the initial prediction value is determined; obtaining the historical data of the next moment, and updating the initial parameter estimation of each predictor according to the historical data to obtain the updated each predictor; according to the updated each predictor and the historical data of the next moment, the prediction value of the next moment of the initial prediction value is determined, by introducing intermediate / recursion variables, the predictor parameters and the residual sum of squares calculation results of the last iteration are passed to the next iteration, and the prediction value of the next moment is quickly determined.

[0167] Embodiment four

[0168] Figure 4 The structure schematic diagram of a distributed laboratory system time keeping device provided by the embodiment four of the application is shown, the embodiment can be applied to the case that each laboratory node in the distributed laboratory system determines the standard time, and the distributed laboratory system time keeping device is configured in each laboratory node in the distributed laboratory system, and the specific structure of the distributed laboratory system time keeping device is as follows:

[0169] The data acquisition module 410 is used to acquire the calibration data of the coordinated universal time tracing module and the clock difference data of at least two time keeping clocks in the distributed laboratory system;

[0170] The local atomic time generation module 420 is used to determine the local atomic time according to each clock difference data and the corresponding stability;

[0171] The tracing judgment module 430 is used to judge whether the coordinated universal time tracing module is normally running;

[0172] The standard time correction module 440 is configured to calibrate the local atomic time according to the calibration data if yes, to obtain the standard time.

[0173] The standard time prediction correction module 450 is configured to obtain reference data of the reference clock and prediction data of the virtual coordinated universal time prediction module and corresponding uncertainties of the prediction data, and calibrate the local atomic time according to the reference data and the prediction data and the corresponding uncertainties, to obtain the standard time if no.

[0174] The technical scheme of the embodiment is configured to obtain calibration data of a coordinated universal time tracing module and clock difference data of at least two time-keeping clocks in a distributed laboratory system, determine the local atomic time according to each clock difference data and corresponding stabilities, determine whether the coordinated universal time tracing module is normally running, calibrate the local atomic time according to the calibration data if yes, to obtain the standard time, otherwise, obtain reference data of a reference clock and prediction data of a virtual coordinated universal time prediction module and corresponding uncertainties of the prediction data, and calibrate the local atomic time according to the reference data and the prediction data and the corresponding uncertainties, to obtain the standard time. When the coordinated universal time is abnormal, the frequency and phase are corrected in time through the reference data and the prediction data of the coordinated universal time, to determine the standard time. Meanwhile, each laboratory node in the distributed laboratory system has the standard time generation correction capability and the time data resource, the decentralization of the standard time generation is realized, and the disaster recovery performance and the safety of the standard time determination are improved. Therefore, the technical scheme solves the problems of the relatively concentrated standard time resource of each laboratory, the poor disaster recovery performance, and the low safety, and achieves the effect of improving the disaster recovery performance and the safety of the standard time generation of the distributed laboratory system.

[0175] Optionally, the local atomic time generation module 420 comprises:

[0176] The fusion weight determination unit is configured to determine the fusion weight of each clock difference data according to the corresponding stability of each clock difference data.

[0177] The local atomic time determination unit is configured to determine the local atomic time according to each clock difference data and the corresponding fusion weight.

[0178] Optionally, the standard time prediction correction module 450 comprises:

[0179] The data acquisition unit is configured to obtain reference data of a reference clock and prediction data of a virtual coordinated universal time prediction module.

[0180] The weight determination unit is configured to determine a reference weight of the reference data and a virtual weight of the prediction data according to the uncertainties of the reference data and the prediction data.

[0181] The fusion reference determination unit is configured to determine fusion reference data by weighted averaging of the reference data and the prediction data based on a reference weight and a virtual weight.

[0182] Optionally, the standard time prediction calibration module 450 further comprises:

[0183] The virtual weight reduction unit is configured to reduce the virtual weight to zero by using a gradual zeroing strategy after a first preset time length is reached or the virtual weight is less than a first weight threshold.

[0184] Optionally, the prediction data of the virtual coordinated universal time prediction module is generated based on the following manner:

[0185] The historical data of a current time, initial parameter estimates of each predictor and residual sum of squares are obtained.

[0186] The initial prediction value and the time complexity are determined according to the initial parameter estimates of each predictor and the residual sum of squares and the historical data of the current time.

[0187] The historical data of a next time is obtained, and the initial parameter estimates of each predictor are updated according to the historical data to obtain updated each predictor.

[0188] The prediction value of the next time is determined according to the updated each predictor and the historical data of the next time.

[0189] Optionally, the initial parameter estimates of each predictor are obtained by least square fitting through an inverse covariance matrix and a cumulative vector, or by least square fitting through a closed formula.

[0190] The distributed laboratory system time-keeping device provided by the embodiments of the present application can execute the distributed laboratory system time-keeping method provided by any of the embodiments of the present application, and has the corresponding function modules and beneficial effects of executing the distributed laboratory system time-keeping method.

[0191] According to the embodiments of the present application, the present application further provides an electronic device, a readable storage medium and a computer program product.

[0192] Embodiment five

[0193] Figure 5 A structural schematic diagram of an electronic device provided by the embodiment five of the present application is shown in Figure 5 The electronic device includes a processor 510, a memory 520, an input device 530 and an output device 540; the number of processors 510 in the electronic device can be one or more, Figure 5 The processor 510 in the electronic device is taken as an example; the processor 510, the memory 520, the input device 530 and the output device 540 in the electronic device can be connected through a bus or other means,Figure 5 Taking the example of a connection between China and Israel via a bus.

[0194] The memory 520, as a computer-readable storage medium, can be used to store software programs, computer-executable programs, and modules, such as the program instructions / modules corresponding to the distributed laboratory system timekeeping method in this embodiment (e.g., data acquisition module 410, local atomic time generation module 420, source tracing judgment module 430, standard time correction module 440, and standard time prediction correction module 450). The processor 510 executes various functional applications and data processing of the electronic device by running the software programs, instructions, and modules stored in the memory 520, thereby implementing the aforementioned distributed laboratory system timekeeping method.

[0195] The memory 520 may primarily include a program storage area and a data storage area. The program storage area may store the operating system and at least one application program required for a given function; the data storage area may store data created based on terminal usage. Furthermore, the memory 520 may include high-speed random access memory and may also include non-volatile memory, such as at least one disk storage device, flash memory device, or other non-volatile solid-state storage device. In some instances, the memory 520 may further include memory remotely located relative to the processor 510, which can be connected to the electronic device via a network. Examples of such networks include, but are not limited to, the Internet, intranets, local area networks, mobile communication networks, and combinations thereof.

[0196] Input device 530 can be used to receive input character information and generate key signal inputs related to user settings and function control of the electronic device. Output device 540 may include display devices such as a display screen.

[0197] Example 6

[0198] Embodiment Six of this application also provides a storage medium containing computer-executable instructions. When executed by a computer processor, the computer-executable instructions are used to execute a timekeeping method for a distributed laboratory system. The method includes: acquiring calibration data from a Coordinated Universal Time (UTC) tracing module and clock difference data and corresponding stability of at least two clocks in the distributed laboratory system; determining local atomic time based on each clock difference data and its corresponding stability; determining whether the UTC tracing module is operating normally; if so, calibrating the local atomic time based on the calibration data to obtain standard time; otherwise, acquiring reference reference data of a reference clock and prediction data and corresponding uncertainties of a virtual UTC prediction module, and calibrating the local atomic time based on the reference reference data and prediction data and their corresponding uncertainties to obtain standard time.

[0199] Of course, the storage medium provided by the embodiment of the present application includes computer executable instructions, which are not limited to the method operations described above, but can also perform related operations in the time-keeping method of the distributed laboratory system provided by any embodiment of the present application.

[0200] Through the above description of the embodiments, those skilled in the art can clearly understand that the present application can be realized by means of software and necessary general hardware, and of course can also be realized by hardware, but in many cases the former is a better embodiment. Based on such understanding, the technical solutions of the present application can be embodied in the form of a software product, and the computer software product can be stored in a computer readable storage medium, such as a floppy disk, a read-only memory (ROM), a random access memory (RAM), a FLASH, a hard disk or an optical disk, etc., including a plurality of instructions to make an electronic device (which can be a personal computer, a server, or a network device, etc.) execute the method described in each embodiment of the present application.

[0201] It is worth noting that the embodiments of the above-mentioned distributed laboratory system time-keeping device include various units and modules, which are only divided according to functional logic, but are not limited to the above-mentioned division, as long as the corresponding functions can be realized; in addition, the specific name of each functional unit is only for easy mutual distinction, and does not limit the protection scope of the present application.

[0202] Note that the above is only the preferred embodiment of the present application and the technical principle applied. Those skilled in the art will understand that the present application is not limited to the specific embodiments described herein, and those skilled in the art can make various obvious changes, readjustments and substitutions without departing from the scope of the present application. Therefore, although the present application has been described in more detail through the above embodiments, the present application is not limited to the above embodiments, and can include more other equivalent embodiments without departing from the concept of the present application, and the scope of the present application is determined by the scope of the appended claims.

Claims

1. A distributed laboratory system timekeeping method, characterized by, A laboratory node applied to a distributed laboratory system, comprising: obtaining calibration data of a coordinated universal time traceability module and clock difference data of at least two time-keeping clocks in the distributed laboratory system; determining a local atomic time according to each clock difference data and its corresponding stability; judging whether the coordinated universal time traceability module is normally running; if yes, calibrating the local atomic time according to the calibration data to obtain a standard time; otherwise, obtaining reference data of a reference clock and prediction data of a virtual coordinated universal time prediction module and their corresponding uncertainties, calibrating the local atomic time according to the reference data and the prediction data and their corresponding uncertainties to obtain the standard time.

2. The method of claim 1, wherein, The determination of the local atomic time according to each clock difference data and its corresponding stability comprises: determining a fusion weight of each clock difference data according to the corresponding stability of the clock difference data; determining the local atomic time according to each clock difference data and its corresponding fusion weight.

3. The method of claim 1, wherein, The otherwise, obtaining reference data of a reference clock and prediction data of a virtual coordinated universal time prediction module and their corresponding uncertainties, calibrating the local atomic time according to the reference data and the prediction data and their corresponding uncertainties to obtain the standard time, comprises: otherwise, obtaining reference data of a reference clock and prediction data of a virtual coordinated universal time prediction module and their corresponding uncertainties; determining a reference weight of the reference data and a virtual weight of the prediction data according to the uncertainties of the reference data and the prediction data; determining a fusion reference data by weighted averaging the reference data and the prediction data based on the reference weight and the virtual weight; calibrating the local atomic time by the fusion reference data to obtain the standard time.

4. The method of claim 3, wherein, After the fusion reference data is determined by weighted averaging the reference data and the prediction data based on the reference weight and the virtual weight, the method further comprises: after a first preset time length is reached or the virtual weight is less than a first weight threshold, adopting a gradual zeroing strategy to reduce the virtual weight to zero.

5. The method of claim 1, wherein, The prediction data of the virtual coordinated universal time prediction module is generated based on the following manner: obtaining historical data at a current time, initial parameter estimates of each predictor, and residual sum of squares; determining an initial prediction value according to the initial parameter estimates of each predictor, the residual sum of squares, and the historical data at the current time; obtaining historical data at a next time of the initial prediction value, and updating the initial parameter estimates of each predictor according to the historical data to obtain updated each predictor; determining a prediction value at the next time of the initial prediction value according to the updated each predictor and the historical data at the next time of the initial prediction value.

6. The method of claim 5, wherein, The initial parameter estimates of each predictor are obtained by least squares fitting through an inverse covariance matrix and a cumulative vector, or by least squares fitting through a closed formula.

7. A distributed laboratory system timekeeping device, characterized by, A laboratory node configured in a distributed laboratory system, comprising: The data acquisition module is configured to acquire calibration data of the coordinated universal time traceability module and clock difference data of at least two time-keeping clocks in the distributed laboratory system. The local atomic time generation module is configured to determine local atomic time according to each clock difference data and a corresponding stability thereof. The traceability judgment module is configured to judge whether the coordinated universal time traceability module is normally operated. The standard time correction module is configured to, if yes, calibrate the local atomic time according to the calibration data to obtain standard time. The standard time prediction correction module is configured to, if no, acquire reference data of a reference clock and prediction data of a virtual coordinated universal time prediction module and a corresponding uncertainty thereof, and calibrate the local atomic time according to the reference data and the prediction data and the corresponding uncertainty thereof to obtain standard time.

8. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor executes the program to implement the distributed laboratory system time-keeping method according to any one of claims 1-6.

9. A computer-readable storage medium having stored thereon a computer program, characterized in that The program is executed by the processor to implement the distributed laboratory system time-keeping method according to any one of claims 1-6.

10. A computer program product, characterised in that, The computer program product comprises a computer program which, when executed by a processor, implements the distributed laboratory system time-keeping method according to any one of claims 1-6.

Images