CN117648523A - 一种有水气藏动态储量及水侵常数计算方法 - Google Patents

一种有水气藏动态储量及水侵常数计算方法 Download PDF

Info

Publication number
CN117648523A
CN117648523A CN202410117121.8A CN202410117121A CN117648523A CN 117648523 A CN117648523 A CN 117648523A CN 202410117121 A CN202410117121 A CN 202410117121A CN 117648523 A CN117648523 A CN 117648523A
Authority
CN
China
Prior art keywords
water
unit
dimensionless
gas
constant
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
CN202410117121.8A
Other languages
English (en)
Other versions
CN117648523B (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.)
Chengdu Yingwoxin Technology Co ltd
Original Assignee
Chengdu Yingwoxin Technology 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 Chengdu Yingwoxin Technology Co ltd filed Critical Chengdu Yingwoxin Technology Co ltd
Priority to CN202410117121.8A priority Critical patent/CN117648523B/zh
Publication of CN117648523A publication Critical patent/CN117648523A/zh
Application granted granted Critical
Publication of CN117648523B publication Critical patent/CN117648523B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明涉及一种有水气藏动态储量及水侵常数计算方法,其技术方案是考虑已有常规水驱特征曲线高含水期与特高含水期拟合不准确,推导新型水驱特征曲线方程,给定动态储量G和水侵常数B的上下界,通过调整动态储量G和水侵常数B来拟合实际的水气比,得到动态储量G和水侵常数B,并绘制水气比拟合曲线图和气水相对渗透率拟合曲线图。本发明特点在于能反映高含水期水驱特征曲线上翘的现象,拟合更加准确,可用于计算动态储量以及水侵常数。

Description

