CN116702392A - Rainwater pipe network simulation method based on design parameter equivalent setting - Google Patents

Rainwater pipe network simulation method based on design parameter equivalent setting Download PDF

Info

Publication number
CN116702392A
CN116702392A CN202310549990.3A CN202310549990A CN116702392A CN 116702392 A CN116702392 A CN 116702392A CN 202310549990 A CN202310549990 A CN 202310549990A CN 116702392 A CN116702392 A CN 116702392A
Authority
CN
China
Prior art keywords
swmm
pipe network
rainfall
runoff
model
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
CN202310549990.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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN202310549990.3A priority Critical patent/CN116702392A/en
Publication of CN116702392A publication Critical patent/CN116702392A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/14—Pipes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Algebra (AREA)
  • Data Mining & Analysis (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a rainwater pipe network simulation method based on design parameter equivalent setting, which consists of a surface production confluence model and a one-dimensional pipe network confluence model based on the design parameter equivalent setting. Firstly, calculating the surface runoff production process by taking runoff coefficient equivalence as a target and adopting a Hoton infiltration model by combining a phi exponential method; secondly, calculating the surface confluence process by adopting a constant flow time line method, and determining the runoff process line of each sub-catchment area; and finally, coupling the surface production confluence model and the SWMM hydrodynamic module is completed through node inflow attribute setting, and the SWMM dynamic link library is called to drive simulation calculation, so that a simulation result is output. The method solves the problem of equivalence of the design conditions of the model technology when the model technology is applied in the planning and design stage of the rainwater pipe network, avoids design deviation caused by undefined model parameters, and ensures the reliability of the model in the planning and design application of the rainwater pipe network.

Description

Rainwater pipe network simulation method based on design parameter equivalent setting
Technical Field
The invention relates to a rainwater pipe network simulation method based on design parameter equivalent setting, and belongs to the technical intersection field of urban rainwater pipe networks and hydraulic models.
Background
At present, a constant uniform reasoning formula is adopted to design and calculate a rainwater pipe network, and when the method is applied to large-scale rainwater pipe network design, deviation generated by non-uniform flow state of actual water flow in a pipeline can be gradually accumulated, so that the design result is influenced. Therefore, the outdoor drainage design Standard (GB 50014-2021) requires that when the catchment area of a rainwater network exceeds 2km 2 And determining the design flow of the rainwater by adopting a mathematical model method.
With the popularization of model technology, a lot of hydraulic model software is widely applied in China, such as EPASWMM software developed by the United states environmental protection agency, inforworks series software of Wallingford company in England, MIKE series software of Denmark DHI company, and the like, wherein SWMM model is widely applied to the management of urban drainage systems due to the open source characteristic. The SWMM model calculates the surface production flow and the confluence process by adopting a primary damage after-damage method and a nonlinear reservoir respectively, and the required parameters are numerous and only provide a reasonable value range. In the planning and design stage of the rainwater pipe network in China, the runoff coefficient is generally adopted to reflect the water loss in the surface production and confluence process, actual measurement data are not available, accurate construction of a model is difficult to achieve, and blind application is difficult to guide the planning and design of the rainwater pipe network. Therefore, in order to ensure the reliability of the model technology in the rainwater pipe network planning design application, the core lies in that the setting of the model parameters is consistent or equivalent to the pipe network design parameters.
Disclosure of Invention
Aiming at the reliability problem of the current model technology in the rainwater pipe network planning and design, the invention provides a rainwater pipe network simulation method based on the equivalent setting of design parameters, ensures the reliability of the model in the rainwater pipe network planning and design application, and provides technical support for the application of the model technology in the rainwater pipe network planning and design.
In order to achieve the above object, the present invention adopts the following technical scheme:
a rainwater pipe network simulation method based on design parameter equivalent setting comprises the following steps:
s1: basic data of a rainwater pipe network design scheme are obtained, and setting parameters of a model, the attribute and the space information of the pipe network are written into the SWMM. Inp file according to the SWMM. Inp file format requirement;
s2: collecting runoff coefficients, water collecting areas, water collecting time and design rainfall process lines of all sub-water collecting areas of a pipe network, and calculating the runoff process lines of all the sub-water collecting areas based on a surface production confluence model which is set by design parameters in an equivalent mode;
s3: according to the corresponding relation between the sub-catchment area and the inspection well, inflow attribute information of the inspection well is set, the coupling of the surface production confluence model and the pipe network confluence model is completed, and inflow information of the inspection well and a runoff process line of the sub-catchment area are written into swmm.inp according to the format requirement of the inp file, so that the manufacturing of swmm.inp file is completed;
s4: utilizing the prepared swmm.inp file to call swmm_open, swmm_run, swmm_close and other functions in SWMM.dll, simulating a rainwater pipe network design scheme, and outputting swmm.rpt and swmm.out files;
s5: and calling OpenSwmmOutFile and GetGlobalValue functions in SWMM. Dll to open swmm.out files, and then completing the reading of simulation results by using functions such as SWgetnodeindex, SWgetlinkindex and GetSwmmResult.
Further, the specific operation of S1 is:
1) Determining the total duration of a simulation process according to a design rainfall process line, and writing model setting parameters into a swmm.inp file according to the format requirement of the [ OPTIONS ] part of contents in the inp file;
2) According to the attribute information and space coordinate information of the pipe network, writing the format requirements of partial contents of [ JUNCTINS ], [ OUTFALLS ], [ CONDUITS ], [ XSECTIONS ], [ COORDINATES ] and [ VERTICES ] in the swmm.inp file in sequence;
further, the specific operation of S2 is:
1) The design rainfall process line is taken as input, and the flow and rainfall runoff coefficient equivalence is taken as target (formula (1) and formula (2)), based onAn exponential method (formula (3)), and determining relevant parameters of a Huoton infiltration model (formula (4)) in a trial calculation mode, and calculating a net rain process line;
f i ＝f c +(f 0 -f c )e -ki (4)
wherein: i and j are rainfall periods, i=1, 2, …, n, the invention takes 15min as peak flow duration, so j+14 is set; r is R i Net rainfall (mm) for the i-th minute; p (P) i Rainfall (mm) for the i-th minute; c is a rainfall runoff coefficient; c (C) max Is the flow runoff coefficient; f (f) 0 Is the initial infiltration rate (mm/min); f (f) c To stabilize the infiltration rate (mm/min); k is the attenuation coefficient; f (f) i The infiltration intensity (mm) was measured at the i-th period.
(1) Flow runoff coefficient equivalence criterion: the average runoff coefficient of the net rainfall of 15min maximum in the net rainfall process line is equal to the flow runoff coefficient;
(2) rainfall runoff coefficient equivalence criterion: the ratio of the total net rainfall to the total rainfall of the whole rainfall is equal to the rainfall runoff coefficient.
2) Setting a linear confluence curve type by taking the water collecting time, the water collecting area and the rain purifying process line of each sub-catchment area as inputs, and calculating the runoff process line of each sub-catchment area by using a constant flow time line method according to a formula (5);
wherein: i. j is the rainfall moment; t is t d Collecting water for the sub-catchment area (min); q (Q) i Flow rate (L/s) at time i; r is R i-j The net rain strength (mm) at the moment i-j is R when i-j is less than or equal to 0 i-j =0; f is the area of regional confluence (m) 2 ) The method comprises the steps of carrying out a first treatment on the surface of the Alpha is a unit conversion coefficient, and the value is 0.167.
Further, the specific operation of S3 is:
1) Numbering the runoff process lines of the corresponding sub-catchment areas according to the ID of the inspection well according to the corresponding relation between each sub-catchment area and the inspection well, such as InlowID;
2) According to the requirement of the inp file on the part of [ INFLOWS ], writing all inspection well IDs and corresponding runoff process line IDs into swmm. Inp files, setting 'Constient' and 'Type' as 'FLOW' attributes, and finishing the coupling setting of the model;
3) And writing the runoff process lines of all the sub-catchment areas into the swmm.inp file according to the format requirement of the inp file on the [ TIMESERIES ] part of content, and completing the manufacture of the swmm.inp file.
The specific requirements of the format of the inp file should meet the relevant specifications in STORM WATER MANAGEMENT MODEL USER' S MUNUALVERSION 5.1.
The beneficial effects of the invention are mainly as follows:
1. the invention takes 15min as the peak flow duration calculated by the flow runoff coefficient, adopts the runoff coefficient and the flow runoff coefficientThe method is characterized in that the surface flow process is calculated in an exponential method combined mode, and the surface confluence process is calculated by adopting an equal flow time line method, so that the equivalence of design conditions can be better realized, and the reliability of the model in the planning and design application of a rainwater pipe network is ensured under the condition that all flow confluence parameters cannot be acquired through coupling of the surface flow confluence model and the SWMM hydrodynamic module。
2. The method can be applied to the verification of the design scheme of the rainwater pipe network, and the original planning and design scheme can be adjusted according to the combination of model verification, so that the design precision of the rainwater pipe network is improved, and technical support is provided for engineering designers.
Drawings
FIG. 1 is a schematic flow chart of the present invention;
FIG. 2 is a schematic diagram of the coupling of the surface production confluence model and the one-dimensional pipe network confluence model according to the present invention;
FIG. 3 is a schematic plan view of a rainwater pipe network system design;
the model-coupled [ INFLOWS ] content in the swmm.inp file created in the example of FIG. 4;
FIG. 5 is a calculation of the surface production flow of a sub-catchment area in an example;
FIG. 6 is a calculation of the surface confluence of a sub-catchment area in an example;
FIG. 7 is a pipe fullness profile of an example pipe network design.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the following details are set forth in conjunction with the Beijing example, as shown in FIG. 1, and the specific steps are as follows:
taking a rainwater pipe network in a certain Beijing area as an example, as shown in fig. 3, a complete diversion system is adopted. The total area of catchment in the drainage areas of the pipe network sheet areas is 2.48km 2 The total of 37 sub-catchment areas, the pipe network system comprises 38 pipe sections, the total length of 10.41km,38 inspection wells and 1 water outlet, the design standard adopts 3 years to meet, the flow runoff coefficient of each sub-catchment area is 0.60-0.80, and the rainfall runoff coefficient is correspondingly 0.48-0.68. The rainfall process line is calculated by the Beijing city hydrologic total station based on a II region heavy rain intensity formula, as shown in a formula (6), by adopting the same frequency amplification method, the rainfall duration is 180min, and the time step is 1min.
By collecting water from each subThe regional flow runoff coefficient and the rainfall runoff coefficient are equivalent as targets, and the method adoptsAnd (3) an exponential method, fitting and determining the parameters of a Huton hypotonic equation, and calculating a net rain process line of a water collecting area of the pipe network drainage basin of the sheet area. And (3) using an equal flow time line method, taking a rain purification process line as input, simulating a surface confluence process, and completing calculation of runoff process lines of all sub-catchment areas. The water collecting time in the pipe network is 15min, and the area of the water collecting area is 15.61hm 2 For example, the runoff peak flow of the sub-catchment area calculated by the method is 2.569m 3 And/s, compared with the result of the reasoning formula, the relative error is-1.28%, and the equivalent setting of the design parameters can be better realized. The calculation process is shown in fig. 5 and 6.
After the calculation of runoff process lines of all sub-catchment areas is completed, as shown in fig. 2, the coupling of the surface production confluence model and the pipe network confluence model is completed by applying the method of the invention, as shown in fig. 4, simulation analysis is carried out on the design scheme of the rainwater pipe network of the embodiment, and the maximum fullness simulation results of all pipelines are counted, as shown in fig. 7. Through the simulation analysis of the application model to the design scheme, whether all pipe sections reach the design standard is determined, wherein the pipe sections with the fullness of 1 need to be adjusted by a designer so as to meet the design standard, and the design reliability of the whole pipe network is improved.

