CN113626992A - Method and device for pre-estimating characteristic parameters of Venturi tube and computer equipment - Google Patents

Method and device for pre-estimating characteristic parameters of Venturi tube and computer equipment Download PDF

Info

Publication number
CN113626992A
CN113626992A CN202110821881.3A CN202110821881A CN113626992A CN 113626992 A CN113626992 A CN 113626992A CN 202110821881 A CN202110821881 A CN 202110821881A CN 113626992 A CN113626992 A CN 113626992A
Authority
CN
China
Prior art keywords
differential pressure
flow
water supply
cycle
venturi tube
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202110821881.3A
Other languages
Chinese (zh)
Inventor
陈永伟
索凌平
李东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China General Nuclear Power Corp
CGN Power Co Ltd
China Nuclear Power Operation Co Ltd
Original Assignee
China General Nuclear Power Corp
CGN Power Co Ltd
China Nuclear Power Operation 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 China General Nuclear Power Corp, CGN Power Co Ltd, China Nuclear Power Operation Co Ltd filed Critical China General Nuclear Power Corp
Priority to CN202110821881.3A priority Critical patent/CN113626992A/en
Publication of CN113626992A publication Critical patent/CN113626992A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

The application relates to a method and a device for pre-estimating characteristic parameters of a Venturi tube and computer equipment. The method comprises the following steps: the method comprises the steps of obtaining heat balance flow and main water supply flow of a Venturi tube, respectively carrying out flow differential pressure conversion on the heat balance flow and the main water supply flow to obtain heat balance flow differential pressure and main water supply flow differential pressure, obtaining erosion degree and operation time of the Venturi tube through the heat balance flow differential pressure and the main water supply flow differential pressure, and obtaining the erosion rate of the Venturi tube through the erosion degree and the operation time. By adopting the method, the accurate heat balance flow and main water supply flow can be obtained, the reliability of the reference standard is improved, and the erosion rate of the Venturi tube can be accurately estimated in advance through the accurate heat balance flow and main water supply flow.

Description