一种有水气藏动态储量及水侵常数计算方法
技术领域
本发明涉及一种有水气藏动态储量及水侵常数计算方法,属于气藏工程领域。
背景技术
四川盆地天然气开发历史悠久,资源丰富,其中绝大部分为有水气藏,地质情况极为复杂,分别呈现出边水,底水,无统一气水界面等多种水体分布现象,气藏开发过程中水侵危害时常严重,水驱气藏气井见水后,产气产水规律复杂,使得水驱气藏开发难度显著增大。就目前水驱气藏而言,缺乏相应成熟有效的动态分析方法,需要更加准确的进行水侵危害预测,优化开发策略。
目前,已有常规水驱特征曲线在高含水期与特高含水期拟合不准确,只有在含水饱和度处于中间阶段时才能拟合准确。该新型水驱特征曲线方程能反映高含水期水驱特征曲线上翘的现象,拟合更加准确,可用于计算动态储量以及水侵常数,这对现在开发有一定指导意义。
发明内容
本发明的目的是:本发明针对现有技术存在的问题推导新型水驱特征曲线方程,该方程能反映高含水期水驱特征曲线上翘的现象使得生产历史拟合更加准确,也可以更加准确的预测水侵危害以及计算动态储量和水侵常数。
为达到上述目的,本发明提供了一种有水气藏动态储量及水侵常数计算方法,该方法包括以下步骤:
S100、推导新型水驱特征曲线方程:
S101、将存水体积系数的定义方程和/>方程依次代入存水体积与天然气地下体积之比方程/>推导含水饱和度方程:
式中,B为水侵常数,定义为,单位为无量纲;S w为含水饱和度,单位为无量纲;S wi为原始条件下含水饱和度,单位为无量纲;R为采出程度,单位为无量纲;B gi为原始条件下天然气体积系数,单位为无量纲;W e为水侵量,单位为方;W p为累产水,单位为方;/>为存水体积系数,单位为无量纲;B w是地层水体积系数,单位为无量纲;G是动态储量,单位为万方;R B为采出程度的水侵常数次方,根据定义有采出程度的水侵常数次方R B等于存水体积系数/>
S102、将S101中推导出的含水饱和度方程和气水相对渗透率比值方程/>依次代入水气比方程/>推导新型水驱特征曲线方程:
其中,B g为天然气的体积系数,定义为,单位为无量纲;p sc为地面压力,单位为MPa;p为地层压力,单位为MPa;T为地层温度,单位为K;T sc为地面温度,单位为K;Z为地层条件下的气体偏差因子,单位为无量纲;Z sc为地面条件下的气体偏差因子,单位为无量纲;S w为含水饱和度,单位为无量纲;S wi为原始条件下含水饱和度,单位为无量纲;R为采出程度,单位为无量纲;WGR是水气比,单位为方/万方;/>是天然气黏度,单位为mPa·s;是地层水黏度,单位为mPa·s;B w是地层水体积系数,单位为无量纲;K rw是岩心水相相对渗透率,单位为无量纲;K rg是岩心气相相对渗透率,单位为无量纲;nb为与储层和流体物性有关的常数,由气水相对渗透率曲线拟合得到,单位为无量纲;R B为采出程度的水侵常数次方,根据定义有采出程度的水侵常数次方R B等于存水体积系数/>
S200、准备静态资料和生产历史数据,静态资料包括气水相对渗透率曲线、原始含水饱和度、天然气黏度、地层水黏度、地层温度、地面温度、地面压力;生产历史数据包括日产气量、日产水量、地层压力;
S300、根据S102中的新型水驱特征曲线方程,进行生产历史拟合的思路为:
S301、给定动态储量G,水侵常数B两个参数的上下界限,其中动态储量G的上界为20000万方,动态储量G的下界为1000万方,水侵常数B的上界为20,水侵常数B的下界为1,上下界限内随机给定初值;采出程度R由累产G p除以动态储量G计算,其中的动态储量G为可变参数;
S302、结合S301中的动态储量G和水侵常数B代入新型水驱特征曲线方程可以算出预测的水气比;其中,在计算地层条件下的天然气体积系数B g时,地层条件下的气体偏差因子Z由地层压力根据DAK法计算;
S303、再重复S301~S302的过程,以生产数据中的一天为一个周期,计算得到整个生产阶段的水气比;
S304、在参数的上下限内调整参数,采用自动拟合方法拟合实际的水气比,寻求预测值与实际值的最佳拟合,水气比拟合的收敛条件表示为:
其中,E为偏差值;m为迭代次数;lnWGR i(G,B)为预测的水气比取对数的值,单位为方/万方;lnWGR i为实际的水气比取对数的值,单位为方/万方;
S305、对有水气藏生产历史进行拟合,拟合成功后,停止调整动态储量G、水侵常数B两个参数,以拟合最优时的动态储量G、水侵常数B为计算结果。
与现有方法相比,本发明具有以下有益效果:(1)推导新型水驱特征曲线公式,能够反映高含水期水驱特征曲线上翘的现象,拟合更加准确;(2)能够更加准确的计算有水气藏的动态储量以及水侵常数;
附图说明
图1是本方法技术路线图;
图2是某有水气藏气水相对渗透率拟合曲线图;
图3是某有水气藏水气比拟合曲线图。
具体实施方式
本发明提供了一种有水气藏动态储量及水侵常数计算方法,图1为本方法技术路线图,针对国内某有水气藏,实施以下步骤:
第一步:推导新型水驱特征曲线方程:将存水体积系数的定义方程和/>方程依次带入存水体积与天然气地下体积之比方程/>推导含水饱和度方程:
式中,B为水侵常数,定义为,单位为无量纲;S w为含水饱和度,单位为无量纲;S wi为原始条件下含水饱和度,单位为无量纲,取0.4996;R为采出程度,单位为无量纲;B gi为原始条件下天然气体积系数,单位为无量纲;W e为水侵量,单位为方;W p为累产水,单位为方;/>为存水体积系数,单位为无量纲;B w是地层水体积系数,单位为无量纲;G是动态储量,单位为万方;R B为采出程度的水侵常数次方,根据定义有采出程度的水侵常数次方R B等于存水体积系数/>;再将推导出的含水饱和度方程和气水相对渗透率比值方程/>依次代入水气比方程推导新型水驱特征曲线方程:
其中,B g为天然气的体积系数,定义为,单位为无量纲;p sc为地面压力,单位为MPa,取0.101;p为地层压力,单位为MPa;T为地层温度,单位为K,取323.15;T sc为地面温度,单位为K,取293.15;Z为地层条件下的气体偏差因子,单位为无量纲;Z sc为地面条件下的气体偏差因子,单位为无量纲,取1;S w为含水饱和度,单位为无量纲;S wi为原始条件下含水饱和度,单位为无量纲,取0.4996;R为采出程度,单位为无量纲;WGR是水气比,单位为方/万方;/>是天然气黏度,单位为mPa·s,取0.02;/>是地层水黏度,单位为mPa·s,取0.4;B w是地层水体积系数,单位为无量纲,取1;K rw是岩心水相相对渗透率,单位为无量纲;K rg是岩心气相相对渗透率,单位为无量纲;nb为与储层和流体物性有关的常数,由气水相对渗透率曲线拟合得到,单位为无量纲;n为2200,b为6,拟合效果如图2所示;
第二步:准备静态资料和生产历史数据,静态资料包括气水相对渗透率曲线、原始含水饱和度、天然气黏度、地层水黏度、地层温度、地面温度、地面压力;生产历史数据包括日产气量、日产水量、地层压力;
第三步:对生产数据进行历史拟合:取动态储量G为100000万方,水侵常数B为1,其中动态储量G的上界为200000万方,动态储量G的下界为1000万方,水侵常数B的上界为20,水侵常数B的下界为1;结合上述的动态储量G和水侵常数B代入新型水驱特征曲线方程算出整个生产阶段预测的水气比;其中,在计算地层条件下的天然气体积系数B g时,地层条件下的气体偏差因子Z由地层压力根据DAK法计算;在参数的上下限内调整参数,采用自动拟合方法拟合实际的水气比,寻求预测值与实际值的最佳拟合,水气比拟合的收敛条件表示为:
其中,E为偏差值;m为迭代次数;lnWGR i(G,B)为预测的水气比取对数的值,单位为方/万方;lnWGR i为实际的水气比取对数的值,单位为方/万方;
拟合成功后,停止调整动态储量G、水侵常数B两个参数,拟合最优时的动态储量G为91800万方、水侵常数B为1.7,以此时的值作为计算结果,拟合效果如图3所示。
最后所应说明的是:以上实施例仅用以说明而非限制本发明的技术方案,尽管参照上述实施例对本发明进行了详细说明,本领域的普通技术人员应该理解:依然可以对本发明进行修改或者等同替换,而不脱离本发明的精神和范围的任何修改或局部替换,其均应涵盖在本发明的权利要求范围当中。