Claims (5)

1. The rainwater pipe network simulation method based on design parameter equivalent setting is characterized by comprising the following steps of:
s1: basic data of a rainwater pipe network design scheme are obtained, and setting parameters of a model, the attribute and the space information of the pipe network are written into the SWMM. Inp file according to the SWMM. Inp file format requirement;
s2: collecting runoff coefficients, water collecting areas, water collecting time and design rainfall process lines of all sub-water collecting areas of a pipe network, and calculating the runoff process lines of all the sub-water collecting areas based on a surface production confluence model which is set by design parameters in an equivalent mode;
s3: according to the corresponding relation between the sub-catchment area and the inspection well, inflow attribute information of the inspection well is set, the coupling of the surface production confluence model and the pipe network confluence model is completed, and inflow information of the inspection well and a runoff process line of the sub-catchment area are written into swmm.inp according to the format requirement of the inp file, so that the manufacturing of swmm.inp file is completed;
s4: the function in SWMM. Dll is called by using the prepared swmm. Inp file, simulation of the design scheme of the rainwater pipe network is carried out, and swmm. Rpt and swmm. Out files are output;
s5: and calling OpenSwmmOutFile and GetGlobalValue functions in SWMM. Dll to open swmm. Out files, and then completing reading of simulation results by using the functions.
2. The rainwater pipe network simulation method based on design parameter equivalent setting according to claim 1, wherein S1 comprises the following specific substeps:
1) Determining the total duration of a simulation process according to a design rainfall process line, and writing model setting parameters into a swmm.inp file according to the format requirement of the [ OPTIONS ] part of contents in the inp file;
2) According to the attribute information and space coordinate information of pipe network, according to the format requirements of the inp file pair [ JUNCTINS ], [ OUTFALLS ], [ CONDUITS ], [ XSECTIONS ], [ COORDINATES ] and [ VERTICES ] partial contents, writing into swmm. Inp file.
3. The rainwater pipe network simulation method based on design parameter equivalent setting according to claim 1, wherein the step of calculating runoff process lines of each sub-catchment area based on the surface production confluence model based on the design parameter equivalent setting in S2 comprises the following specific sub-steps:
1) The design rainfall process line is taken as input, and the flow and rainfall runoff coefficient is equivalent as a target, namely a formula (1) and a formula (2), based onAn exponential method is a formula (3), the related parameters of a Huoton infiltration model is a formula (4) are determined in a trial calculation mode, and net rain is calculatedA process line;
f i ＝f c +(f 0 -f c )e -ki (4)
wherein: i and j are both rainfall periods, i=1, 2, …, n, 15min is taken as the peak flow duration, so j+14 is set; r is R i Net rainfall at i min, mm; p (P) i Rainfall in mm for the i-th minute; c is a rainfall runoff coefficient; c (C) max Is the flow runoff coefficient; f (f) 0 The initial infiltration rate is mm/min; f (f) c To stabilize the infiltration rate, mm/min; k is the attenuation coefficient; f (f) i The infiltration strength is the ith period of time, mm;
(1) flow runoff coefficient equivalence criterion: the average runoff coefficient of the net rainfall of 15min maximum in the net rainfall process line is equal to the flow runoff coefficient;
(2) rainfall runoff coefficient equivalence criterion: the ratio of the total net rainfall to the total rainfall of the whole rainfall is equal to the rainfall runoff coefficient;
2) Setting a linear confluence curve type by taking the water collecting time, the water collecting area and the rain purifying process line of each sub-catchment area as inputs, and calculating the runoff process line of each sub-catchment area by using a constant flow time line method according to a formula (5);
wherein: i. j is the rainfall moment; t is t d Collecting water for the sub-catchment area, and min;Q i the flow at the moment i is L/s; r is R i-j For the net rain strength at i-j, mm, R when i-j is less than or equal to 0 i-j =0; f is the area of regional confluence, m 2 The method comprises the steps of carrying out a first treatment on the surface of the Alpha is a unit conversion coefficient, and the value is 0.167.
4. The rainwater pipe network simulation method based on design parameter equivalent setting as claimed in claim 1, wherein S3 comprises the following sub-steps:
1) Numbering the runoff process lines of the corresponding sub-catchment areas according to the corresponding relation between each sub-catchment area and the inspection well and the ID of the inspection well;
2) According to the requirement of the inp file on the part of [ INFLOWS ], writing all inspection well IDs and corresponding runoff process line IDs into swmm. Inp files, setting 'Constient' and 'Type' as 'FLOW' attributes, and finishing the coupling setting of the model;
3) And writing the runoff process lines of all the sub-catchment areas into the swmm.inp file according to the format requirement of the inp file on the [ TIMESERIES ] part of content, and completing the manufacture of the swmm.inp file.
5. The method for simulating a rainwater pipe network based on design parameter equivalent setting according to claim 2 or 4, wherein the specific requirements of the manufacturing format of the inp file are in accordance with the relevant regulations in "STORMWATERMANAGEMENT MODELUSER' smunualvlersion 5.1".
CN202310549990.3A 2023-05-16 2023-05-16 Rainwater pipe network simulation method based on design parameter equivalent setting Pending CN116702392A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310549990.3A CN116702392A (en) 2023-05-16 2023-05-16 Rainwater pipe network simulation method based on design parameter equivalent setting

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310549990.3A CN116702392A (en) 2023-05-16 2023-05-16 Rainwater pipe network simulation method based on design parameter equivalent setting