Method and device for pre-estimating characteristic parameters of Venturi tube and computer equipment
Technical Field
The application relates to the technical field of venturi tube characteristic parameter estimation, in particular to a method and a device for estimating venturi tube characteristic parameters and computer equipment.
Background
Nuclear power belongs to low-carbon environment-friendly clean energy, is high-density energy, has large single-machine capacity, and can effectively ensure the electric energy quality, but nuclear power units of nuclear power plants have very high technical requirements. The main water supply flow of the nuclear power unit is used as an important process parameter of the nuclear power unit and participates in water level adjustment of a steam generator of an important control system and a reactor protection logic steam-water mismatch signal.
At present, the main water supply flow of a nuclear power unit is mainly measured by adopting a venturi tube differential pressure type measuring principle. However, the venturi tube inevitably suffers from the erosion action of the fluid under the high-temperature and high-pressure operating environment, and the characteristics of the venturi tube also change, so that the measurement accuracy of the main feed water flow is affected. Therefore, how to predict the erosion rate of the venturi tube in advance to ensure the safe and stable normal operation of the nuclear power unit becomes a problem which needs to be solved urgently in the existing nuclear power plant.
Disclosure of Invention
In view of the above, it is necessary to provide a method, an apparatus and a computer device for estimating characteristic parameters of a venturi tube.
A method of estimating a characteristic parameter of a venturi, the method comprising:
acquiring heat balance flow and main water supply flow of a Venturi tube, wherein the heat balance flow is the water supply flow at an inlet of a steam generator;
respectively carrying out flow differential pressure conversion on the heat balance flow and the main water supply flow to obtain a heat balance flow differential pressure and a main water supply flow differential pressure;
acquiring the erosion degree and the operation time of the Venturi tube according to the heat balance flow differential pressure and the main feed water flow differential pressure;
and acquiring the erosion rate of the venturi tube according to the erosion degree and the running time.
In one embodiment, the performing flow differential pressure conversion on the heat balance flow and the main feed water flow respectively to obtain a heat balance flow differential pressure and a main feed water flow differential pressure includes:
and respectively carrying out flow differential pressure conversion on the heat balance flow and the main water supply flow through a preset flow threshold value to obtain the heat balance flow differential pressure and the main water supply flow differential pressure.
In one embodiment, if the operation cycle of the nuclear power plant is single cycle, the obtaining the erosion degree and the operation time of the venturi tube through the heat balance flow differential pressure and the main feed water flow differential pressure includes:
calculating a main water supply differential pressure credible value through the main water supply flow differential pressure;
obtaining the single-cycle erosion degree of the Venturi tube by making a difference between the main water supply differential pressure credible value and the heat balance flow differential pressure;
and calculating the single-cycle running time of the Venturi tube according to the interval time of the two cross comparisons, the refueling overhaul time and the low-load running time.
In one embodiment, said obtaining the erosion rate of the venturi by the erosion degree and the operating time comprises:
and carrying out quotient on the single-cycle erosion degree and the single-cycle operation time to obtain the single-cycle erosion rate of the Venturi tube.
In one embodiment, if the operation cycle of the nuclear power generating unit is multi-cycle, the obtaining the erosion rate of the venturi tube through the erosion degree and the operation time includes:
and acquiring the multi-cycle erosion rate of the venturi tube through the single-cycle erosion degree corresponding to each single cycle in the multi-cycle and the single-cycle running time corresponding to each single cycle.
In one embodiment, the obtaining the multi-cycle erosion rate of the venturi tube by the single-cycle erosion degree corresponding to each single cycle in the multi-cycle and the single-cycle running time corresponding to each single cycle comprises:
and summing the single-cycle erosion degrees corresponding to the single cycles, summing the single-cycle running times corresponding to the single cycles, and then making a quotient to obtain the multi-cycle erosion rate.
In one embodiment, the main feedwater flow differential pressure comprises: a first main water supply flow differential pressure corresponding to the upper pipe wall of the Venturi tube and a second main water supply flow differential pressure corresponding to the lower pipe wall of the Venturi tube; the main water supply differential pressure credible value is calculated through the main water supply flow differential pressure, and the method comprises the following steps:
the first main water supply flow differential pressure and the second main water supply flow differential pressure are differenced, and the absolute value of the difference value is compared with a preset rejection threshold value;
and if the absolute value of the difference is smaller than a preset rejection threshold, taking the average value of the first main water supply flow differential pressure and the second main water supply flow differential pressure as the main water supply differential pressure credible value.
An apparatus for estimating characteristic parameters of a venturi, the apparatus comprising:
the flow acquisition module is used for acquiring heat balance flow and main water supply flow of the Venturi tube, wherein the heat balance flow is the water supply flow at the inlet of the steam generator;
the flow differential pressure conversion module is used for respectively carrying out flow differential pressure conversion on the heat balance flow and the main water supply flow to obtain a heat balance flow differential pressure and a main water supply flow differential pressure;
the intermediate parameter acquisition module is used for acquiring the erosion degree and the running time of the Venturi tube according to the heat balance flow differential pressure and the main water supply flow differential pressure;
and the erosion rate acquisition module is used for acquiring the erosion rate of the Venturi tube according to the erosion degree and the running time.
A computer device comprising a memory and a processor, the memory storing a computer program, the processor implementing the following steps when executing the computer program:
acquiring heat balance flow and main water supply flow of a Venturi tube, wherein the heat balance flow is the water supply flow at an inlet of a steam generator;
respectively carrying out flow differential pressure conversion on the heat balance flow and the main water supply flow to obtain a heat balance flow differential pressure and a main water supply flow differential pressure;
acquiring the erosion degree and the operation time of the Venturi tube according to the heat balance flow differential pressure and the main feed water flow differential pressure;
and acquiring the erosion rate of the venturi tube according to the erosion degree and the running time.
A storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of:
acquiring heat balance flow and main water supply flow of a Venturi tube, wherein the heat balance flow is the water supply flow at an inlet of a steam generator;
respectively carrying out flow differential pressure conversion on the heat balance flow and the main water supply flow to obtain a heat balance flow differential pressure and a main water supply flow differential pressure;
acquiring the erosion degree and the operation time of the Venturi tube according to the heat balance flow differential pressure and the main feed water flow differential pressure;
and acquiring the erosion rate of the venturi tube according to the erosion degree and the running time.
According to the method and the device for predicting the characteristic parameters of the Venturi tube and the computer equipment, the computer equipment can obtain the heat balance flow and the main water supply flow of the Venturi tube, respectively carry out flow differential pressure conversion on the heat balance flow and the main water supply flow to obtain the heat balance flow differential pressure and the main water supply flow differential pressure, obtain the erosion degree and the running time of the Venturi tube through the heat balance flow differential pressure and the main water supply flow differential pressure, and obtain the erosion rate of the Venturi tube through the erosion degree and the running time; according to the method, the relatively accurate heat balance flow and main water supply flow can be obtained, the reliability of the reference standard is improved, and the erosion rate of the Venturi tube can be accurately estimated in advance through the accurate heat balance flow and main water supply flow.
Drawings
FIG. 1 is a diagram of an exemplary embodiment of an environment for estimating characteristic parameters of a venturi;
FIG. 2 is a diagram of a portion of an environment in which the method for estimating characteristic parameters of a venturi tube according to another embodiment may be used;
FIG. 3 is a schematic flow chart illustrating a method for estimating characteristic parameters of a venturi in one embodiment;
FIG. 4 is a schematic flow chart illustrating a method for obtaining the erosion rate and operating time of a venturi according to another embodiment;
FIG. 5 is a schematic flow chart of a method for calculating a main feed water differential pressure confidence value through a main feed water flow differential pressure in another embodiment;
FIG. 6 is a block diagram of an apparatus for estimating characteristic parameters of a venturi in one embodiment;
FIG. 7 is a diagram illustrating an internal structure of a computer device according to an embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
The estimation method of the characteristic parameters of the venturi tube can be applied to the application environment shown in fig. 1, and the estimation system of the characteristic parameters of the venturi tube comprises a steam generator, the venturi tube, a sensor and computer equipment. The steam generator can be a vertical natural circulation type saturated steam generating device, during reactor hot shut-down (nuclear power unit state) and power operation, a reactor coolant flows in a heat transfer pipe of the steam generator, heat is transferred to secondary loop water outside the heat transfer pipe, the secondary loop water naturally circulates in the steam generator, and when the secondary loop water flows through the heat transfer pipe, part of the water is changed into saturated steam. The feed water of the steam generator flows downwards to the descending channel through the inverted J-shaped pipe and enters the bottom of the ascending channel, and in the ascending channel, the two-loop water absorbs the heat transferred by the U-shaped pipe in the primary loop, the temperature of the two-loop water is continuously increased, and various forms of boiling are generated. The steam generator water supply main pipe is respectively provided with a throttling element Venturi tube which can be a main water supply flow channel.
In this embodiment, the estimation system may include three steam generators, a venturi, two sensors and a computer device, and the sensors and the computer device communicate with each other through a network. Only one steam generator and one sensor are shown in fig. 1, and fig. 2 is a schematic view of the connections between the venturi tube, the two sensors and the computer device. Optionally, the two sensors are both located on the outer wall of the venturi tube, and the two sensors are located opposite to each other, one on the upper wall and one on the lower wall. The sensor can be a water flow sensor, a differential pressure sensor, a temperature sensor, a flow velocity sensor and the like. The computer device may be, but not limited to, various personal computers, notebook computers, smart phones, tablet computers, and portable wearable devices, and may also be implemented by an independent server or a server cluster composed of a plurality of servers.
In one embodiment, as shown in fig. 3, a method for estimating characteristic parameters of a venturi is provided, which is illustrated by applying the method to the computer device in fig. 1, and comprises the following steps:
and S100, acquiring heat balance flow and main water supply flow of a Venturi tube, wherein the heat balance flow is the water supply flow at the inlet of the steam generator.
Specifically, the computer device may obtain a feed water flow at the steam generator inlet measured by the sensor, i.e., a heat balance flow QeAnd the main feed water flow of the venturi. Optionally, the main feedwater flow of the venturi may include: first main water supply flow Q measured by a sensor corresponding to the upper pipe wall of the Venturi tubeASecond main feed water flow Q measured by a sensor corresponding to the lower pipe wall of the venturi tubeB
In this embodiment, the main feed water flow rate can be approximated as the main feed waterThe real flow, and regard main water supply flow as reference standard, it is higher in order to improve the precision of measured main water supply flow, this embodiment can install high accuracy orifice plate at steam generator entrance, simultaneously, can be through installing temperature sensor and pressure sensor in venturi's the outside, measure the temperature and the pressure of main water supply in venturi through temperature sensor and pressure sensor, use standard temperature threshold value and standard pressure threshold value as the basis, constantly revise temperature and pressure, so that the heat balance flow that obtains and the precision of main water supply flow are higher in the messenger of measurement, therefore, heat balance flow QeFirst main water supply flow QAAnd a second main feed water flow rate QBMay be used as a reference standard for the present embodiment.
And S200, respectively carrying out flow differential pressure conversion on the heat balance flow and the main water supply flow to obtain a heat balance flow differential pressure and a main water supply flow differential pressure.
In particular, the computer device may balance the flow Q against heateCarrying out flow differential pressure conversion to obtain a heat balance flow differential pressure delta peAlso, the first main water supply flow Q can be respectively adjustedAAnd a second main feed water flow rate QBCarrying out flow differential pressure conversion to obtain a first main water supply flow differential pressure delta pAAnd a second main feed water flow differential pressure Δ pB. Optionally, the processing method of the flow differential pressure conversion may be to balance the heat flow QeThe arithmetic processing is performed, and the arithmetic processing may be arithmetic operation, and may also be square-open, logarithmic, power operation, or the like, and of course, may also be a combination operation of these operations.
And S300, acquiring the erosion degree and the operation time of the Venturi tube through the heat balance flow differential pressure and the main water supply flow differential pressure.
Specifically, the computer device can balance the flow differential pressure Δ p through heateFirst main feed water flow differential pressure delta pAAnd a second main feed water flow differential pressure Δ pBAnd carrying out a series of operation processing to obtain the erosion degree delta p of the Venturi tube and the running time T of the Venturi tube. Alternatively, the arithmetic operation may be an arithmetic operation, an open square, a logarithm, a power operation, or the like, of courseMay be a combination of these operations.
And S400, acquiring the erosion rate of the venturi tube through the erosion degree and the running time.
Specifically, the computer device can perform operation processing through the erosion degree Δ p' of the venturi tube and the operating time T of the venturi tube to obtain the erosion rate τ of the venturi tube. Alternatively, the arithmetic processing may be arithmetic operation, square-open operation, logarithm operation, power operation, or the like, and of course, may be combined operation of these operations.
In the method for estimating the characteristic parameters of the venturi tube, computer equipment can acquire heat balance flow and main water supply flow of the venturi tube, flow differential pressure conversion is respectively carried out on the heat balance flow and the main water supply flow to obtain heat balance flow differential pressure and main water supply flow differential pressure, the erosion degree and the running time of the venturi tube are acquired through the heat balance flow differential pressure and the main water supply flow differential pressure, and the erosion rate of the venturi tube is acquired through the erosion degree and the running time; the method can acquire relatively accurate heat balance flow and main water supply flow, improves the reliability of a reference standard, can accurately estimate the erosion rate of the Venturi tube in advance through the accurate heat balance flow and main water supply flow, further provides an accurate measurement standard for evaluating the overall operation condition and the service life of the Venturi tube, and can be used as a basis for formulating a corresponding maintenance strategy; meanwhile, the reliability of important key sensitive equipment of the nuclear power unit can be effectively improved, the safe and stable operation of the nuclear power unit is ensured, the transient state of the unit or the automatic shutdown of a reactor caused by abnormal measurement parameters in an important process is reduced, and the nuclear safety is ensured; in addition, the method can accurately estimate the erosion rate of the venturi tube only through a series of operation processing, so that the method is simple and has strong feasibility.
As an embodiment, in the step S200, the step of performing flow differential pressure conversion on the heat balance flow and the main water supply flow respectively to obtain a heat balance flow differential pressure and a main water supply flow differential pressure may be specifically implemented by the following steps: and respectively carrying out flow differential pressure conversion on the heat balance flow and the main water supply flow through a preset flow threshold value to obtain a heat balance flow differential pressure and a main water supply flow differential pressure.
Specifically, the preset flow threshold Q is setmaxThe flow threshold value can be set arbitrarily, but each nuclear power unit has a corresponding flow threshold value. Optionally, the computer device may use a preset flow threshold as a reference parameter for the heat balance flow Q respectivelyeFirst main water supply flow QAAnd a second main feed water flow rate QBCarrying out flow differential pressure conversion to obtain a heat balance flow differential pressure delta peFirst main feed water flow differential pressure delta pAAnd a second main feed water flow differential pressure Δ pB. Optionally, the processing method of the flow-to-differential pressure conversion may be a combination operation of division and square, and the combination manner may not be limited. In this embodiment, the processing method of the flow differential pressure conversion may specifically refer to the following formula:
Figure BDA0003169481310000071
Figure BDA0003169481310000072
Figure BDA0003169481310000073
further, if the operation cycle of the nuclear power plant is single cycle, as shown in fig. 4, the step of obtaining the erosion degree and the operation time of the venturi tube by using the heat balance flow differential pressure and the main feed water flow differential pressure in S300 may be implemented by the following steps:
and S310, calculating a main water supply differential pressure credible value through the main water supply flow differential pressure.
Specifically, the computer device may pass a first main feedwater flow differential pressure Δ pAAnd a second main feed water flow differential pressure Δ pBCalculating the main water supply differential pressure credible value pz. Optionally, the calculation method can be that the first main water supply flow differential pressure p is calculatedAAnd a second principal feed waterDifferential pressure p of flowBThe arithmetic operation is performed, and may be a combined arithmetic operation. Alternatively, the single cycle may be understood as a fuel cycle of a nuclear power plant.
And S320, obtaining the single-cycle erosion degree of the Venturi tube by making a difference between the main water supply differential pressure credible value and the heat balance flow differential pressure.
In particular, the computer device may be responsive to a mains supply differential pressure confidence value pzDifferential pressure Δ p with heat balance floweDifference is made to obtain the single-cycle erosion degree delta p of the Venturi tubei', i.e. Δ pi'=pz-Δpe
S330, calculating the single-cycle running time of the Venturi tube according to the interval time of the two-time cross comparison, the refueling overhaul time and the low-load running time.
It will be appreciated that during the calculation of the single-cycle runtime TiIn the (year), the nuclear power unit is in a shutdown state during the refueling overhaul period, so that the phenomenon of erosion of a venturi tube cannot be caused, the main feed water erosion only occurs during the operation of the nuclear power unit, namely the single-cycle operation time can eliminate the time consumption of the refueling overhaul. Meanwhile, because the erosion effect of the fluid is weakened in the low-load operating state, the time of the low-load operating state equivalent to the full-power operating state needs to be considered.
In this embodiment, the computer device can compare the interval T between two adjacent cross comparisonse(day), dressing change overhaul time TOLow-load operation time T of (sky) and nuclear power unitd(days), calculating the single cycle runtime T of the VenturiiIn (year), the calculation method may be a combination operation in an arithmetic operation, and the combination method is not limited. However, in this embodiment, the venturi has a single cycle runtime TiThe calculation formula of (a) is as follows:
Figure BDA0003169481310000081
in the formula, k is an equivalent coefficient. The cross-comparison may be understood as the process of comparing the measured main feedwater flow to a reference heat balance flow.
Wherein, the step of obtaining the erosion rate of the venturi tube through the erosion degree and the operating time in S400 may specifically include: and (4) carrying out quotient on the single-cycle erosion degree and the single-cycle operation time to obtain the single-cycle erosion rate of the Venturi tube.
In this embodiment, the size, characteristics, fluid density, etc. of the venturi affect the change in the main feedwater flow differential pressure Δ p, and thus, the main feedwater flow differential pressure Δ p as a function of the erosion rate τ may be expressed as:
Figure BDA0003169481310000091
in the formula (5), qmRepresenting the mass flow rate, C representing the outflow coefficient (i.e. the ratio of the actual flow rate to the theoretical flow rate), epsilon representing the expansion coefficient used to correct for the change in coefficient caused by the change in density of the fluid passing through the orifice, D representing the internal diameter of the pipe, D representing the internal diameter of the throat of the venturi, p being the density of the fluid, and tau being the erosion rate.
As can be seen, there is an approximately linear relationship between the main feedwater flow differential pressure Δ p and the erosion rate τ, and thus, in this embodiment, the computer device can calculate the single cycle erosion rate Δ pi' with Single cycle run time TiThe quotient is obtained to obtain the single-cycle erosion rate tau of the Venturi tubei. The erosion degree Δ p' can be obtained by subtracting the main water supply flow differential pressure Δ p corresponding to two moments.
The method for estimating the characteristic parameters of the Venturi tube can calculate the single-cycle erosion degree and the single-cycle running time of the Venturi tube, and further calculate the single-cycle erosion rate of the Venturi tube according to the single-cycle erosion degree and the single-cycle running time; the method can obtain relatively accurate heat balance flow and main water supply flow, so that the single-cycle erosion rate of the Venturi tube can be accurately estimated.
In some scenarios, to reduce the calculation error and improve the accuracy of the erosion rate, the present embodiment further calculates the multi-cycle erosion rate. As an embodiment, if the operation cycle of the nuclear power generating unit is multi-cycle, the step of obtaining the erosion rate of the venturi tube according to the erosion degree and the operation time in S400 may include: and acquiring the multi-cycle erosion rate of the Venturi tube through the single-cycle erosion degree corresponding to each single cycle in the multi-cycle and the single-cycle running time corresponding to each single cycle.
It can be understood that, if the operation cycle of the nuclear power plant is multi-cycle, the computer device may pass through the single-cycle erosion degree Δ p corresponding to each single cycle in the multi-cycleiAnd single cycle run time TiThe venturi's multi-cycle erosion rate τ, i is calculated to represent the running serial number of each single cycle in the multi-cycle, and if the multi-cycle includes n single cycles, i may be equal to 1, 2. Optionally, the calculation method at this time may be any one of arithmetic operations, and may also be a combined operation in the arithmetic operations, which is not limited herein.
The step of obtaining the multi-cycle erosion rate of the venturi tube through the single-cycle erosion degree corresponding to each single cycle in the multi-cycle and the single-cycle operation time corresponding to each single cycle may specifically include: and summing the single-cycle erosion degrees corresponding to the single cycles, summing the single-cycle running times corresponding to the single cycles, and then making a quotient to obtain the multi-cycle erosion rate.
In this embodiment, the computer device may map the single-cycle erosion degree Δ p corresponding to each single cycle of the multiple cyclesiSumming the single-cycle running times T corresponding to the individual single cyclesiSumming and then making a quotient to obtain the multicycle erosion rate tau, wherein the specific implementation formula can be as follows:
Figure BDA0003169481310000101
in the method for estimating the characteristic parameters of the Venturi tube, the computer equipment can further calculate the multi-cycle erosion rate by calculating the single-cycle erosion degree and the single-cycle running time, so that the accuracy of estimating the characteristic parameters of the Venturi tube is improved, an accurate measurement standard is provided for estimating the overall running condition and the service life of the Venturi tube, and the method can be used as a basis for formulating a corresponding maintenance strategy; meanwhile, the reliability of important key sensitive equipment of the nuclear power unit can be effectively improved, the safe and stable operation of the nuclear power unit is ensured, the transient state of the unit or the automatic shutdown of a reactor caused by the abnormal measurement parameters of an important process is reduced, and the nuclear safety is ensured.
As an embodiment, in order to reduce the influence of measurement errors and pipeline deviations, data with large redundant measurement deviations needs to be removed, and the reliability of such data on analysis and processing of subsequent data is relatively low, so in this embodiment, as shown in fig. 5, the step of calculating the reliability value of the main feedwater differential pressure through the main feedwater flow differential pressure in the above-mentioned S310 may be specifically implemented by the following steps:
s311, making a difference between the first main water supply flow differential pressure and the second main water supply flow differential pressure, and comparing the absolute value of the difference with a preset rejection threshold value.
In this embodiment, the data with large redundancy measurement deviation can be rejected by using the ralida criterion, the grubbs criterion, and the like. However, in this embodiment, the method for removing data with large deviation of redundant measurement may be a method for removing data with cross-comparison degradation of redundant measurement. Alternatively, the computer device may first measure the first main feedwater flow differential pressure pADifferential pressure p with second main water supply flowBAnd (4) making a difference, and comparing the absolute value of the difference with a preset rejection threshold value C (mv). Optionally, the preset rejection threshold C may be referred to as a rejection standard, the rejection standard may employ an uncertain propagation rate model and a fault probability distribution model, the technical precision of the verified sensor is calculated first, mainly considering factors such as inherent precision of the sensor, stability of the meter, temperature influence quantity, static pressure influence quantity, uncertain factors of measurement fluctuation, and the like, and then a degradation standard is calculated by the technical precision of the sensor, the degradation standard being a function of the technical precision of the flowmeter and a real-time flow measurement value.
C(mv)=fmv(QA,QBAB) (7);
Wherein f represents an uncertain propagation rate model orFailure probability distribution model, εAAnd εBRespectively representing the uncertainty of the main feed water flow of the A row and the B row of the Venturi tube measured by the sensor.
And S312, if the absolute value of the difference value is smaller than a preset rejection threshold value, taking the average value of the first main water supply flow differential pressure and the second main water supply flow differential pressure as a main water supply differential pressure credible value.
Specifically, if the computer device determines that the absolute value of the difference is less than a predetermined culling threshold (i.e. | p)A-pBWhen | ≦ C (mv), the first main feed water flow differential pressure p satisfying the condition can be setAAnd a second main feed water flow differential pressure pBMaintaining a first main feedwater flow differential pressure p that does not satisfy this conditionAAnd a second main feed water flow differential pressure pBDeleting and calculating the retained first main water supply flow differential pressure pAAnd a second main feed water flow differential pressure pBTaking the average value as a main water supply differential pressure credible value pz
According to the method for predicting the characteristic parameters of the venturi tube, data with large redundancy measurement deviation in the first main water supply flow differential pressure and the second main water supply flow differential pressure can be eliminated, so that a main water supply differential pressure credible value with higher accuracy can be obtained, and the accuracy of predicting the characteristic parameters of the venturi tube can be improved.
It should be understood that although the various steps in the flow charts of fig. 3-5 are shown in order as indicated by the arrows, the steps are not necessarily performed in order as indicated by the arrows. The steps are not performed in the exact order shown and described, and may be performed in other orders, unless explicitly stated otherwise. Moreover, at least some of the steps in fig. 3-5 may include multiple steps or multiple stages, which are not necessarily performed at the same time, but may be performed at different times, which are not necessarily performed in sequence, but may be performed in turn or alternately with other steps or at least some of the other steps.
In one embodiment, as shown in fig. 6, there is provided an estimation device for characteristic parameters of a venturi tube, comprising: flow acquisition module 11, flow differential pressure conversion module 12, intermediate parameter acquisition module 13 and erosion rate acquisition module 14, wherein:
the flow acquiring module 11 is used for acquiring a thermal balance flow and a main water supply flow of the venturi tube, wherein the thermal balance flow is the water supply flow at an inlet of the steam generator;
the flow differential pressure conversion module 12 is used for respectively performing flow differential pressure conversion on the heat balance flow and the main water supply flow to obtain a heat balance flow differential pressure and a main water supply flow differential pressure;
the intermediate parameter acquisition module 13 is used for acquiring the erosion degree and the operation time of the venturi tube through the heat balance flow differential pressure and the main water supply flow differential pressure;
and the erosion rate acquisition module 14 is used for acquiring the erosion rate of the Venturi tube according to the erosion degree and the running time.
The estimation device for the characteristic parameters of the venturi tube provided by the embodiment can implement the method embodiments, and the implementation principle and the technical effect are similar, and are not described herein again.
In one embodiment, the flow differential pressure conversion module 12 is specifically configured to perform flow differential pressure conversion on the heat balance flow and the main water supply flow through a preset flow threshold, so as to obtain a heat balance flow differential pressure and a main water supply flow differential pressure.
The estimation device for the characteristic parameters of the venturi tube provided by the embodiment can implement the method embodiments, and the implementation principle and the technical effect are similar, and are not described herein again.
In one embodiment, if the operation period of the nuclear power plant is a single cycle, the intermediate parameter obtaining module 13 includes: a confidence value calculation unit, a single-cycle aggressiveness calculation unit, and a single-cycle runtime calculation unit, wherein,
the credibility value calculating unit is used for calculating a main water supply differential pressure credibility value through a main water supply flow differential pressure;
the single-cycle erosion degree calculation unit is used for obtaining the single-cycle erosion degree of the Venturi tube by making a difference between the main water supply differential pressure credible value and the heat balance flow differential pressure;
and the single-cycle running time calculation unit is used for calculating the single-cycle running time of the Venturi tube according to the interval time of the two cross comparisons, the refueling overhaul time and the low-load running time.
The estimation device for the characteristic parameters of the venturi tube provided by the embodiment can implement the method embodiments, and the implementation principle and the technical effect are similar, and are not described herein again.
In one embodiment, the erosion rate acquisition module 14 comprises a single-cycle erosion rate acquisition unit, wherein,
and the single-cycle erosion rate acquisition unit is used for making a quotient between the single-cycle erosion degree and the single-cycle operation time to obtain the single-cycle erosion rate of the Venturi tube.
The estimation device for the characteristic parameters of the venturi tube provided by the embodiment can implement the method embodiments, and the implementation principle and the technical effect are similar, and are not described herein again.
In one embodiment, if the operation cycle of the nuclear power plant is multi-cycle, the erosion rate obtaining module 14 further includes: a multi-cycle erosion rate acquisition unit, wherein,
and the multi-cycle erosion rate acquisition unit is used for acquiring the multi-cycle erosion rate of the Venturi tube through the single-cycle erosion degree corresponding to each single cycle in the multi-cycle and the single-cycle running time corresponding to each single cycle.
The estimation device for the characteristic parameters of the venturi tube provided by the embodiment can implement the method embodiments, and the implementation principle and the technical effect are similar, and are not described herein again.
In one embodiment, the multi-cycle erosion rate obtaining unit is specifically configured to sum the single-cycle erosion degrees corresponding to the single cycles, sum the single-cycle running times corresponding to the single cycles, and then perform quotient calculation to obtain the multi-cycle erosion rate.
The estimation device for the characteristic parameters of the venturi tube provided by the embodiment can implement the method embodiments, and the implementation principle and the technical effect are similar, and are not described herein again.
In one embodiment, the main feedwater flow differential pressure comprises: a first main water supply flow differential pressure corresponding to the upper pipe wall of the Venturi tube and a second main water supply flow differential pressure corresponding to the lower pipe wall of the Venturi tube; the credibility value calculation unit comprises: a first operation subunit and a second operation subunit, wherein,
the first operation subunit is used for making a difference between the first main water supply flow differential pressure and the second main water supply flow differential pressure and comparing the absolute value of the difference with a preset rejection threshold value;
and the second operation subunit is used for taking the average value of the first main water supply flow differential pressure and the second main water supply flow differential pressure as a main water supply differential pressure credible value when the absolute value of the difference value is smaller than a preset rejection threshold value.
The estimation device for the characteristic parameters of the venturi tube provided by the embodiment can implement the method embodiments, and the implementation principle and the technical effect are similar, and are not described herein again.
The specific definition of the estimation device for the characteristic parameters of the venturi tube can be referred to the definition of the estimation method for the characteristic parameters of the venturi tube in the above, and will not be described in detail herein. The modules in the estimation device of the characteristic parameters of the venturi tube can be wholly or partially realized by software, hardware and a combination thereof. The modules can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the modules.
In one embodiment, a computer device is provided, which may be a server, the internal structure of which may be as shown in fig. 7. The computer device includes a processor, a memory, and a network interface connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program, and a database. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The database of the computer device is used for storing the heat balance flow and the main water supply flow. The network interface of the computer device is used for communicating with an external terminal through a network connection. The computer program is executed by a processor to implement a method of estimating a characteristic parameter of a venturi tube.
Those skilled in the art will appreciate that the architecture shown in fig. 7 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In one embodiment, a computer device is provided, comprising a memory and a processor, the memory having a computer program stored therein, the processor implementing the following steps when executing the computer program:
acquiring heat balance flow and main water supply flow of a Venturi tube, wherein the heat balance flow is the water supply flow at an inlet of a steam generator;
respectively carrying out flow differential pressure conversion on the heat balance flow and the main water supply flow to obtain a heat balance flow differential pressure and a main water supply flow differential pressure;
acquiring the erosion degree and the operation time of the Venturi tube through the heat balance flow differential pressure and the main water supply flow differential pressure;
and acquiring the erosion rate of the Venturi tube according to the erosion degree and the running time.
In one embodiment, a storage medium is provided having a computer program stored thereon, the computer program when executed by a processor implementing the steps of:
acquiring heat balance flow and main water supply flow of a Venturi tube, wherein the heat balance flow is the water supply flow at an inlet of a steam generator;
respectively carrying out flow differential pressure conversion on the heat balance flow and the main water supply flow to obtain a heat balance flow differential pressure and a main water supply flow differential pressure;
acquiring the erosion degree and the operation time of the Venturi tube through the heat balance flow differential pressure and the main water supply flow differential pressure;
and acquiring the erosion rate of the Venturi tube according to the erosion degree and the running time.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database or other medium used in the embodiments provided herein can include at least one of non-volatile and volatile memory. Non-volatile Memory may include Read-Only Memory (ROM), magnetic tape, floppy disk, flash Memory, optical storage, or the like. Volatile Memory can include Random Access Memory (RAM) or external cache Memory. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM), among others.
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above-mentioned embodiments only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.

