CN112926275A - Method for constructing compact sandstone reservoir water yield prediction model - Google Patents

Method for constructing compact sandstone reservoir water yield prediction model Download PDF

Info

Publication number
CN112926275A
CN112926275A CN202110069991.9A CN202110069991A CN112926275A CN 112926275 A CN112926275 A CN 112926275A CN 202110069991 A CN202110069991 A CN 202110069991A CN 112926275 A CN112926275 A CN 112926275A
Authority
CN
China
Prior art keywords
water
gas
relative permeability
prediction model
constructing
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
CN202110069991.9A
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 National Petroleum Corp
China Petroleum Logging Co Ltd
Original Assignee
China National Petroleum Corp
China Petroleum Logging 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 National Petroleum Corp, China Petroleum Logging Co Ltd filed Critical China National Petroleum Corp
Priority to CN202110069991.9A priority Critical patent/CN112926275A/en
Publication of CN112926275A publication Critical patent/CN112926275A/en
Pending legal-status Critical Current

Links

Images

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
    • 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
    • G06Q—INFORMATION 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/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION 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
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Fluid Mechanics (AREA)
  • Evolutionary Computation (AREA)
  • Operations Research (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Game Theory and Decision Science (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Development Economics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Quality & Reliability (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Agronomy & Crop Science (AREA)
  • Animal Husbandry (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Mining & Mineral Resources (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

The invention discloses a method for constructing a compact sandstone reservoir water production rate prediction model, which comprises the steps of calculating a gas-water relative permeability ratio through a core analysis gas-water relative permeability experiment, constructing a new parameter containing a water saturation variable, establishing a fitting model of the gas-water relative permeability ratio and the new parameter containing the water saturation variable, and substituting a simplified flow splitting equation to construct a prediction model of new water production rate and water saturation. Through analysis of actual data, the new model has higher accuracy of water yield predicted by the physical method of the traditional oil layer, is suitable for different oil and gas reservoirs (different oil and gas reservoirs have different oil/gas-water viscosity ratios), has good universality and has good popularization and application values.

Description

Method for constructing compact sandstone reservoir water yield prediction model
Technical Field
The invention belongs to the technical field of reservoir evaluation, and particularly relates to a method for constructing a compact sandstone reservoir water yield prediction model.
Background
The water yield evaluation is a basic content in the oil field reservoir evaluation work, and the water yield is predicted by mainly utilizing the relative permeability analysis data of two-phase fluid and adopting the traditional oil layer physics, an incomplete Beta function and other fitting calculation flow splitting equations in the prior art. Due to the non-Darcy seepage characteristic of the tight sandstone reservoir, the water yield prediction precision of the traditional oil layer physics fitting model is low under the conditions of medium-low water content and high water content, and the fine evaluation of the water production condition of the reservoir is not facilitated; the incomplete Beta function method requires a solid mathematical work base, understands computer programming, has relatively large application difficulty, and is not beneficial to popularization and application of the method.
Disclosure of Invention
The invention aims to overcome the defects and provides a method for constructing a compact sandstone reservoir water yield prediction model, which is simple in model method, high in predicted water yield precision and suitable for different oil and gas reservoirs.
In order to achieve the above object, the present invention comprises the steps of:
performing a gas-water relative permeability experiment on a core of a target reservoir in a research area to obtain data of water saturation, gas-phase relative permeability and water-phase relative permeability;
calculating a gas-water relative permeability ratio value by utilizing a core gas-water relative permeability experiment, and constructing a parameter containing a water saturation variable;
analyzing the data obtained in the step one and the parameters obtained in the step two, and establishing a fitting model of the gas-water relative permeability ratio and new parameters containing water saturation variables;
and step four, substituting the fitting model obtained in the step three into a simplified flow rate equation to complete the construction of a new prediction model of the water production rate and the water saturation.
In the first step, the experimental data of the gas-water relative permeability of the rock core is obtained according to the flow specified by the Standard of the method for measuring the relative permeability of the two-phase fluid in GB/T28912-2012 rock.
In the second step, the calculation formula of the gas-water relative permeability ratio is as follows:
y＝krg/krw
in the formula, krgIs gas relative permeability, krwRelative water permeability.
In the second step, a new parameter calculation formula containing the water saturation variable is as follows:
Figure BDA0002905430740000021
in the formula, SwThe water saturation.
In the third step, the expression of the fitting model is as follows:
Figure BDA0002905430740000022
wherein y is the gas-water relative permeability ratio krg/krwX is a new parameter containing a water saturation variable
Figure BDA0002905430740000023
SwFor water saturation, a and b are model fitting parameters, a>0，b<0。
In the fourth step, the simplified equation expression of the flow rate is as follows:
Figure BDA0002905430740000024
in the formula: f. ofwTo yield water, muwIs the viscosity of water, mugIs the viscosity of gas, krgIs gas relative permeability, krwRelative water permeability.
In the fourth step, the expression of the new prediction model of the water yield and the water saturation is as follows:
Figure BDA0002905430740000025
in the formula: mu.swIs the viscosity of water, mugIs the viscosity of gas, a and b are model fitting parameters, a>0，b<0，SwThe water saturation.
Compared with the prior art, the method calculates the gas-water relative permeability ratio through a core analysis gas-water relative permeability experiment, constructs a new parameter containing a water saturation variable, establishes a fitting model of the gas-water relative permeability ratio and the new parameter containing the water saturation variable, substitutes a simplified flow splitting equation, and constructs a new prediction model of the water production rate and the water saturation. Through analysis of actual data, the new model has higher accuracy of water yield predicted by the physical method of the traditional oil layer, is suitable for different oil and gas reservoirs (different oil and gas reservoirs have different oil/gas-water viscosity ratios), has good universality and has good popularization and application values.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a graph of gas-water relative permeability versus water saturation for an example of the present invention;
FIG. 3 is a graph of gas-water relative permeability versus new parameters including water saturation variables for an embodiment of the present invention;
FIG. 4 is a comparison of water production rates predicted by the new model and the conventional reservoir physics methods in accordance with an embodiment of the present invention;
FIG. 5 is a comparison of water production rates predicted by the new model and the conventional reservoir physics methods for different gas-water viscosity ratios in an embodiment of the present invention; wherein (a) the gas-water viscosity ratio [ mu ] isw/μgIs 100, (b) the gas-water viscosity ratio muw/μgIs 10, (c) the gas-water viscosity ratio muw/μ g1, the gas-water viscosity ratio mu in (d)w/μgIs 0.1, and (e) the gas-water viscosity ratio muw/μgIs 0.01.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
Referring to fig. 1, the present invention comprises the steps of:
performing a gas-water relative permeability experiment on a core of a target reservoir in a research area to obtain data of water saturation, gas-phase relative permeability and water-phase relative permeability;
calculating a gas-water relative permeability ratio value by utilizing a core gas-water relative permeability experiment, and constructing a parameter containing a water saturation variable;
analyzing the data obtained in the step one and the parameters obtained in the step two, and establishing a fitting model of the gas-water relative permeability ratio and new parameters containing water saturation variables;
and step four, substituting the fitting model obtained in the step three into a simplified flow rate equation to complete the construction of a new prediction model of the water production rate and the water saturation.
Example (b):
referring to fig. 1, an embodiment of the present invention provides a method for constructing a tight sandstone reservoir water production rate prediction model, including the following steps:
step one, performing a gas-water relative permeability experiment on a rock core of a target reservoir in a research area to obtain water saturation, gas-phase relative permeability and water-phase relative permeability data. The experimental data of the gas-water relative permeability of the core are obtained according to the flow specified by the standard GB/T28912-2012 method for determining the relative permeability of the two-phase fluid in the rock.
And step two, calculating a gas-water relative permeability ratio value by using the gas-water relative permeability experiment of the core analysis in the step one, and constructing a new parameter containing a water saturation variable.
The calculation formula of the gas-water relative permeability value is as follows:
y＝krg/krw
in the formula, krgGas relative permeability, dimensionless; k is a radical ofrwIs the relative permeability of water, and is dimensionless.
The calculation formula for the new parameter containing the water saturation variable is as follows:
Figure BDA0002905430740000041
in the formula, SwIs the water saturation, decimal.
The results of the gas-water relative permeability experimental analysis of the core analysis are shown in table 1.
TABLE 1 gas-water relative permeability experiment analysis results table
Figure BDA0002905430740000042
Figure BDA0002905430740000051
Analyzing the two parameters, and establishing a fitting model of the gas-water relative permeability ratio and new parameters containing water saturation variables, wherein the model expression is as follows:
Figure BDA0002905430740000052
wherein y is the gas-water relative permeability ratio krg/krwNo dimension; x is a new parameter containing a water saturation variable
Figure BDA0002905430740000053
Dimensionless; swWater saturation, decimal; a. b is a model fitting parameter (a)>0，b<0)。
The water saturation S in Table 1wAnd a new parameter comprising a water saturation variable
Figure BDA0002905430740000054
Relative permeability k to gas-water respectivelyrg/krwStatistical fitting, both of which are fitted with exponential functions, with the highest correlation, 0.9488 and 0.9901, respectively, as shown in fig. 2 and 3, indicates that the new parameters including water saturation variables predict gas-water relative permeability values with higher accuracy than the gas-water relative permeability values of the water saturation predictions, and the fitted model expressions are:
Figure BDA0002905430740000055
R2＝0.9488
Figure BDA0002905430740000056
R2＝0.9901
and step four, substituting the fitting model into a simplified flow rate equation to construct a prediction model of new water production rate and water saturation. The simplified split flow equation is:
Figure BDA0002905430740000061
in the formula: f. ofwWater production rate, decimal; mu.swViscosity of water, mpa.s; mu.sgIs the viscosity of gas, mpa.s; k is a radical ofrgGas relative permeability, decimal; k is a radical ofrwRelative water permeability, decimal.
According to the traditional oil layer physics method, the established water yield model is as follows:
Figure BDA0002905430740000062
the water production rate model established according to the new parameters containing the water saturation variables, which is constructed according to the patent, has the following model expression:
Figure BDA0002905430740000063
in the actual data processing, the water yield curves calculated by the two methods are compared with the actual experimental results, and as shown in fig. 4, the gas-water viscosity ratio of the research area is 0.0157. It can be seen from the figure that the water yield predicted by the traditional fitting model of reservoir physics has a large error with the actual analysis data point, the predicted water yield is lower than the actual water yield under the condition of extra-high water content, and the predicted water yield is higher than the actual water yield under the condition of medium water content. The water yield calculated by the prediction model established by the method is well consistent with the analysis data points, the coincidence degree is obviously higher than that of the traditional oil layer physics method, and the prediction precision of the water yield is effectively improved.
Due to different oil and gas reservoirs with different water viscosity and oil/gas viscosity, the model method established by the method is contrastively analyzed with the water yield predicted by the traditional oil reservoir physics under the condition of different oil/gas-water viscosity ratios, as shown in figure 5. Drawing (A)Gas-water viscosity ratio muw/μgThe water content of the model method established by the method is 100, 10, 1, 0.1 and 0.01 in sequence, and as can be seen from the figure, the water content of the model method established by the method is well consistent with the analysis data points, and the coincidence degree is higher than that of the traditional oil layer physics method, so that the model method established by the method has good universality.
According to the method provided by the embodiment of the invention, the gas-water relative permeability ratio is calculated through a core analysis gas-water relative permeability experiment, a new parameter containing a water saturation variable is constructed, a fitting model of the gas-water relative permeability ratio and the new parameter containing the water saturation variable is established, and then a simplified flow splitting equation is substituted, so that a new prediction model of the water yield and the water saturation is constructed. Compared with the traditional oil layer physics method, the new model method effectively improves the prediction precision of the water yield, is suitable for different oil and gas reservoirs (different oil and gas-water viscosity ratios of the oil and gas reservoirs), has good universality and has good popularization and application prospects.

Claims (7)

1. A method for constructing a compact sandstone reservoir water yield prediction model is characterized by comprising the following steps:
performing a gas-water relative permeability experiment on a core of a target reservoir in a research area to obtain data of water saturation, gas-phase relative permeability and water-phase relative permeability;
calculating a gas-water relative permeability ratio value by utilizing a core gas-water relative permeability experiment, and constructing a parameter containing a water saturation variable;
analyzing the data obtained in the step one and the parameters obtained in the step two, and establishing a fitting model of the gas-water relative permeability ratio and new parameters containing water saturation variables;
and step four, substituting the fitting model obtained in the step three into a simplified flow rate equation to complete the construction of a new prediction model of the water production rate and the water saturation.
2. The method for constructing the compact sandstone reservoir water yield prediction model according to claim 1, wherein in the first step, the experimental data of the gas-water relative permeability of the core is obtained according to a flow specified in GB/T28912 + two-phase fluid relative permeability determination method in 2012 rocks.
3. The method for constructing the tight sandstone reservoir water yield prediction model according to claim 1, wherein in the second step, the calculation formula of the gas-water relative permeability ratio is as follows:
y＝krg/krw
in the formula, krgIs gas relative permeability, krwRelative water permeability.
4. The method for constructing the tight sandstone reservoir water production rate prediction model according to claim 1, wherein in the second step, a new parameter calculation formula containing a water saturation variable is as follows:
Figure FDA0002905430730000011
in the formula, SwThe water saturation.
5. The method for constructing the tight sandstone reservoir water yield prediction model according to claim 1, wherein in the third step, the expression of the fitting model is as follows:
Figure FDA0002905430730000021
wherein y is the gas-water relative permeability ratio krg/krwX is a new parameter containing a water saturation variable
Figure FDA0002905430730000022
SwFor water saturation, a and b are model fitting parameters, a>0，b<0。
6. The method for constructing the tight sandstone reservoir water yield prediction model according to claim 1, wherein in the fourth step, the simplified equation expression of the flow rate is as follows:
Figure FDA0002905430730000023
in the formula: f. ofwTo yield water, muwIs the viscosity of water, mugIs the viscosity of gas, krgIs gas relative permeability, krwRelative water permeability.
7. The method for constructing the tight sandstone reservoir water yield prediction model according to claim 1, wherein in the fourth step, the expression of the prediction model with the new water yield and the new water saturation is as follows:
Figure FDA0002905430730000024
in the formula: mu.swIs the viscosity of water, mugIs the viscosity of gas, a and b are model fitting parameters, a>0，b<0，SwThe water saturation.
CN202110069991.9A 2021-01-19 2021-01-19 Method for constructing compact sandstone reservoir water yield prediction model Pending CN112926275A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110069991.9A CN112926275A (en) 2021-01-19 2021-01-19 Method for constructing compact sandstone reservoir water yield prediction model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110069991.9A CN112926275A (en) 2021-01-19 2021-01-19 Method for constructing compact sandstone reservoir water yield prediction model

Publications (1)

Publication Number Publication Date
CN112926275A true CN112926275A (en) 2021-06-08

Family

ID=76163433

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110069991.9A Pending CN112926275A (en) 2021-01-19 2021-01-19 Method for constructing compact sandstone reservoir water yield prediction model

Country Status (1)

Country Link
CN (1) CN112926275A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114139470A (en) * 2021-11-29 2022-03-04 长江大学 Method for calculating gas relative permeability of tight sandstone gas reservoir

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104504457A (en) * 2014-12-04 2015-04-08 中国石油大港油田勘探开发研究院 Water production gas well productivity prediction method
CN107725044A (en) * 2017-09-27 2018-02-23 中国石油天然气集团公司 The method of Sandstone Gas Bearing reservoir producing water ratio prediction based on array induction, laterolog
CN110895637A (en) * 2018-09-07 2020-03-20 中国石油化工股份有限公司 Method for controlling reserves of single well of water producing gas well
CN111425191A (en) * 2020-03-25 2020-07-17 西南石油大学 Method and device for calculating gas-water relative permeability of tight reservoir
CN111476448A (en) * 2019-01-24 2020-07-31 中国石油化工股份有限公司 Dynamic yield splitting method for multilayer commingled production vertical well

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104504457A (en) * 2014-12-04 2015-04-08 中国石油大港油田勘探开发研究院 Water production gas well productivity prediction method
CN107725044A (en) * 2017-09-27 2018-02-23 中国石油天然气集团公司 The method of Sandstone Gas Bearing reservoir producing water ratio prediction based on array induction, laterolog
CN110895637A (en) * 2018-09-07 2020-03-20 中国石油化工股份有限公司 Method for controlling reserves of single well of water producing gas well
CN111476448A (en) * 2019-01-24 2020-07-31 中国石油化工股份有限公司 Dynamic yield splitting method for multilayer commingled production vertical well
CN111425191A (en) * 2020-03-25 2020-07-17 西南石油大学 Method and device for calculating gas-water relative permeability of tight reservoir

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
杨冬;张海荣;吴一雄;张恒荣;谭伟;: "电缆地层测试资料在低渗气藏产能预测中的应用", 中国海上油气, vol. 30, no. 06, pages 109 - 114 *
王谦;谭茂金;石玉江;李高仁;程相志;罗伟平;: "径向基函数神经网络法致密砂岩储层相对渗透率预测与含水率计算", 石油地球物理勘探, vol. 55, no. 04, pages 864 - 872 *
窦祥骥;林恺;何岩峰;王相;刘成;赵虹宇;: "束缚水饱和度应力敏感致密气井动态特征研究", 常州大学学报(自然科学版), vol. 30, no. 05, pages 75 - 80 *
胡俊坤;李晓平;宋代诗雨;: "水驱气藏动态储量计算新方法", 天然气地球科学, vol. 24, no. 03, pages 628 - 632 *
苟燕;刘华勋;高树生;胡志明;薛蕙;周慧智;: "水驱气藏的水驱特征曲线与应用效果分析", 石油钻采工艺, vol. 35, no. 03, pages 63 - 65 *
鹿克峰;蔡华;王理;石美雪;何贤科;范洪军;: "中国东海气区初始产水评价图版的建立", 天然气工业, vol. 39, no. 05, pages 63 - 70 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114139470A (en) * 2021-11-29 2022-03-04 长江大学 Method for calculating gas relative permeability of tight sandstone gas reservoir

Similar Documents

Publication Publication Date Title
CN103452547B (en) The analysis and processing method of afterflow data and system in well test data
CN104134101A (en) Low-seepage reservoir natural gas productivity prediction method
CN104899411B (en) A kind of reservoir productivity prediction model method for building up and system
CN108661631B (en) Yield prediction method
CN110992200A (en) Shale gas well staged fracturing effect evaluation and yield prediction method based on random forest
CN111425193A (en) Reservoir compressibility evaluation method based on clustering analysis logging rock physical facies division
CN107808032B (en) Shale gas reservoir water saturation calculation method and system
CN114201932B (en) Compact oil reservoir fracturing well test simulation method under complex condition
US10385659B2 (en) Evaluation of production performance from a hydraulically fractured well
CN112392464A (en) Method for calculating reservoir water production rate based on conventional logging information
CN112926275A (en) Method for constructing compact sandstone reservoir water yield prediction model
CN112145165B (en) Microcrack-pore type reservoir dynamic and static permeability conversion method
CN111912761B (en) Method for accurately representing oil-water phase permeability relation curve in ultra-high water-cut period
CN110909466B (en) Method for improving calculation accuracy of average capillary pressure curve of differential pore sorting reservoir
CN111950111B (en) Dynamic analysis method suitable for open-bottom carbonate reservoir
CN112257349A (en) Method for judging whether compact sandstone movable water-gas reservoir gas well has development value
CN108664678B (en) Yield prediction method
CN112035993A (en) Method for testing and evaluating carbonate reservoir with constant pressure at bottom
CN111488666A (en) Gas reservoir horizontal well productivity prediction model establishing and predicting method and device
CN111827995B (en) Nuclear magnetic resonance-based calculation method for permeability of overpressure conglomerate reservoir
CN110500082B (en) Method for determining gas invasion dangerous time of well cementation cement slurry
CN114047087A (en) Screen pipe erosion life prediction method based on reference well and experimental data
CN111241651A (en) Method for making novel water-drive plate
CN112727433A (en) Drilling parameter optimization method
CN110909311A (en) Method for calculating gas content of thin coal seam

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