Claims (1)

1.一种有水气藏动态储量及水侵常数计算方法,其特征在于,该方法包括以下步骤:
S100、推导新型水驱特征曲线方程,步骤如下:
S101、将存水体积系数的定义方程和/>方程依次代入存水体积与天然气地下体积之比方程/>推导含水饱和度方程:
式中,B为水侵常数,定义为,单位为无量纲;S w为含水饱和度,单位为无量纲;S wi为原始条件下含水饱和度,单位为无量纲;R为采出程度,单位为无量纲;B gi为原始条件下天然气体积系数,单位为无量纲;W e为水侵量,单位为方;W p为累产水,单位为方;/>为存水体积系数,单位为无量纲;B w是地层水体积系数,单位为无量纲;G是动态储量,单位为万方;R B为采出程度的水侵常数次方,根据定义有采出程度的水侵常数次方R B等于存水体积系数/>
S102、将S101中推导出的含水饱和度方程和气水相对渗透率比值方程/>依次代入水气比方程/>推导新型水驱特征曲线方程:
其中,B g为天然气的体积系数,定义为,单位为无量纲;p sc为地面压力,单位为MPa;p为地层压力,单位为MPa;T为地层温度,单位为K;T sc为地面温度,单位为K;Z为地层条件下的气体偏差因子,单位为无量纲;Z sc为地面条件下的气体偏差因子,单位为无量纲;S w为含水饱和度,单位为无量纲;S wi为原始条件下含水饱和度,单位为无量纲;R为采出程度,单位为无量纲;WGR是水气比,单位为方/万方;/>是天然气黏度,单位为mPa·s;是地层水黏度,单位为mPa·s;B w是地层水体积系数,单位为无量纲;K rw是岩心水相相对渗透率,单位为无量纲;K rg是岩心气相相对渗透率,单位为无量纲;nb为与储层和流体物性有关的常数,由气水相对渗透率曲线拟合得到,单位为无量纲;R B为采出程度的水侵常数次方,根据定义有采出程度的水侵常数次方R B等于存水体积系数/>
S200、准备静态资料和生产历史数据,静态资料包括气水相对渗透率曲线、原始含水饱和度、天然气黏度、地层水黏度、地层温度、地面温度、地面压力;生产历史数据包括日产气量、日产水量、地层压力;
S300、根据S102中的新型水驱特征曲线方程,进行生产历史拟合的思路为:
S301、给定动态储量G,水侵常数B两个参数的上下界限,其中动态储量G的上界为20000万方,动态储量G的下界为1000万方,水侵常数B的上界为20,水侵常数B的下界为1,上下界限内随机给定初值;采出程度R由累产G p除以动态储量G计算,其中的动态储量G和水侵常数B均为可变参数;
S302、结合S301中的动态储量G和水侵常数B代入新型水驱特征曲线方程可以算出预测的水气比;其中,在计算地层条件下的天然气体积系数B g时,地层条件下的气体偏差因子Z由地层压力根据DAK法计算;
S303、再重复S301~S302的过程,以生产数据中的一天为一个周期,计算得到整个生产阶段的水气比;
S304、在参数的上下限内调整参数,采用自动拟合方法拟合实际的水气比,寻求预测值与实际值的最佳拟合,水气比拟合的收敛条件表示为:
其中,E为偏差值;m为迭代次数;lnWGR i(G,B)为预测水气比取对数的值,单位为方/万方;lnWGR i为实际的水气比取对数的值,单位为方/万方;
S305、对有水气藏生产历史进行拟合,拟合成功后,停止调整动态储量G、水侵常数B两个参数,以拟合最优的动态储量G、水侵常数B为计算结果。
CN202410117121.8A 2024-01-29 2024-01-29 一种有水气藏动态储量及水侵常数计算方法 Active CN117648523B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410117121.8A CN117648523B (zh) 2024-01-29 2024-01-29 一种有水气藏动态储量及水侵常数计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410117121.8A CN117648523B (zh) 2024-01-29 2024-01-29 一种有水气藏动态储量及水侵常数计算方法

Publications (2)

Publication Number Publication Date
CN117648523A true CN117648523A (zh) 2024-03-05
CN117648523B CN117648523B (zh) 2024-04-05

Family

ID=90049942

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410117121.8A Active CN117648523B (zh) 2024-01-29 2024-01-29 一种有水气藏动态储量及水侵常数计算方法

Country Status (1)

Country Link
CN (1) CN117648523B (zh)

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104951636A (zh) * 2014-03-24 2015-09-30 中国石油化工股份有限公司 特高含水期甲型水驱曲线非线性关系表征方法
CN105225019A (zh) * 2015-11-11 2016-01-06 四川建筑职业技术学院 特高含水期累积产油量及采收率预测方法
RU2704400C1 (ru) * 2018-06-14 2019-10-28 Петрочайна Компани Лимитед Способ и устройство для прогнозирования изменения скорости нарастания обводненности в нефтяном пласте с водонапорным режимом
CN110939432A (zh) * 2019-12-06 2020-03-31 中海石油(中国)有限公司 一种特高含水期油藏含水率预测方法
CN110969307A (zh) * 2019-12-06 2020-04-07 中海石油(中国)有限公司 一种特高含水期油藏累产油量预测方法
CN111927411A (zh) * 2020-09-24 2020-11-13 西南石油大学 一种智能化有水气藏水侵跟踪和预警方法
CN112307634A (zh) * 2020-11-04 2021-02-02 中国石油大学(华东) 一种特高含水期油田可采储量标定方法、系统及存储介质
CN112696194A (zh) * 2019-10-22 2021-04-23 中国石油化工股份有限公司 一种确定特高含水油藏可动油饱和度的方法
CN113027399A (zh) * 2021-03-18 2021-06-25 中国石油大学(华东) 一种基于微观流动模拟获取高含水区块水驱曲线方法
CN113236207A (zh) * 2021-07-13 2021-08-10 西南石油大学 一种强非均质性储层中产水气井的定产量递减预测方法
CN114742330A (zh) * 2022-06-13 2022-07-12 西南石油大学 一种高含硫有水气藏水封气量预测方法
CN115422493A (zh) * 2022-09-13 2022-12-02 西南石油大学 一种含砂体叠置有水气藏动态储量计算方法
CN115906677A (zh) * 2021-08-18 2023-04-04 中国石油天然气股份有限公司 一种非均质性边水气藏水体活跃程度预测方法
CN116776041A (zh) * 2023-07-13 2023-09-19 西南石油大学 一种超高压边底水凝析气藏水侵程度的图版计算方法
CN117291059A (zh) * 2023-11-23 2023-12-26 成都英沃信科技有限公司 一种高温高压有水气藏可采储量及水侵参数计算方法

