CN107766985A - A kind of river mouth salinity Forecasting Methodology - Google Patents

A kind of river mouth salinity Forecasting Methodology Download PDF

Info

Publication number
CN107766985A
CN107766985A CN201711038900.5A CN201711038900A CN107766985A CN 107766985 A CN107766985 A CN 107766985A CN 201711038900 A CN201711038900 A CN 201711038900A CN 107766985 A CN107766985 A CN 107766985A
Authority
CN
China
Prior art keywords
river mouth
mrow
salinity
day
homohalin
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.)
Granted
Application number
CN201711038900.5A
Other languages
Chinese (zh)
Other versions
CN107766985B (en
Inventor
王青
孔俊
汪玉平
张梦茹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN201711038900.5A priority Critical patent/CN107766985B/en
Publication of CN107766985A publication Critical patent/CN107766985A/en
Application granted granted Critical
Publication of CN107766985B publication Critical patent/CN107766985B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Marketing (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a kind of river mouth salinity Forecasting Methodology, including:River mouth footpath flow data and river mouth Salinity Data are gathered, establishes database;Average daily data processing is carried out to river mouth footpath flow data and river mouth Salinity Data, obtains the average daily run-off in river mouth and river mouth day homohalin value;The average daily run-off statistical model of river mouth day homohalin and river mouth is estimated, establishes the fit correlation formula of the average daily run-off of river mouth day homohalin and river mouth;River mouth day homohalin and memory salinity statistical model are estimated, establishes the fit correlation formula of river mouth day homohalin and memory salinity;Consider the average daily run-off in river mouth and remember the influence of salinity, establish river mouth salinity prediction statistical model, statistical model prediction river mouth salinity is predicted according to river mouth salinity.The salt angle value that the present invention dynamically predicts next day, prediction result are accurate;Forecast model only has a unknown parameter, simplifies salinity forecasting problem.The present invention can fetch water for Hekou Area and provide foundation and technical support.

Description

A kind of river mouth salinity Forecasting Methodology
Technical field
The present invention relates to a kind of river mouth salinity Forecasting Methodology, belongs to river mouth salty tide technical field.
Background technology
Hekou Area is the area that crosses in ocean and river, is the area that economy is most active in the world, population is most dense One.And as social economy is fast-developing, more and more deep to the resources development and utilization of Hekou Area, river mouth salty tide traces back problem Also it is following, directly affect the fresh water supply of the resident living of Hekou Area, field irrigation, urban industry production etc..Therefore, Accurate Prediction river mouth salty tide, which is traced back, to have important practical significance.
The content of the invention
It is an object of the invention to overcome deficiency of the prior art, there is provided a kind of river mouth salinity Forecasting Methodology, Neng Gouzhun Really prediction Hekou Area salinity, foundation is provided for Hekou Area water intaking.
To reach above-mentioned purpose, the technical solution adopted in the present invention is:A kind of river mouth salinity Forecasting Methodology, including it is as follows Step:River mouth footpath flow data and river mouth Salinity Data are gathered, establishes database;
Average daily data processing is carried out to river mouth footpath flow data and river mouth Salinity Data, obtains the average daily run-off in river mouth and river mouth Day homohalin value;
The average daily run-off statistical model of river mouth day homohalin and river mouth is estimated, establishes the average daily footpath of river mouth day homohalin and river mouth The fit correlation formula of flow;
River mouth day homohalin and memory salinity statistical model are estimated, the fitting for establishing river mouth day homohalin and memory salinity is closed It is formula;
Consider the average daily run-off in river mouth and remember the influence of salinity, establish river mouth salinity prediction statistical model, according to River mouth salinity prediction statistical model prediction river mouth salinity.
The survey station arranged by Hekou Area to be predicted gathers river mouth footpath flow data and river mouth Salinity Data.
Acquisition time was at intervals of 1 hour.
River mouth day homohalin and the average daily run-off in river mouth fit correlation formula it is as follows:
SQ=aQexp(-bQQ) (1)
In formula:SQOnly to consider the day homohalin predicted value under the average daily run-off Q in river mouth;Q is the average daily runoff in river mouth;aQWith bQFor fitting coefficient.
The specific method for establishing the fit correlation formula of river mouth day homohalin and memory salinity is as follows:
Current time day homohalin measured value is calculated using autoregressive method:
St=F (St-1,St-2,St-3) (2)
In formula:StFor the day homohalin measured value of the t days;St-iFor the day homohalin measured value of the t-i days;
Using dynamic model, and combine gamma distribution function, day homohalin predictor formula it is as follows:
In formula:ωt-iFor important coefficient;α is form parameter;β is scale parameter;To consider the of important coefficient The day homohalin measured value of t-i days;
In formula:To consider the day homohalin predicted value of the t+1 days of important coefficient;It is for fitting Number.
The river mouth salinity prediction statistical model is as follows:
In formula:St+1For the day homohalin predicted value of the t+1 days;Qt+1For the average daily runoff of the t+1 days;A, b, c are fitting Coefficient;E is constant term.
Compared with prior art, the beneficial effect that is reached of the present invention is:The salt angle value of next day is dynamically predicted, and And prediction result is very accurate;Forecast model only has an average daily runoff of main unknown parameter, and salinity is predicted challenge letter Dan Hua, there is simple and practical outstanding feature;Foundation and technical support are provided for Hekou Area water intaking.
Brief description of the drawings
Fig. 1 is the flow chart of the present invention;
Fig. 2 is the geographical position figure of embodiment survey region;
Fig. 3 is the average daily run-off in river mouth and the Day-to-day variability figure of day homohalin;
Fig. 4 is that the different fitting forms of the average daily run-off in river mouth and day homohalin compare figure;
Wherein, a represents exponential function;B represents power function;
Fig. 5 is only to consider in September, 2011~2013 year flat hilllock day homohalin prediction in December under the influence of the average daily run-off in river mouth Value and day homohalin measured value comparison diagram;
Fig. 6 is to consider the Nian8Yue Ping hilllocks day of the average daily run-off in river mouth and memory Effects of Salinity lower in September, 2011~2012 Homohalin predicted value and day homohalin measured value comparison diagram;
Fig. 7 is to consider the average daily run-off in river mouth and the lower flat hilllock in September, 2012~2013 year December of memory Effects of Salinity Day homohalin predicted value and day homohalin measured value comparison diagram.
Embodiment
The invention will be further described below in conjunction with the accompanying drawings.Following examples are only used for clearly illustrating the present invention Technical scheme, and can not be limited the scope of the invention with this.
As shown in figure 1, being the flow chart of salinity Forecasting Methodology in river mouth of the present invention, comprise the following steps:
Step 1:Estimation range is treated using the survey station arranged in Hekou Area to be monitored, and gathers river mouth footpath flow data With river mouth Salinity Data, database is established.Acquisition time was at intervals of 1 hour.
Step 2:Average daily data processing is carried out to river mouth footpath flow data and river mouth Salinity Data, obtains the average daily runoff in river mouth Amount and river mouth day homohalin value;
Average daily data processing refers to:According to the river mouth runoff measured data and river mouth salinity measured data gathered per hour, Arithmetic average is made to both respectively in units of day, accordingly, you can obtain the average daily run-off in river mouth and river mouth day homohalin value.
Step 3:The average daily run-off statistical model of river mouth day homohalin and river mouth is estimated, establishes river mouth day homohalin and river The fit correlation formula of the average daily run-off of mouth;
Exemplified by Pearl River Estuary, the regional geography position is as shown in Figure 2.Wherein upstream runoff survey station is Ma Kou stations and three water To stand, salinity survey station is flat hilllock station, wherein, calculated to simplify, Ma Kou is stood and is added with the footpath flow valuve of three water stations, with addition Obtained runoff and studied.Upstream runoff and Salinity Data are all the actual measurement per hour in September, 2011~2013 year December Data.
From the figure 3, it may be seen that there is obvious dependency relation between average daily runoff and day homohalin, with subtracting for average daily runoff Small, the day homohalin of river mouth section significantly increases;With the increase of average daily runoff, the day homohalin of river mouth section is substantially reduced, even Close to 0.By the observation being distributed to day homohalin and the scatterplot of average daily runoff, comparison index functional form and power function shape Formula, as shown in Figure 4, it is known that goodness of fit index R corresponding to exponential function2For 0.29, goodness of fit index R corresponding to power function2 For 0.27, exponential function form is more excellent, is closed so as to establish the exponential fitting of the average daily run-off of river mouth day homohalin and river mouth System:
SQ=aQexp(-bQQ) (1)
In formula:SQOnly to consider the day homohalin predicted value under the average daily run-off Q in river mouth, unit psu;Q is river mouth day Equal runoff, unit m3/s;aQAnd bQFor fitting coefficient.Prediction result is as shown in Figure 5.
Step 4:River mouth day homohalin and memory salinity statistical model are estimated, establishes river mouth day homohalin and memory salinity Fit correlation formula;
Consider the influence of the memory effect of history salinity, current time day homohalin measured value is calculated using autoregressive method:
St=F (St-1,St-2,St-3) (2)
In formula:StFor the day homohalin measured value of the t days, unit pus;St-iFor the day homohalin actual measurement of the t-i days Value, unit pus;
Using dynamic model, and combine gamma distribution function, day homohalin predictor formula it is as follows:
In formula:ωt-iFor important coefficient;α is form parameter;β is scale parameter;To consider the of important coefficient The day homohalin measured value of t-i days, unit psu;
In formula:To consider the day homohalin predicted value of the t+1 days of important coefficient, unit pus; For fitting coefficient.
According to formula (5), by the day homohalin measured value of t-1, t-2 and t-3 days for introducing consideration important coefficientWithIt is fitted the day homohalin measured value S of the t daystObtain fitting coefficientAgain by fitting coefficient band Enter formula (6), by the day homohalin measured value of t, t-1 and t-2 days of consideration important coefficientWithPredict t The day homohalin value of+1 day, so far complete a circulation.
Step 5:Consider the average daily run-off in river mouth and remember the influence of salinity, establish river mouth salinity prediction statistics mould Type, statistical model prediction river mouth salinity is predicted according to river mouth salinity.
The river mouth salinity prediction statistical model is as follows:
In formula:St+1For the day homohalin predicted value of the t+1 days, unit psu;Qt+1It is single for the average daily runoff of the t+1 days Position is m3/s;A, b, c are fitting coefficient;E is constant term.Prediction result is as shown in Figure 6 and Figure 7.
The present invention is established based on average daily by providing a kind of river mouth salinity Forecasting Methodology based on runoff and memory salinity The river mouth salt solution forecast model of runoff and memory salinity, the salt angle value for dynamically predicting next day, and prediction result is very Accurately;Forecast model only has an average daily runoff of main unknown parameter, and salinity prediction challenge is simplified, and has easy to be real Outstanding feature;Foundation and technical support are provided for Hekou Area water intaking.
Described above is only the preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art For member, without departing from the technical principles of the invention, some improvement and deformation can also be made, these are improved and deformation Also it should be regarded as protection scope of the present invention.

Claims (6)

1. a kind of river mouth salinity Forecasting Methodology, it is characterised in that comprise the following steps:
River mouth footpath flow data and river mouth Salinity Data are gathered, establishes database;
Average daily data processing is carried out to river mouth footpath flow data and river mouth Salinity Data, the average daily run-off in river mouth is obtained and river mouth is average daily Salt angle value;
The average daily run-off statistical model of river mouth day homohalin and river mouth is estimated, establishes the average daily run-off of river mouth day homohalin and river mouth Fit correlation formula;
River mouth day homohalin and memory salinity statistical model are estimated, establishes the fit correlation of river mouth day homohalin and memory salinity Formula;
Consider the average daily run-off in river mouth and remember the influence of salinity, river mouth salinity prediction statistical model is established, according to river mouth Salinity prediction statistical model prediction river mouth salinity.
2. salinity Forecasting Methodology in river mouth according to claim 1, it is characterised in that arranged by Hekou Area to be predicted Survey station gathers river mouth footpath flow data and river mouth Salinity Data.
3. salinity Forecasting Methodology in river mouth according to claim 2, it is characterised in that acquisition time was at intervals of 1 hour.
4. salinity Forecasting Methodology in river mouth according to claim 1, it is characterised in that river mouth day the homohalin and average daily footpath in river mouth The fit correlation formula of flow is as follows:
SQ=aQexp(-bQQ) (1)
In formula:SQOnly to consider the day homohalin predicted value under the average daily run-off Q in river mouth;Q is the average daily runoff in river mouth;aQAnd bQFor Fitting coefficient.
5. salinity Forecasting Methodology in river mouth according to claim 4, it is characterised in that establish river mouth day homohalin and memory salt The specific method of the fit correlation formula of degree is as follows:
Current time day homohalin measured value is calculated using autoregressive method:
St=F (St-1,St-2,St-3) (2)
In formula:StFor the day homohalin measured value of the t days;St-iFor the day homohalin measured value of the t-i days;
Using dynamic model, and combine gamma distribution function, day homohalin predictor formula it is as follows:
<mrow> <msub> <mi>&amp;omega;</mi> <mrow> <mi>t</mi> <mo>-</mo> <mi>i</mi> </mrow> </msub> <mo>=</mo> <mi>G</mi> <mi>a</mi> <mi>m</mi> <mi>m</mi> <mi>a</mi> <mrow> <mo>(</mo> <mi>&amp;alpha;</mi> <mo>,</mo> <mi>&amp;beta;</mi> <mo>,</mo> <mi>i</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>&amp;beta;</mi> <mi>&amp;alpha;</mi> </msup> <mfrac> <mn>1</mn> <mrow> <mi>&amp;Gamma;</mi> <mrow> <mo>(</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <msup> <mrow> <mo>(</mo> <mi>i</mi> <mo>)</mo> </mrow> <mrow> <mi>&amp;alpha;</mi> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&amp;beta;</mi> <mi>i</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msubsup> <mi>S</mi> <mrow> <mi>t</mi> <mo>-</mo> <mi>i</mi> </mrow> <mo>*</mo> </msubsup> <mo>=</mo> <msub> <mi>&amp;omega;</mi> <mrow> <mi>t</mi> <mo>-</mo> <mi>i</mi> </mrow> </msub> <msub> <mi>S</mi> <mrow> <mi>t</mi> <mo>-</mo> <mi>i</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
In formula:ωt-iFor important coefficient;α is form parameter;β is scale parameter;To consider the t-i days of important coefficient Day homohalin measured value;
In formula:To consider the day homohalin predicted value of the t+1 days of important coefficient;For fitting coefficient.
6. salinity Forecasting Methodology in river mouth according to claim 5, it is characterised in that the river mouth salinity predicts statistical model It is as follows:
<mrow> <msub> <mi>S</mi> <mrow> <mi>t</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mi>a</mi> <mi> </mi> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mi>bQ</mi> <mrow> <mi>t</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>cS</mi> <mrow> <mi>t</mi> <mo>+</mo> <mn>1</mn> </mrow> <mo>*</mo> </msubsup> <mo>+</mo> <mi>e</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
In formula:St+1For the day homohalin predicted value of the t+1 days;Qt+1For the average daily runoff of the t+1 days;A, b, c are fitting coefficient; E is constant term.
CN201711038900.5A 2017-10-31 2017-10-31 Estuary salinity prediction method Active CN107766985B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711038900.5A CN107766985B (en) 2017-10-31 2017-10-31 Estuary salinity prediction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711038900.5A CN107766985B (en) 2017-10-31 2017-10-31 Estuary salinity prediction method

Publications (2)

Publication Number Publication Date
CN107766985A true CN107766985A (en) 2018-03-06
CN107766985B CN107766985B (en) 2021-09-17

Family

ID=61271802

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711038900.5A Active CN107766985B (en) 2017-10-31 2017-10-31 Estuary salinity prediction method

Country Status (1)

Country Link
CN (1) CN107766985B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109840313A (en) * 2019-01-22 2019-06-04 河海大学 A kind of Hekou Area less salt water intaking determination method
CN110110942A (en) * 2019-05-20 2019-08-09 水利部交通运输部国家能源局南京水利科学研究院 A kind of salinity effect method of riverine getting water from water head site mouth
CN110287624A (en) * 2019-06-28 2019-09-27 宁波市气象台 A method of for cultivating salinity effect model foundation in bay during typhoon influence
CN115081740A (en) * 2022-07-20 2022-09-20 中山大学 Estuary area salinity prediction method, system and equipment based on cellular automaton
CN116776539A (en) * 2023-05-09 2023-09-19 珠江水利委员会珠江水利科学研究院 Salt tide forecasting method and system based on cross wavelet analysis

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288229A (en) * 2011-05-11 2011-12-21 中国水利水电科学研究院 Runoff quantity simulating and predicting method
CN103617306A (en) * 2013-10-23 2014-03-05 广东省水文局佛山水文分局 River mouth salt tide forecasting method based on tides and runoff

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288229A (en) * 2011-05-11 2011-12-21 中国水利水电科学研究院 Runoff quantity simulating and predicting method
CN103617306A (en) * 2013-10-23 2014-03-05 广东省水文局佛山水文分局 River mouth salt tide forecasting method based on tides and runoff

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
孙昭华等: "长江口北支倒灌影响区盐度预测经验模型", 《水科学进展》 *
官明开等: "基于一维盐度模型的珠江河网咸潮上溯距离分析", 《水运工程》 *
诸裕良等: "珠江口盐水入侵预测模式研究", 《水利学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109840313A (en) * 2019-01-22 2019-06-04 河海大学 A kind of Hekou Area less salt water intaking determination method
CN109840313B (en) * 2019-01-22 2023-08-01 河海大学 Low-salt water intake judging method for estuary area
CN110110942A (en) * 2019-05-20 2019-08-09 水利部交通运输部国家能源局南京水利科学研究院 A kind of salinity effect method of riverine getting water from water head site mouth
CN110287624A (en) * 2019-06-28 2019-09-27 宁波市气象台 A method of for cultivating salinity effect model foundation in bay during typhoon influence
CN115081740A (en) * 2022-07-20 2022-09-20 中山大学 Estuary area salinity prediction method, system and equipment based on cellular automaton
CN116776539A (en) * 2023-05-09 2023-09-19 珠江水利委员会珠江水利科学研究院 Salt tide forecasting method and system based on cross wavelet analysis
CN116776539B (en) * 2023-05-09 2023-12-12 珠江水利委员会珠江水利科学研究院 Salt tide forecasting method and system based on cross wavelet analysis

Also Published As

Publication number Publication date
CN107766985B (en) 2021-09-17

Similar Documents

Publication Publication Date Title
CN107766985A (en) A kind of river mouth salinity Forecasting Methodology
CN113379109B (en) Runoff forecasting method based on prediction model self-adaption
CN110991687B (en) Water resource scheduling optimization method based on empirical model
Xing et al. Projection of future runoff change using climate elasticity method derived from Budyko framework in major basins across China
Du et al. Optimal sensor placement and measurement of wind for water quality studies in urban reservoirs
Wang et al. Simulating the impact of climate change on runoff in a typical river catchment of the Loess Plateau, China
Wu et al. Runoff forecasting using convolutional neural networks and optimized bi-directional long short-term memory
CN108596998A (en) A kind of rainfall runoff correlation drawing drawing method based on Copula functions
CN105323321A (en) Water networking system
Huang et al. Dominant climatic factors driving annual runoff changes at the catchment scale across China
CN105320827A (en) Rainfall runoff forecasting system
CN108830419B (en) Cascade reservoir group-entering flow joint prediction method based on ECC post-processing
CN102567635A (en) Method for quantificationally distinguishing contributions of different factors in water cycle evolution process
CN113204865B (en) Urban design rainstorm type calculation method and system based on intrinsic orthogonal decomposition
Guo et al. VIC distributed hydrological model to predict climate change impact in the Hanjiang basin
Zhang et al. Calculation of ecological water requirements of urban rivers using a hydrological model: A case study of Beiyun River
CN112686426A (en) Incoming water quantity early warning method and system based on hydropower station basin key points
CN111984700B (en) Moon runoff prediction method based on hydrothermal coupling balance principle
Copetti et al. Impacts of climate change on water quality
Ahmad et al. An efficient Bayesian modelling of extreme winds in the favour​ of energy generation in Pakistan
Wang et al. Assessment of future socioeconomic drought based on CMIP6: evolution, driving factors and propagation
Kreuwel et al. Forecasting day-ahead 1-minute irradiance variability from numerical weather predictions
Liu et al. Effect of the quality of streamflow forecasts on the operation of cascade hydropower stations using stochastic optimization models
CN106355264A (en) Combined prediction method of reference crop evapotranspiration
Zhang et al. Effects of climate change and land use on runoff in the Huangfuchuan Basin, China

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant