CN109960861A - Cienega source dividing method based on runoff yield under saturated storage mode - Google Patents
Cienega source dividing method based on runoff yield under saturated storage mode Download PDFInfo
- Publication number
- CN109960861A CN109960861A CN201910195865.0A CN201910195865A CN109960861A CN 109960861 A CN109960861 A CN 109960861A CN 201910195865 A CN201910195865 A CN 201910195865A CN 109960861 A CN109960861 A CN 109960861A
- Authority
- CN
- China
- Prior art keywords
- basin
- moisture content
- capacity
- evaporation
- unit
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 229920006395 saturated elastomer Polymers 0.000 title claims abstract description 11
- 230000008020 evaporation Effects 0.000 claims abstract description 44
- 238000001704 evaporation Methods 0.000 claims abstract description 44
- 239000002689 soil Substances 0.000 claims abstract description 41
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000010998 test method Methods 0.000 claims description 4
- 230000002045 lasting effect Effects 0.000 abstract description 2
- 230000008595 infiltration Effects 0.000 description 7
- 238000001764 infiltration Methods 0.000 description 7
- 238000000691 measurement method Methods 0.000 description 3
- 238000001556 precipitation Methods 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 238000005303 weighing Methods 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 239000004411 aluminium Substances 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000005527 soil sampling Methods 0.000 description 1
- 238000013022 venting Methods 0.000 description 1
Classifications
-
- 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
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)
- Sewage (AREA)
- Measuring Volume Flow (AREA)
Abstract
The cienega source division methods based on runoff yield under saturated storage mode that the invention discloses a kind of, the following steps are included: determining field capacity first, calculating basin saturation moisture content, basin average moisture content, reservoir capacity exponent, basin reservoir storage maximum value, determine the average initial soil moisture content in basin, the maximum value of basin soil moisture content, basin Pluvial period evaporation capacity and basin evaporation coefficient, secondly judge that the average initial soil moisture content in basin and basin are averaged the relationship of reservoir storage, evaporation capacity is calculated, determines rainfall;Finally judge the size of rainfall and evaporation capacity, calculates rainfall and average the sum of the initial soil moisture content in basin, the relationship between sum of the two and saturation moisture content and field capacity is judged, to calculate rainwash and interflow subsurface drainage.The method of the invention more simply, conveniently, is counted without lasting to every net rainfall, can be expressed with computer programming language, can be preferably integrated among hydrological model.
Description
Technical field
The present invention relates to cienega source dividing method, especially a kind of cienega based on runoff yield under saturated storage mode
Source division methods.
Background technique
Basin Outlet Section flow is made of the runoff component at different water sources, and is regulated and stored work because of its motion path and by basin
Difference, make mutually variant in Outlet Section discharge process feature: surface runoff directly imports the network of waterways, movement velocity by slope surface
Fastly, the concentration time is short, and length is fallen the discharge process of formation suddenly suddenly;Interflow subsurface drainage is flowed out by penetrating into the sluggish flow of phreatic surface,
Movement velocity is slow, and the concentration time is long, and variation is gentle, is the bulk composition of flood water-break end segment flow, often extends to subsequent flood
During water.In actual operation, the different runoff components and its composition ratio for generally requiring analysis, research depth of runoff, often need
It is calculated separately or is simulated by various runoff components, thus need to carry out runoff yield water source division, so that it is determined that basin outlet is disconnected
Face discharge process, regulating and storing, control flood etc. for river provides foundation.Wherein, most basic, the most commonly used is divide direct runoff and underground
Runoff generally divides water source using basin ground water recession and measured discharge graph, the bifurcation of the two often by
It is considered as direct runoff terminating point, connects to rise in discharge curve and rise a little with runoff terminating point, top is direct runoff, and lower part is ground
Lower runoff.In natural storage model, direct runoff and the interflow subsurface drainage in runoff yield are generally divided by stablizing infiltration rate,
Time segments division is carried out to rainfall, the rainfall intensity in the period and the relationship between infiltration rate are analyzed, to inquire into direct runoff and ground
Lower runoff, its calculation formula is:
Wherein, RG is total interflow subsurface drainage;RS is total rainwash;I is i-th of period;PE is precipitation;FC is steady seeps
Rate;R is run-off;fcTo stablize infiltration rate.
In natural storage model, different stabilization infiltration rates is selected, counted runoff component is different.It is general to stablize
Infiltration rate fcIt is calculated using following formula:
fc=Rg/T
Wherein, Rg is that interflow subsurface drainage is deep, and T lasts for net rainfall.
But among practical application, net rainfall lasts T and is difficult directly to count, therefore above formula is difficult to be utilized just to calculate under stabilization
Infiltration rate fc.It is complex by the way of discharge curve when dividing direct runoff and interflow subsurface drainage, it can not be directly by the original
It ought to be used among computer program.It would therefore be highly desirable to solve the above problems.
Summary of the invention
Goal of the invention: the object of the present invention is to provide a kind of cienega source division side based on runoff yield under saturated storage mode
Method, this method are produced and are flowed using saturation moisture content and field capacity as boundary line delimitation earth's surface and underground, calculate rainwash, underground diameter
Flow is easy to be expressed with computer programming language, more simply, conveniently.
Technical solution: in order to achieve the above object, the cienega source of the present invention based on runoff yield under saturated storage mode is drawn
Point method the following steps are included:
Step 1, field capacity is measured in watershed, determines field capacity Wm;
Step 2, basin saturation moisture content W is calculatedb, calculation formula are as follows:
Wb=Wm*1.2
Step 3, basin average moisture content W is determinedM;
Step 4, reservoir capacity exponent B is determined;
Step 5, the maximum value of reservoir capacity curve, i.e. basin reservoir storage maximum value W are calculatedMM, calculation formula are as follows:
WMM=WM*(1+B)
Step 6, field measuring is carried out using soil nmoisture content analyser, obtains soil moisture content data in basin, determine stream
The average initial soil moisture content W in domain;
Step 7, the maximum value A of basin soil moisture content is calculated;
Step 8, in heavy, field is measured water surface evaporation, determines period in basin Pluvial evaporation capacity E;
Step 9, basin evaporation coefficient K is determined;
Step 10, judge the average initial soil moisture content W in basin and basin is averaged reservoir storage WMRelationship, calculate evaporation
Measure EE.
Step 11, it is measured using rainfall gauge or according to meteorological data, determines rainfall P;
Step 12, judge the size of rainfall P Yu evaporation capacity EE, calculate the rainfall P initial soil average with basin and contain
The sum of water W judges sum of the two and saturation moisture content WbWith field capacity WmBetween relationship, calculate rainwash RS and
Interflow subsurface drainage RG.
Further, it is measured using field capacity in field test method or indoor measurement method watershed.
Further, the calculation formula of the maximum value A of basin soil moisture content are as follows:
Further, using the ratio of E601 type evaporator and 20cm diameter evaporation pan evaporation capacity as evaporation coefficient, really
Constant current domain evaporation coefficient K.
Further, in step 10, the calculation method of evaporation capacity EE are as follows:
If W > WM,
EE=K*E
If W≤WM,
Further, the calculation method of rainwash RS and interflow subsurface drainage RG are calculated in step 12 are as follows:
When P-EE≤0
RS=0
RG=0
As P+W > Wb,
RS=P-EE- (Wb-W)
RG=Wb--Wm
Work as Wm≤P+W≤Wb,
RS=P-EE- (Wb-W)
RG=P-EE+W-Wm
As P+W < Wm,
In formula, WbFor basin saturation moisture content, unit mm;
WmFor field capacity, unit mm;
WMFor basin average moisture content, unit mm;
B is storage capacity curve of a river basin index, is constant, no unit;
WMMFor basin reservoir storage maximum value, unit mm;
W is the average initial soil moisture content in basin, unit mm;
A is the maximum value of basin soil moisture content, unit mm;
E is basin Pluvial period evaporation capacity, unit mm;
K is evaporator reduction coefficient, is constant, no unit;
EE is basin evaporation capacity, unit mm;
P is rainfall, unit mm;
RS is rainwash, unit mm;
RG is interflow subsurface drainage, unit mm.
The utility model has the advantages that compared with prior art, the present invention has following remarkable advantage: being divided with original cienega source
Mode is compared, and the water source division methods based on runoff yield under saturated storage more simply, conveniently, are counted without lasting to every net rainfall,
The index for stablizing infiltration rate without calculating, but easily measuring, calculating from saturated soil water content, field capacity etc. is right
Cienega source is divided.And the cienega source division methods in the present invention can be carried out with computer programming language
Expression, can preferably be integrated among hydrological model.
Detailed description of the invention
Fig. 1 is the calculation flow chart of water source division methods of the present invention.
Specific embodiment
Technical solution of the present invention is described further with reference to the accompanying drawing.
As shown in Figure 1, the cienega source division methods based on runoff yield under saturated storage mode the following steps are included:
Step 1, it is measured using field capacity in field test method or indoor measurement method watershed, determines that field is held
Water Wm;
Step 2, basin saturation moisture content W is calculatedb
Wb=Wm*1.2
Step 3, the basin average moisture content is determined by calculation (usually in the rainfall runoff data for searching institute's survey region
Between 120-200mm);Specifically basin average moisture content W can be determined according to the basin feature in research areaM, based on experience value,
The basin average moisture content of moist area of southern China is about 120-150mm, and Semi-humid area is about 150-200mm;
Step 4, based on experience value, determine reservoir capacity exponent B (usually between 0.1-0.4);In Runoff Model
In, B reflects that the degree of irregularity of basin water content, value are related to watershed unit, area: if drainage area is flat less than 100
Square km, then take B=0.1;If drainage area is 100-1000 sq-km, B=0.2-0.3;If drainage area is greater than
1000 sq-kms, then B=0.4;
Step 5, the maximum value of reservoir capacity curve, i.e. basin reservoir storage maximum value W are calculatedMM, calculation formula are as follows:
WMM=WM*(1+B)
Step 6, field measuring is carried out using soil nmoisture content analyser at rainfall initial stage, analyzes soil moisture content in the basin
Data determine the average initial soil moisture content W in basin;
Step 7, the maximum value A of basin soil moisture content, calculation formula are calculated are as follows:
Step 8, in heavy, field is measured water surface evaporation, determines period in basin Pluvial evaporation capacity E;
Step 9, stream is determined as evaporation coefficient using the ratio of E601 type evaporator and 20cm diameter evaporation pan evaporation capacity
Domain evaporation coefficient K;
Step 10, judge the average initial soil moisture content W in basin and basin is averaged reservoir storage WMRelationship, calculate evaporation
Measure EE;
Step 11, using rainfall gauge measurement or meteorological data, rainfall P is determined;
Step 12, judge the size of rainfall P Yu evaporation capacity EE, calculate the rainfall P initial soil average with basin and contain
The sum of water W judges sum of the two and saturation moisture content WbWith field capacity WmBetween relationship, calculate rainwash RS and
Interflow subsurface drainage RG.
In step 1, field test method specifically: choose representational location, the cell of certain area is impaled, through overcharging
Divide and pour water, after the extra bulk water of venting, measures the maximum kept in soil layer and poise water.Indoor measurement method specifically: ring
Skill in using a kitchen knife in cookery field soil sampling, soil sample air-dry, the weighing of aluminium box, original state tool changing soil sample impregnate, seep under moisture, weigh wet soil soil sample, constant temperature
It dries, be cooled to room temperature, weighing dry soil samples, calculating achievement.
The step 3 includes the following steps:
Step 301, according to special drought is in early period in arid grade selection basin, (continuously without rainfall number of days, spring was at 61 days
Above, summer in 46 days or more, autumn and winter at 91 days or more) or weight drought (continuously without rainfall number of days, spring 46-60 days or more,
Summer is in 36-45 days or more, autumn and winter at 71-90 days or more) in the period of, generation, rainfall grade are sudden and violent after choosing drouth
More than rain (24 hourly rainfall depth 50-100mm) and heavy rain rainfall is king-sized, provides the history rainfall runoff of full flow anomaly
Material, physical resource can search acquisition, arid, rainfall in China Meteorological data network, Water Year Book, meteorological department, Hydrology department
Quantitative description determined according to the arid grade of national standard, rainfall grade;
Step 302, it is calculated according to the water balance equation of a precipitation: reservoir storage W in basin at the end of rainfall2Just
Beginning reservoir storage W1=rainfall P- evaporation capacity E- run-off R;
Step 303, basin average moisture content W is determinedM: since rainfall is very arid early period, W1≈ 0, it is believed that stream after rainfall
Domain has stored completely, therefore W2≈WM。
Further, the step 10 includes the following steps:
Step 1001, the average initial soil moisture content W > basin average moisture content W in basinM, evaporation capacity EE=evaporation folding
Calculate COEFFICIENT K * period in Pluvial evaporation capacity E;
Step 1002, the average initial soil moisture content W in basin≤basin average moisture content WM, evaporation capacity EE=evaporation folding
Calculate the average basin the initial soil moisture content W/ average moisture content W in the basin COEFFICIENT K * period in Pluvial evaporation capacity E*M;
Further, the step 12 includes the following steps:
Step 1201: rainfall P is less than evaporation capacity EE, rainwash RS=0, interflow subsurface drainage RG=0;
Step 1202: the average initial soil moisture content W > < basin saturation moisture content W in the basin rainfall P+b, earth's surface diameter
Flow RS=P-EE- (Wb- W), interflow subsurface drainage RG=Wb-Wm;
Step 1203: basin field capacity WmThe basin≤rainfall P+ initial time soil moisture content W < is saturated aqueous
Measure Wb, rainwashInterflow subsurface drainage RG=P-
EE-R;
Step 1204: rainfall P+ initial time soil moisture content W < basin field capacity Wm, rainwash RS=0,
Interflow subsurface drainage
Claims (6)
1. a kind of cienega source division methods based on runoff yield under saturated storage mode, it is characterised in that: the following steps are included:
Step 1, field capacity is measured in watershed, determines field capacity Wm;
Step 2, basin saturation moisture content W is calculatedb, calculation formula are as follows:
Wb=Wm*1.2
Step 3, basin average moisture content W is determinedM;
Step 4, reservoir capacity exponent B is determined;
Step 5, the maximum value of reservoir capacity curve, i.e. basin reservoir storage maximum value W are calculatedMM, calculation formula are as follows:
WMM=WM*(1+B)
Step 6, field measuring is carried out using soil nmoisture content analyser, obtains soil moisture content data in basin, determine that basin is flat
Equal initial soil moisture content W;
Step 7, the maximum value A of basin soil moisture content is calculated;
Step 8, in heavy, field is measured water surface evaporation, determines period in basin Pluvial evaporation capacity E;
Step 9, basin evaporation coefficient K is determined;
Step 10, judge the average initial soil moisture content W in basin and basin is averaged reservoir storage WMRelationship, calculate evaporation capacity EE.
Step 11, it is measured using rainfall gauge or according to meteorological data, determines rainfall P;
Step 12, judge the size of rainfall P Yu evaporation capacity EE, calculate rainfall family and the average initial soil moisture content in basin
The sum of W judges sum of the two and saturation moisture content WbWith field capacity WmBetween relationship, calculate rainwash RS and underground
Runoff RG.
2. cienega source according to claim 1 division methods, it is characterised in that: use field test method or interior
Field capacity is measured in measuring method watershed.
3. cienega source according to claim 1 division methods, it is characterised in that: the maximum of basin soil moisture content
The calculation formula of value A are as follows:
4. cienega source according to claim 1 division methods, it is characterised in that: using E601 type evaporator with
The ratio of 20cm diameter evaporation pan evaporation capacity determines basin evaporation coefficient K as evaporation coefficient.
5. cienega source according to claim 1 division methods, it is characterised in that: in step 10, evaporation capacity EE's
Calculation method are as follows:
If W > WM,
EE=K*E
If W≤WM,
6. cienega source according to claim 1 division methods, it is characterised in that: calculate rainwash in step 12
The calculation method of RS and interflow subsurface drainage RG are as follows:
When P-EE≤0
RS=0
RG=0
As P+W > Wb,
RS=P-EE- (Wb-W)
RG=Wb-Wm
Work as Wm≤P+W≤Wb,
RS=P-EE- (Wb-W)
RG=P-EE+W-Wm
As P+W < Wm,
In formula, WbFor basin saturation moisture content, unit mm;
WmFor field capacity, unit mm;
WMFor basin average moisture content, unit mm;
B is storage capacity curve of a river basin index, is constant, no unit;
WMMFor basin reservoir storage maximum value, unit mm;
W is the average initial soil moisture content in basin, unit mm;
A is the maximum value of basin soil moisture content, unit mm;
E is basin Pluvial period evaporation capacity, unit mm;
K is evaporator reduction coefficient, is constant, no unit;
EE is basin evaporation capacity, unit mm;
P is rainfall, unit mm;
RS is rainwash, unit mm;
RG is interflow subsurface drainage, unit mm.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910195865.0A CN109960861B (en) | 2019-03-14 | 2019-03-14 | Earth surface and underground water source segmentation method based on full-accumulation flow-producing mode |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910195865.0A CN109960861B (en) | 2019-03-14 | 2019-03-14 | Earth surface and underground water source segmentation method based on full-accumulation flow-producing mode |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109960861A true CN109960861A (en) | 2019-07-02 |
CN109960861B CN109960861B (en) | 2023-05-09 |
Family
ID=67024413
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910195865.0A Active CN109960861B (en) | 2019-03-14 | 2019-03-14 | Earth surface and underground water source segmentation method based on full-accumulation flow-producing mode |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109960861B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113656745A (en) * | 2021-08-26 | 2021-11-16 | 中国水利水电科学研究院 | Method for calculating runoff production reference groundwater burial depth reflecting rainfall runoff relation |
CN113657047A (en) * | 2021-08-13 | 2021-11-16 | 河海大学 | Production flow calculation method and device based on double free reservoirs and storage medium |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102314554A (en) * | 2011-08-08 | 2012-01-11 | 大唐软件技术股份有限公司 | Land-atmosphere coupling-based method and system for flood forecast of minor watersheds |
-
2019
- 2019-03-14 CN CN201910195865.0A patent/CN109960861B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102314554A (en) * | 2011-08-08 | 2012-01-11 | 大唐软件技术股份有限公司 | Land-atmosphere coupling-based method and system for flood forecast of minor watersheds |
Non-Patent Citations (2)
Title |
---|
宋丹丹 等: "江苏省南水北调受水区水资源配置", 《南水北调与水利科技》 * |
杨树滩 等: "苏北四市应用单层蓄满产流模型分割旱地地下水研究", 《治淮》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113657047A (en) * | 2021-08-13 | 2021-11-16 | 河海大学 | Production flow calculation method and device based on double free reservoirs and storage medium |
CN113656745A (en) * | 2021-08-26 | 2021-11-16 | 中国水利水电科学研究院 | Method for calculating runoff production reference groundwater burial depth reflecting rainfall runoff relation |
CN113656745B (en) * | 2021-08-26 | 2024-03-12 | 中国水利水电科学研究院 | Calculation method for runoff yield reference underground water burial depth reflecting rainfall runoff relation |
Also Published As
Publication number | Publication date |
---|---|
CN109960861B (en) | 2023-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Brocca et al. | On the estimation of antecedent wetness conditions in rainfall–runoff modelling | |
Yu et al. | Soil moisture variations at different topographic domains and land use types in the semi-arid Loess Plateau, China | |
Winchell et al. | On the simulation of infiltration‐and saturation‐excess runoff using radar‐based rainfall estimates: Effects of algorithm uncertainty and pixel aggregation | |
Feddes | Simulation of field water use and crop yield | |
Kohfahl et al. | Determining groundwater recharge and vapor flow in dune sediments using a weighable precision meteo lysimeter | |
CN109960861A (en) | Cienega source dividing method based on runoff yield under saturated storage mode | |
Ramírez et al. | Tropical Montane Cloud Forests in the Orinoco river basin: The role of soil organic layers in water storage and release | |
Ding et al. | Effects of forest cover type and ratio changes on runoff and its components | |
Plattner et al. | Spatial variability of snow accumulation on Vernagtferner, Austrian Alps, in winter 2003/2004 | |
Wellings et al. | Physical controls of water movement in the unsaturated zone | |
Silva et al. | Nonlinear models for soil moisture sensor calibration in tropical mountainous soils | |
Venkatesh et al. | Modelling soil moisture under different land covers in a sub-humid environment of Western Ghats, India | |
Alivio et al. | Investigating the reduction of rainfall intensity beneath an urban deciduous tree canopy | |
Scholl | Soil moisture flux and evapotranspiration determined from soil hydraulic properties in a chaparral stand | |
Faridah et al. | Infiltration model of Mediterranean soil with clay texture. | |
Durães et al. | Groundwater recharge behavior based on surface runoff hydrographs in two basins of the Minas Gerais State | |
Hume | Determination of infiltration characteristics by volume balance for border check irrigation | |
Marc et al. | Application of the deuterium tracing method for the estimation of tree sap flow and stand transpiration of a beech forest (Fagus silvatica L.) in a mountainous Mediterranean region | |
Peper et al. | Comparison of four foliar and woody biomass estimation methods applied to open-grown deciduous trees | |
McCown et al. | Assessment of available water storage capacity of soils with restricted subsoil permeability | |
Jin et al. | Evapotranspiration of xerophytic shrub Salsola passerina and Reaumuria soongorica in an arid desert ecosystem of NW China | |
Wilcox et al. | Soil Moisture Studies: II. Some Relationships Between Moisture Measurements and Mechanical Analysis | |
Wopereis et al. | Comparing the performance of a soil-water balance model using measured and calibrated hydraulic conductivity data: a case study for dryland rice | |
CN111077273A (en) | Method for measuring vegetation growth index and scale effect of hydrological element slope | |
Deng et al. | Litter fall production in west-african forests and plantations |
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 |