Claims (10)

1. A method for estimating characteristic parameters of a Venturi tube is characterized by comprising the following steps:
acquiring heat balance flow and main water supply flow of a Venturi tube, wherein the heat balance flow is the water supply flow at an inlet of a steam generator;
respectively carrying out flow differential pressure conversion on the heat balance flow and the main water supply flow to obtain a heat balance flow differential pressure and a main water supply flow differential pressure;
acquiring the erosion degree and the operation time of the Venturi tube according to the heat balance flow differential pressure and the main feed water flow differential pressure;
and acquiring the erosion rate of the venturi tube according to the erosion degree and the running time.
2. The method of claim 1, wherein the converting the thermal balance flow and the main feed water flow into a differential flow pressure to obtain a differential thermal balance flow and a differential main feed water flow comprises:
and respectively carrying out flow differential pressure conversion on the heat balance flow and the main water supply flow through a preset flow threshold value to obtain the heat balance flow differential pressure and the main water supply flow differential pressure.
3. The method of claim 1, wherein obtaining the erosion rate and the operating time of the venturi tube through the thermal balance flow differential pressure and the main feedwater flow differential pressure if the operating cycle of the nuclear power plant is single cycle comprises:
calculating a main water supply differential pressure credible value through the main water supply flow differential pressure;
obtaining the single-cycle erosion degree of the Venturi tube by making a difference between the main water supply differential pressure credible value and the heat balance flow differential pressure;
and calculating the single-cycle running time of the Venturi tube according to the interval time of the two cross comparisons, the refueling overhaul time and the low-load running time.
4. The method of claim 3, wherein said obtaining the venturi's erosion rate from said erosion level and said runtime comprises:
and carrying out quotient on the single-cycle erosion degree and the single-cycle operation time to obtain the single-cycle erosion rate of the Venturi tube.
5. The method of claim 4, wherein obtaining the erosion rate of the venturi tube from the erosion rate and the operating time if the operating cycle of the nuclear power plant is multi-cycle comprises:
and acquiring the multi-cycle erosion rate of the venturi tube through the single-cycle erosion degree corresponding to each single cycle in the multi-cycle and the single-cycle running time corresponding to each single cycle.
6. The method of claim 5, wherein obtaining the multi-cycle erosion rate of the venturi by a single-cycle erosion level for each single cycle of the multi-cycle and a single-cycle run time for each single cycle comprises:
and summing the single-cycle erosion degrees corresponding to the single cycles, summing the single-cycle running times corresponding to the single cycles, and then making a quotient to obtain the multi-cycle erosion rate.
7. The method of claim 3, wherein the main feedwater flow differential pressure comprises: a first main water supply flow differential pressure corresponding to the upper pipe wall of the Venturi tube and a second main water supply flow differential pressure corresponding to the lower pipe wall of the Venturi tube; the main water supply differential pressure credible value is calculated through the main water supply flow differential pressure, and the method comprises the following steps:
the first main water supply flow differential pressure and the second main water supply flow differential pressure are differenced, and the absolute value of the difference value is compared with a preset rejection threshold value;
and if the absolute value of the difference is smaller than a preset rejection threshold, taking the average value of the first main water supply flow differential pressure and the second main water supply flow differential pressure as the main water supply differential pressure credible value.
8. An estimation device of characteristic parameters of a venturi tube, characterized in that the device comprises:
the flow acquisition module is used for acquiring heat balance flow and main water supply flow of the Venturi tube, wherein the heat balance flow is the water supply flow at the inlet of the steam generator;
the flow differential pressure conversion module is used for respectively carrying out flow differential pressure conversion on the heat balance flow and the main water supply flow to obtain a heat balance flow differential pressure and a main water supply flow differential pressure;
the intermediate parameter acquisition module is used for acquiring the erosion degree and the running time of the Venturi tube according to the heat balance flow differential pressure and the main water supply flow differential pressure;
and the erosion rate acquisition module is used for acquiring the erosion rate of the Venturi tube according to the erosion degree and the running time.
9. A computer device comprising a memory and a processor, the memory storing a computer program, wherein the processor implements the steps of the method of any one of claims 1 to 7 when executing the computer program.
10. A storage medium having a computer program stored thereon, the computer program, when being executed by a processor, realizing the steps of the method according to any one of claims 1 to 7.
CN202110821881.3A 2021-07-19 2021-07-19 Method and device for pre-estimating characteristic parameters of Venturi tube and computer equipment Pending CN113626992A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110821881.3A CN113626992A (en) 2021-07-19 2021-07-19 Method and device for pre-estimating characteristic parameters of Venturi tube and computer equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110821881.3A CN113626992A (en) 2021-07-19 2021-07-19 Method and device for pre-estimating characteristic parameters of Venturi tube and computer equipment

Publications (1)

Publication Number Publication Date
CN113626992A true CN113626992A (en) 2021-11-09

Family

ID=78380379

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110821881.3A Pending CN113626992A (en) 2021-07-19 2021-07-19 Method and device for pre-estimating characteristic parameters of Venturi tube and computer equipment

Country Status (1)

Country Link
CN (1) CN113626992A (en)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201034647Y (en) * 2007-02-01 2008-03-12 长春锅炉仪表程控设备股份有限公司 Combined type damp steam flow quantity, dryness measurement mechanism
CN102175561A (en) * 2011-01-21 2011-09-07 大连海事大学 Cavitation erosion test equipment for testing material performance and test method thereof
CN205879978U (en) * 2016-07-14 2017-01-11 南京擎能自动化设备有限公司 Venturi type wind speed differential pressure measuring device
CN109192343A (en) * 2018-07-11 2019-01-11 岭澳核电有限公司 Reduce the measured deviation method and device of pressurized water reactor out-pile nuclear measurement system
CN109192338A (en) * 2018-07-11 2019-01-11 岭澳核电有限公司 A kind of scaling method of nuclear power unit start-up course centre range protection definite value
CN109215820A (en) * 2018-08-14 2019-01-15 中广核核电运营有限公司 Nuclear power plant evaporator Steam-water Flow Measurement channel comparative approach and system
CN109215821A (en) * 2018-08-14 2019-01-15 中广核核电运营有限公司 Nuclear power plant evaporator main feedwater flow redundant measurement channel comparative approach and system
CN210671309U (en) * 2019-09-10 2020-06-05 航天(北京)食品技术研究院 Liquid manure circulating device
CN111928925A (en) * 2020-07-16 2020-11-13 中广核核电运营有限公司 Venturi tube fault correction method and device, computer equipment and storage medium
CN111951988A (en) * 2020-07-13 2020-11-17 中广核核电运营有限公司 Method and device for detecting abnormity of main water supply flow and computer equipment
CN212871319U (en) * 2020-04-28 2021-04-02 中广核核电运营有限公司 Venturi tube for nuclear power station water supply system
CN112964316A (en) * 2021-02-01 2021-06-15 深圳市联恒星科技有限公司 Long-throat-neck venturi tube-based moisture measurement and uncertainty prediction method and system

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201034647Y (en) * 2007-02-01 2008-03-12 长春锅炉仪表程控设备股份有限公司 Combined type damp steam flow quantity, dryness measurement mechanism
CN102175561A (en) * 2011-01-21 2011-09-07 大连海事大学 Cavitation erosion test equipment for testing material performance and test method thereof
CN205879978U (en) * 2016-07-14 2017-01-11 南京擎能自动化设备有限公司 Venturi type wind speed differential pressure measuring device
CN109192343A (en) * 2018-07-11 2019-01-11 岭澳核电有限公司 Reduce the measured deviation method and device of pressurized water reactor out-pile nuclear measurement system
CN109192338A (en) * 2018-07-11 2019-01-11 岭澳核电有限公司 A kind of scaling method of nuclear power unit start-up course centre range protection definite value
CN109215820A (en) * 2018-08-14 2019-01-15 中广核核电运营有限公司 Nuclear power plant evaporator Steam-water Flow Measurement channel comparative approach and system
CN109215821A (en) * 2018-08-14 2019-01-15 中广核核电运营有限公司 Nuclear power plant evaporator main feedwater flow redundant measurement channel comparative approach and system
CN210671309U (en) * 2019-09-10 2020-06-05 航天(北京)食品技术研究院 Liquid manure circulating device
CN212871319U (en) * 2020-04-28 2021-04-02 中广核核电运营有限公司 Venturi tube for nuclear power station water supply system
CN111951988A (en) * 2020-07-13 2020-11-17 中广核核电运营有限公司 Method and device for detecting abnormity of main water supply flow and computer equipment
CN111928925A (en) * 2020-07-16 2020-11-13 中广核核电运营有限公司 Venturi tube fault correction method and device, computer equipment and storage medium
CN112964316A (en) * 2021-02-01 2021-06-15 深圳市联恒星科技有限公司 Long-throat-neck venturi tube-based moisture measurement and uncertainty prediction method and system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CHEN YONGWEI, XIE YONGJING: "\"NPP main feedwater flow Venturi tube erosion characteristics study and improvement\"", 《ANNALS OF NUCLEAR ENERGY》, pages 828 - 834 *
YONGWEI CHENA, YONGGANG LI, ET AL.: ""A multi-level assessment and correction method for Venturi tube flow measurements"", 《NUCLEAR ENGINEERING AND DESIGN》, pages 1 - 7 *
陈永伟、付敬强等: ""核电站蒸汽发生器水位控制主/旁路阀切换改进与优化"", 《原子能科学技术》, vol. 49, no. 3, pages 523 - 528 *