Patent Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104951636A (zh) * 2014-03-24 2015-09-30 中国石油化工股份有限公司 特高含水期甲型水驱曲线非线性关系表征方法
CN105225019A (zh) * 2015-11-11 2016-01-06 四川建筑职业技术学院 特高含水期累积产油量及采收率预测方法
RU2704400C1 (ru) * 2018-06-14 2019-10-28 Петрочайна Компани Лимитед Способ и устройство для прогнозирования изменения скорости нарастания обводненности в нефтяном пласте с водонапорным режимом
CN112696194A (zh) * 2019-10-22 2021-04-23 中国石油化工股份有限公司 一种确定特高含水油藏可动油饱和度的方法
CN110939432A (zh) * 2019-12-06 2020-03-31 中海石油(中国)有限公司 一种特高含水期油藏含水率预测方法
CN110969307A (zh) * 2019-12-06 2020-04-07 中海石油(中国)有限公司 一种特高含水期油藏累产油量预测方法
US20220010674A1 (en) * 2020-09-24 2022-01-13 Southwest Petroleum University Intelligent water invasion tracking and early warning method for water-gas reservoirs
CN111927411A (zh) * 2020-09-24 2020-11-13 西南石油大学 一种智能化有水气藏水侵跟踪和预警方法
CN112307634A (zh) * 2020-11-04 2021-02-02 中国石油大学(华东) 一种特高含水期油田可采储量标定方法、系统及存储介质
CN113027399A (zh) * 2021-03-18 2021-06-25 中国石油大学(华东) 一种基于微观流动模拟获取高含水区块水驱曲线方法
CN113236207A (zh) * 2021-07-13 2021-08-10 西南石油大学 一种强非均质性储层中产水气井的定产量递减预测方法
CN115906677A (zh) * 2021-08-18 2023-04-04 中国石油天然气股份有限公司 一种非均质性边水气藏水体活跃程度预测方法
CN114742330A (zh) * 2022-06-13 2022-07-12 西南石油大学 一种高含硫有水气藏水封气量预测方法
US20230417138A1 (en) * 2022-06-13 2023-12-28 Southwest Petroleum University Method for predicting amount of water-sealed gas in high-sulfur water-bearing gas reservoir
CN115422493A (zh) * 2022-09-13 2022-12-02 西南石油大学 一种含砂体叠置有水气藏动态储量计算方法
CN116776041A (zh) * 2023-07-13 2023-09-19 西南石油大学 一种超高压边底水凝析气藏水侵程度的图版计算方法
CN117291059A (zh) * 2023-11-23 2023-12-26 成都英沃信科技有限公司 一种高温高压有水气藏可采储量及水侵参数计算方法

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
TAN XIAOHUA 等: "water-drive characteristic curve", 《NATURAL GAS INDUSTRY B》, vol. 08, no. 04, 31 August 2021 (2021-08-31), pages 353 - 358 *
XIAOBING HAN 等: "Water invasion performance of complex fracture-vuggy gas reservoirs based on classification modeling", 《ADVANCES IN GEO-ENERGY RESEARCH》, vol. 5, no. 2, 11 February 2021 (2021-02-11), pages 222 - 232 *
刘俊丰;徐晓玲;柴小颖;韩道明;张勇年: "气藏动态储量计算新方法", 《化工管理》, no. 13, 1 May 2020 (2020-05-01), pages 89 - 90 *
李珂;杨莉;张迎春;皮建;: "一种新型水驱特征曲线的推导及应用", 断块油气田, no. 06, 26 April 2019 (2019-04-26), pages 13 - 16 *
王继强;石成方;纪淑红;李冠林;陈映桥;: "特高含水期新型水驱特征曲线", 石油勘探与开发, no. 06, 21 August 2017 (2017-08-21), pages 955 - 960 *
胡俊坤;李晓平;宋代诗雨;: "水驱气藏动态储量计算新方法", 天然气地球科学, no. 03, 10 June 2013 (2013-06-10), pages 628 - 632 *
范海军;朱学谦;: "高含水期油田新型水驱特征曲线的推导及应用", 断块油气田, no. 01, 29 January 2016 (2016-01-29), pages 105 - 108 *
赵翔: "S边水气藏动态储量及水侵特征研究", 《中国优秀硕士学位论文全文数据库-工程科技I辑》, no. 07, 15 July 2019 (2019-07-15), pages 019 - 370 *
陈存良;王相;何芬;刘超;崔龙涛;: "高含水期新型水驱特征曲线研究与应用", 常州大学学报(自然科学版), no. 05, 28 September 2020 (2020-09-28), pages 76 - 80 *