Publications (1)

Publication Number Publication Date
CN116702392A true CN116702392A (en) 2023-09-05

Family

ID=87833128

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310549990.3A Pending CN116702392A (en) 2023-05-16 2023-05-16 Rainwater pipe network simulation method based on design parameter equivalent setting

Country Status (1)

Country Link
CN (1) CN116702392A (en)

Similar Documents

Publication Publication Date Title
Aronica et al. Studying the hydrological response of urban catchments using a semi-distributed linear non-linear model
CN105138707A (en) Input file Inp file generation method of SWMM hydraulic model based on GIS model database
CN108984823B (en) Method for determining scale of combined overflow storage tank
CN105117573B (en) Auto hydraulic model building method based on CAD drainage pipeline networks drawing informations
CN110909485B (en) SWMM model parameter self-calibration method based on BP neural network
CN111695305B (en) Water surface line calculation method for rain source type river under condition of no actual measurement hydrological data
CN101692273A (en) Modeling method of on-line hydraulic model of municipal drainage pipe network
CN106897516B (en) Urban underground drainage pipe network reconstruction and reconstruction evaluation method
CN106650053A (en) USRAMS-based optimization design method of rainwater pipe network
Hou et al. Numerical simulation for runoff regulation in rain garden using 2D hydrodynamic model
Peng et al. Simulating rainfall runoff and assessing low impact development (LID) facilities in sponge airport
CN113626959B (en) Sewage pipe network hydraulic model building method based on three-dimensional geographic information
CN104615823B (en) A kind of method circulated using SWMM software simulating riverways water
CN116702392A (en) Rainwater pipe network simulation method based on design parameter equivalent setting
Bu et al. A method for determining reasonable water area ratio based on flood risk and cost-effectiveness in Rainy City
CN117764001A (en) CFD and Modelica-based building siphon drainage pipe network optimization method, device and medium
CN116383950A (en) Multi-platform collaborative outdoor rainwater pipe network parameterization design method and system
Zhang et al. Analysis of regional flooding in the urbanization expansion process based on the SWMM model
CN110298075A (en) Graph model processing method is opened based on ArcGIS+InfoWorks ICM Municipal flowing line one
Venigalla et al. GIS-based engineering management service functions: taking GIS beyond mapping for municipal governments
CN110866215B (en) Method for calculating water collection amount of urban precipitation area
CN113806901A (en) Method for optimizing rainwater pipe system designed based on rainwater management model and storage device
CN112580131A (en) Rapid search connection method for sub-catchment areas based on GIS (geographic information System) and SWMM (Single wire Meter-based Module)
CN109837862B (en) River step development design and scheduling application method based on ecological water level connection
CN112861449A (en) Multi-river-section combined correction model based on flood forecast error inversion

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