Similar Documents

Publication Publication Date Title
CN111951988B (en) Method and device for detecting abnormity of main water supply flow and computer equipment
CN104536290A (en) Soft measuring method and system based on kernel principal component analysis and radial basis function neural network
CN108446465B (en) Method for measuring and calculating steam quantity for thermal power plant on line through working medium decomposition
CN111982245B (en) Main steam flow calibration method, device, computer equipment and storage medium
WO2022067308A3 (en) Two and three-dimensional model based correction of elbow tap flow measurement
Dong et al. Online reliability assessment of energy systems based on a high-order extended-state-observer with application to nuclear reactors
CN109709911B (en) On-line measuring method and system for leakage of circulating working medium of thermal power generating unit
RU2621422C2 (en) System and method for steam turbine parameter testing
CN112037946B (en) Leakage detection method and device for steam generator of nuclear power station
CN113626992A (en) Method and device for pre-estimating characteristic parameters of Venturi tube and computer equipment
CN112632873B (en) Method for determining position arrangement mode of measuring points of temperature thermocouple of nuclear reactor hot leg
CN111928920B (en) Evaporator liquid level detection method, evaporator liquid level detection device, computer equipment and storage medium
CN114937512A (en) Method and system for flow compensation of coolant of nuclear power unit primary loop
CN114662281A (en) Reactor on-line protection method and system and readable storage medium
CN112163323A (en) State estimation method and system of comprehensive energy system
Szega Using data reconciliation to improve the reliability of the energy evaluation of a gas-and-steam CHP unit
CN113434803B (en) Thermal power calculation method, device, medium and equipment of thermal power conversion system
Zhang et al. Development of online validation and monitoring system for the thermal performance of nuclear power plant in service
Webb et al. Comparison of CECOR algorithm to Lagrange multiplier method to estimate reactor power distributions
CN117709760A (en) Method and platform for supervising thermodynamic performance of unit based on data quality evaluation and coordination
Quiñones-Grueiro et al. Leaks' detection in water distribution networks with demand patterns
Szega Issues of an optimization of measurements location in redundant measurements systems of an energy conversion process–A case study
CN117494453A (en) Systematic parameter estimation method for multi-reactor one-machine nuclear power station
Wang et al. A sensitivity analysis for the adequacy assessment of a multi-state physics modeling approach for reliability analysis
Liu et al. Methods for Predicting the Minimum Temperature of the Outage Loop and the Maximum Power Caused by the Low‐Temperature Coolant

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