Also Published As

Publication number Publication date
CN117648523B (zh) 2024-04-05

Similar Documents

Publication Publication Date Title
CN110334431A (zh) 一种低渗透致密气藏单井控制储量计算及剩余气分析方法
CN109958404B (zh) 底水油藏底水锥进及堵水调剖剂参数计算方法
CN108444884B (zh) 基于Harris模型的油相相对渗透率曲线表征方法
CN111810101B (zh) 一种水驱油藏动态分析方法及装置
CN106150463B (zh) 一种砾岩油藏聚合物驱注入压力升幅确定方法
CN110410045B (zh) 一种针对特高含水期油藏的堵调方法
CN109858177B (zh) 一种基于拟流线法的底水油藏水平井水驱数值模拟方法
CN109958413B (zh) 一种特高含水期油藏动态流动单元划分方法
CN111456709A (zh) 一种基于测井曲线的水平井多级压裂分段分簇方法
CN106321051A (zh) 一种用于优化多段压裂水平井网络裂缝参数的方法
CN110374562B (zh) 一种油藏堵调方法
CN117648523B (zh) 一种有水气藏动态储量及水侵常数计算方法
CN107288595B (zh) 一种注水利用率的评价方法
CN116629165B (zh) 储层压裂改造区和未改造区参数反演方法、系统及设备
CN112443323A (zh) 基于单位压降采出程度的水驱气藏生产动态分析方法
CN113326465B (zh) 油藏开发动态分析方法及装置
CN112949901A (zh) 页岩油藏水平井准自然能量开发方法及系统
CN113468476B (zh) 一种聚合物驱增油量预测方法及装置
CN112796749B (zh) 一种基于毛管数的油藏水驱数值模拟方法
CN104747154A (zh) 一种利用驱油效率比提高蒸汽驱剩余油研究精度的方法
CN111199069B (zh) 基于门限回归理论的面板堆石坝坝顶沉降经验预测方法
CN114169204A (zh) 一种用于海上油气田开发生产的防砂时机确定方法
CN112709547A (zh) 一种产水气井堵水时机的判别方法
CN114429085A (zh) 一种用于分析缝洞型油藏流体势的方法及系统
CN117648888B (zh) 一种有水气藏气水过渡带推水保压延缓气井见水的方法

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