CN110219639B - 试油抽汲作业的自动计量的方法 - Google Patents
试油抽汲作业的自动计量的方法 Download PDFInfo
- Publication number
- CN110219639B CN110219639B CN201910345688.XA CN201910345688A CN110219639B CN 110219639 B CN110219639 B CN 110219639B CN 201910345688 A CN201910345688 A CN 201910345688A CN 110219639 B CN110219639 B CN 110219639B
- Authority
- CN
- China
- Prior art keywords
- time
- liquid level
- period
- time point
- liquid
- 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.)
- Active
Links
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
- E21B49/08—Obtaining fluid samples or testing fluids, in boreholes or wells
Abstract
本发明涉及试油试采技术领域,是一种试油抽汲作业的自动计量的方法,包括S1:根据设置的时间段,采集该时间段内各时间点的液面高度,建立基于时间‑液面高度的第一数组,S2:根据各时间点对应的液面高度,以每抽子液体装满计量罐的时间设置为一个叠加周期;S3:对第一数组中的一个叠加周期T内的液面高度进行降噪处理;S4:计算一个叠加周期内液面上涨斜率;S5:判断计量罐的进液时段;S6:求解该一个叠加周期T内的平均综合含水率、液面高度差、每抽子的出液量和每抽子的产油量。本发明实现自动计算每抽子的产量、综合含水,通过判断抽汲作业每抽子的开始时间、结束时间、每抽汲一次液面准确的上升高度数据以及准确的实时含水率数据,有效提高了检测的准确率。
Description
技术领域
本发明涉及试油试采技术领域,是一种试油抽汲作业的自动计量的方法。
背景技术
目前,抽汲作井业是非自喷井试油求产的一种方式,目的是通过钢丝举升的方式携带井筒内流体,降低井筒内液柱压力,使井筒内液柱压力低于地层压力。在抽汲过程中由于压力差的作用,地层流体不断进入筒,通过抽子将液体抽至地面试井罐,地面计量后评价地层产液能力。
当前计量方式是:判断抽汲作业每抽子的开始时间、结束时间、每抽汲一次液面准确的上升高度数据以及准确的实时含水率数据,但目前存在问题:
1、液面高度存在误差:不进行抽汲作业时,液面高度应不上升,但是由于仪器本身数据漂移及天气因素(刮风)引起不进行抽汲作业时液面高度上升,软件会判断出现了抽汲,导致数据失真,计算结果错误。
2、综合含水存在误差:不进行抽汲作业时,综合含水应为零,但由于出液管线内存在残液引起不进行抽汲作业时综合含水示数不为零,软件判断失真,计算结果不准。
根据液面高度或综合含水确定抽汲作业每抽子的起、止时间存在一定困难,因此亟需一种在试油抽汲作业中能准确判断每抽子的起始时间并且能够准确计量抽汲作业每抽子的抽出液量(出水量、出油量)及其综合含水方法。
发明内容
本发明提供了一种试油抽汲作业的自动计量的方法,克服了上述现有技术之不足,其能有效解决现有的试油层的计算方式存在费工费时以及单点计量结果误差大的问题。
本发明的技术方案之一是通过以下措施来实现的:一种试油抽汲作业的自动计量的方法,包括以下步骤:
S1:设置自动计量的时间段,采集该时间段内各时间点的计量罐的液面高度,建立基于时间-液面高度的第一数组arr1,记为(№11,t11,h11)、(№12,t12,h12)……(№1i,t1i,h1i)……(№1n,t1n,h1n),其中,i取值为(1≤i≤n),№为序号,t为时间,h为液面高度;
S2:根据各时间点对应的计量罐的液面高度变化情况,测量每抽子液体装满计量罐的时间,并设置该时间段为一个叠加周期T;
S3:对第一数组中的一个叠加周期T内的液面高度进行降噪处理:取数组arr1中任意一点(№1i,t1i,h1i),以一个叠加周期T为步长,分别对一个叠加周期T内各时间点的液面的高度数据进行叠加去噪,形成第二数组arr2,记为(№21,t21,h21)、(№22,t22,h22)……(№2i,t2i,h2i)……(№2n,t2n,h2n),具体包括:
因降噪后的时间点t2i与降噪前的时间点t1i是指同一个时间点,则根据t2i=t1i,该一个叠加周期T内任一时间点降噪后的液面高度h2i计算公式如下:
其中,n为t1i-T/8至t1i+T/8之间对应的时间点的个数;
S4:计算液面上升斜率ai:以一个叠加周期T为步长,利用最小二乘法,求一个叠加周期T内各时间点的液面上升斜率ai,根据该一个叠加周期T内各时间点的液面上升斜率ai的计算结果建立第三数组arr3,记为(№31,t31,a1)、(№32,t32,a2)……(№3i,t3i,ai)……(№3n,t3n,an),具体包括:
因时间t3i与时间t2i是指同一个时间点,则根据t3i=t2i,计算一个叠加周期T内任一时间点的液面上涨斜率ai,公式如下:
S5:比较数组arr3的各项,根据公式(2)计算ai-ai-1,判断计量罐的进液时段;
S51:若ai-ai-1>0,则表示液面高度上涨,流体进入计量罐;
S52:若ai-ai-1≤0,则表示液面高度未上涨,无流体进入计量罐;
当ai-ai-1>0时,标记该时间点i对应的液面上升的起始时间为Tksi,标记该时间点i对应的液面高度为Hksi;当ai-ai-1≤0时,该时间点i对应的液面上升的结束点时间Tjsi,标记该时间点i的液面高度Hjsi;
其中,n为Tksi至Tjsi时间段内的时间点的个数,P为Tksi至Tjsi时间段内的各时间点的实时瞬时含水率;
S62:计算Tksi至Tjsi时间段的液面高度差,计算公式如下:
Hgc=Hjsi-Hksi (4)
S63:计算每抽子的出液量,计算公式如下:
V液=Hgc*Vjlg (5)
其中,Vjlg为计量罐的体积系数,根据实际使用的计量罐的大小进行设置;
S64:每抽子的出油量,计算公式如下:
本发明实现自动计算每抽子的产量、综合含水,通过判断抽汲作业每抽子的开始时间、结束时间、每抽汲一次液面准确的上升高度数据以及准确的实时含水率数据。解决了现有的计量检测不准确的问题,有效提高了检测的准确率。
附图说明
附图1为本发明实施例一的方法流程图。
附图2为本发明实施例二的在设定时间段内的液面高度的曲线图。
附图3为本发明实施例二的在设定时间段内的含水率的曲线图。
具体实施方式
本发明不受下述实施例的限制,可根据本发明的技术方案与实际情况来确定具体的实施方式。
下面结合实施例及附图对本发明作进一步描述:
实施例一:如附图1所示,一种试油抽汲作业的自动计量的方法,包括以下步骤:
S1:设置自动计量的时间段,采集该时间段内各时间点的计量罐的液面高度,建立基于时间-液面高度的第一数组arr1,记为(№11,t11,h11)、(№12,t12,h12)……(№1i,t1i,h1i)……(№1n,t1n,h1n),其中,i取值为(1≤i≤n),№为序号,t为时间,h为液面高度;
S2:根据各时间点对应的计量罐的液面高度变化情况,测量每抽子液体装满计量罐的时间,并设置该时间段为一个叠加周期T;
S3:对第一数组中的一个叠加周期T内的液面高度进行降噪处理:取数组arr1中任意一点(№1i,t1i,h1i),以一个叠加周期T为步长,分别对一个叠加周期T内各时间点的液面的高度数据进行叠加去噪,形成第二数组arr2,记为(№21,t21,h21)、(№22,t22,h22)……(№2i,t2i,h2i)……(№2n,t2n,h2n);
S4:计算液面上升斜率ai:以一个叠加周期T为步长,利用最小二乘法,求一个叠加周期T内各时间点的液面上升斜率ai,根据该一个叠加周期T内各时间点的液面上升斜率ai的计算结果建立第三数组arr3,记为(№31,t31,a1)、(№32,t32,a2)……(№3i,t3i,ai)……(№3n,t3n,an);
S5:根据第三数组arr3,判断计量罐的进液时段;
S6:测量一个叠加周期T内各时间点的实时瞬时含水率P,求解该一个叠加周期T内的平均综合含水率、液面高度差、每抽子的出液量、每抽子的出水量和每抽子的出油量。
本发明实现自动计算每抽子的产量、综合含水,通过判断抽汲作业每抽子的开始时间、结束时间、每抽汲一次液面准确的上升高度数据以及准确的实时含水率数据。解决了现有的计量检测不准确的问题,有效提高了检测的准确率。
可根据实际需要,对上述试油抽汲作业的自动计量的方法作进一步优化或/和改进:
如附图1所示,上述在S3中,因降噪后的时间t2i与降噪前的时间t1i是指同一个时间点,根据t2i=t1i,计算一个叠加周期T内任一时间点降噪后的液面高度h2i,公式如下:
其中,n为t1i-T/8至t1i+T/8之间对应的时间点的个数。
如附图1所示,S4中,因时间t3i与时间t2i是指同一个时间点,根据t3i=t2i,计算一个叠加周期T内任一时间点的液面上升斜率ai,公式如下:
如附图1所示,S5中,比较数组arr3的各项,根据公式(2)计算ai-ai-1,判断计量罐的进液时段;
S51:若ai-ai-1>0,则表示液面高度上涨,流体进入计量罐;
S52:若ai-ai-1≤0,则表示液面高度未上涨,无流体进入计量罐。
如附图1所示,S6中,当ai-ai-1>0时,标记该时间点i对应的液面上升的起始时间为Tksi,标记该时间点i对应的液面高度为Hksi;当ai-ai-1≤0时,该时间点i对应的液面上升的结束点时间Tjsi,标记该时间点i的液面高度Hjsi;
S61:计算Tksi至Tjsi时间段内的平均综合含水率,公式如下:
其中n为Tksi至Tjsi时间段内的时间点的个数,P为Tksi至Tjsi时间段内的各时间点的实时瞬时含水率;
S62:计算Tksi至Tjsi时间段的液面高度差,计算公式如下:
Hgc=Hjsi-Hksi (4)
S63:计算每抽子的出液量,计算公式如下:
V液=Hgc*Vjlg (5)
其中,Vjlg为计量罐的体积系数,根据实际使用的计量罐的大小进行设置;
S64:每抽子的出油量,计算公式如下:
实施例二:如图2、3、表1、2、3所示,根据表1、表2采集的某XX井的液面高度和实时瞬时含水率,使用上述试油抽汲作业的自动计量的方法,自动计算抽汲作业中每抽汲的出液量、出油量、综合含水等参数。
如附图2、表1所示,在平面十字坐标系中,以t为横坐标,以液面高度h为纵坐标,绘制关于时间-液面高度的曲线。经数据对比可发现,并经多次现场试验证明:当流体进入计量罐时,液面高度上升至装满的时间为2min左右,将该时间段记为一个叠加周期T(常数)。
如附图3、表2所示,在平面十字坐标系中,以t为横坐标,以含水率为纵坐标,绘制关于时间-含水率的曲线;含水率的计算过程是:计量出产出的液体总体积,记为V总;通过把水和油分离,并分别计量油记为Ⅴ油,水体积记为,V水,含水率的计算公式:
含水率=(Ⅴ水/V总)×100%
若含水率为0且同时液面上升,代表抽出液为纯油;含水率不为0且同时液面上升,判断为油水同出;若液面不上升,不论含水率是否为0,都表示没抽汲。
如表3所示,该一个叠加周期内采用实施例一的方法自动计算出的试油抽汲过程的液面高度差、抽出液量、抽出油量、抽出水量、抽前液面高度和抽后的液面高度。
以上技术特征构成了本发明的实施例,其具有较强的适应性和实施效果,可根据实际需要增减非必要的技术特征,来满足不同情况的需求。
表1某xx井在设定时间内现场采集的数据表一
表2某xx井在设定时间内现场采集的数据表二
表3某xx井在设定时间内现场采集的计算结果示意表
Claims (1)
1.一种试油抽汲作业的自动计量的方法,其特征在于,包括以下步骤:
S1:设置自动计量的时间段,采集该时间段内各时间点的计量罐的液面高度,建立基于时间-液面高度的第一数组arr1,记为(№11,t11,h11)、(№12,t12,h12)……(№1i,t1i,h1i)……(№1n,t1n,h1n),其中,i取值为(1≤i≤n),№为序号,t为时间,h为液面高度;
S2:根据各时间点对应的计量罐的液面高度变化情况,测量每抽子液体装满计量罐的时间,并设置该时间段为一个叠加周期T;
S3:对第一数组中的一个叠加周期T内的液面高度进行降噪处理:取数组arr1中任意一点(№1i,t1i,h1i),以一个叠加周期T为步长,分别对一个叠加周期T内各时间点的液面的高度数据进行叠加去噪,形成第二数组arr2,记为(№21,t21,h21)、(№22,t22,h22)……(№2i,t2i,h2i)……(№2n,t2n,h2n),具体包括:
因降噪后的时间点t2i与降噪前的时间点t1i是指同一个时间点,则根据t2i=t1i,该一个叠加周期T内任一时间点降噪后的液面高度h2i计算公式如下:
其中,n为t1i-T/8至t1i+T/8之间对应的时间点的个数;
S4:计算液面上升斜率ai:以一个叠加周期T为步长,利用最小二乘法,求一个叠加周期T内各时间点的液面上升斜率ai,根据该一个叠加周期T内各时间点的液面上升斜率ai的计算结果建立第三数组arr3,记为(№31,t31,a1)、(№32,t32,a2)……(№3i,t3i,ai)……(№3n,t3n,an),具体包括:
因时间t3i与时间t2i是指同一个时间点,则根据t3i=t2i,计算一个叠加周期T内任一时间点的液面上涨斜率ai,公式如下:
S5:根据第三数组arr3,根据公式(2)计算ai-ai-1,判断计量罐的进液时段;
S51:若ai-ai-1>0,则表示液面高度上涨,流体进入计量罐;
S52:若ai-ai-1≤0,则表示液面高度未上涨,无流体进入计量罐;
S6:测量一个叠加周期T内各时间点的实时瞬时含水率P,求解该一个叠加周期T内的平均综合含水率P、液面高度差Hgc、每抽子的出液量V液、每抽子的出水量V水和每抽子的出油量V油,具体包括:
当ai-ai-1>0时,标记该时间点i对应的液面上升的起始时间为Tksi,标记该时间点i对应的液面高度为Hksi;当ai-ai-1≤0时,该时间点i对应的液面上升的结束点时间Tjsi,标记该时间点i的液面高度Hjsi;
其中,n为Tksi至Tjsi时间段内的时间点的个数,P为Tksi至Tjsi时间段内的各时间点的实时瞬时含水率;
S62:计算Tksi至Tjsi时间段的液面高度差,计算公式如下:
Hgc=Hjsi-Hksi (4)
S63:计算每抽子的出液量,计算公式如下:
V液=Hgc*Vjlg (5)
其中,Vjlg为计量罐的体积系数,根据实际使用的计量罐的大小进行设置;
S64:每抽子的出油量,计算公式如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910345688.XA CN110219639B (zh) | 2019-04-26 | 2019-04-26 | 试油抽汲作业的自动计量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910345688.XA CN110219639B (zh) | 2019-04-26 | 2019-04-26 | 试油抽汲作业的自动计量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110219639A CN110219639A (zh) | 2019-09-10 |
CN110219639B true CN110219639B (zh) | 2023-02-03 |
Family
ID=67820069
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910345688.XA Active CN110219639B (zh) | 2019-04-26 | 2019-04-26 | 试油抽汲作业的自动计量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110219639B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113294139B (zh) * | 2020-08-05 | 2023-09-26 | 中国石油天然气股份有限公司 | 低渗油藏捞油井捞油深度的确定方法及装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA1277229C (en) * | 1986-09-16 | 1990-12-04 | Roy A. Bobo | Gas lift system |
CN2651718Y (zh) * | 2003-10-30 | 2004-10-27 | 浙江大学 | 油井油气水三相流量自动计量装置 |
CN2682203Y (zh) * | 2003-12-29 | 2005-03-02 | 吉林石油集团有限责任公司 | 试油抽吸深度自动计量装置 |
CN201292836Y (zh) * | 2008-09-17 | 2009-08-19 | 中国石油化工股份有限公司胜利油田分公司采油工艺研究院 | 探井试油连续计量装置 |
CN202431263U (zh) * | 2012-01-15 | 2012-09-12 | 洛阳乾禾仪器有限公司 | 一种抽油井自动连续计量装置 |
CN103615216A (zh) * | 2013-11-01 | 2014-03-05 | 延长油田股份有限公司 | 低产井数控油管捞油装置及捞油方法 |
CN108798633A (zh) * | 2018-07-04 | 2018-11-13 | 中国石油集团西部钻探工程有限公司 | 试油单井监测计量方法和系统 |
CN208777977U (zh) * | 2018-07-04 | 2019-04-23 | 中国石油集团西部钻探工程有限公司 | 试油单井原油计量装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN200943482Y (zh) * | 2006-01-09 | 2007-09-05 | 西安凯诺石油科技有限公司 | 抽汲作业测控装置 |
US9394783B2 (en) * | 2011-08-26 | 2016-07-19 | Schlumberger Technology Corporation | Methods for evaluating inflow and outflow in a subterranean wellbore |
CN107587868B (zh) * | 2017-10-16 | 2024-01-26 | 陕西航天泵阀科技集团有限公司 | 一种油井计量集成装置 |
CA2983541C (en) * | 2017-10-24 | 2019-01-22 | Exxonmobil Upstream Research Company | Systems and methods for dynamic liquid level monitoring and control |
CN109098702A (zh) * | 2018-07-04 | 2018-12-28 | 中国石油集团西部钻探工程有限公司 | 试油单井实时监测平台 |
-
2019
- 2019-04-26 CN CN201910345688.XA patent/CN110219639B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA1277229C (en) * | 1986-09-16 | 1990-12-04 | Roy A. Bobo | Gas lift system |
CN2651718Y (zh) * | 2003-10-30 | 2004-10-27 | 浙江大学 | 油井油气水三相流量自动计量装置 |
CN2682203Y (zh) * | 2003-12-29 | 2005-03-02 | 吉林石油集团有限责任公司 | 试油抽吸深度自动计量装置 |
CN201292836Y (zh) * | 2008-09-17 | 2009-08-19 | 中国石油化工股份有限公司胜利油田分公司采油工艺研究院 | 探井试油连续计量装置 |
CN202431263U (zh) * | 2012-01-15 | 2012-09-12 | 洛阳乾禾仪器有限公司 | 一种抽油井自动连续计量装置 |
CN103615216A (zh) * | 2013-11-01 | 2014-03-05 | 延长油田股份有限公司 | 低产井数控油管捞油装置及捞油方法 |
CN108798633A (zh) * | 2018-07-04 | 2018-11-13 | 中国石油集团西部钻探工程有限公司 | 试油单井监测计量方法和系统 |
CN208777977U (zh) * | 2018-07-04 | 2019-04-23 | 中国石油集团西部钻探工程有限公司 | 试油单井原油计量装置 |
Non-Patent Citations (6)
Title |
---|
地层测试三开抽汲若干问题探讨;李民乐等;《油气井测试》;20160831;第25卷(第4期);第36-37页 * |
定流压抽汲井产能预测方法;刘洪等;《断块油气田》;20080131(第1期);第52-54页 * |
油管泵抽汲排液工艺的研究与应用;余东合等;《石油钻采工艺》;20061231(第1期);第25-27页 * |
液位报警器在密闭计量罐上的应用;王继承等;《中国石油和化工标准与质量》;20130228;第34卷(第3期);第97页 * |
致密油水平井试油工艺及产能评价方法研究;于昊等;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20180215(第2期);第B019-78页 * |
试油试采分离计量系统的研制与应用;张东军等;《油气井测试》;20020630(第3期);第64-66页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110219639A (zh) | 2019-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109872018B (zh) | 一种基于示功图的抽油井产量计算方法 | |
CN107121369A (zh) | 一种全自动达西定律实验系统 | |
CN110219639B (zh) | 试油抽汲作业的自动计量的方法 | |
CN113338904B (zh) | 一种海上油井产能评价方法 | |
CN105956938B (zh) | 缝洞型油藏动态储量计算的方法 | |
CN108414343A (zh) | 一种测量粗粒土大型三轴试验中试样橡皮膜嵌入量的方法 | |
CN108920429A (zh) | 一种水位动态监测的异常数据分析方法 | |
CN105181934A (zh) | 基于一维等应变和应力松弛试验估算软土次固结系数方法 | |
CN104405374A (zh) | 一种致密气藏储层应力敏感性的测量方法 | |
CN109184661B (zh) | 用于底水油藏水平井高产液位置识别的监测方法及其系统 | |
CN108593460B (zh) | 基于直剪试验确定土体抗剪强度的动态精确计算方法 | |
CN115856092B (zh) | 基于声发射数据和应力数据确定岩石裂纹起裂应力的方法 | |
CN110909466B (zh) | 提高差孔隙分选储层平均毛管压力曲线计算精度的方法 | |
CN110348040B (zh) | 一种基于均方收敛条件的不均匀沉降函数拟合方法 | |
CN115977586B (zh) | 一种海上气井产能评价新方法 | |
CN106018761A (zh) | 一种建筑物混凝土品质的无损检测系统和方法 | |
CN113153277B (zh) | 一种钻井液溢流或漏失的预警方法 | |
CN114664387A (zh) | 基于累积效应的油井化学堵水性能室内评价方法 | |
CN112196519B (zh) | 一种油井非稳态连续动液面的检测方法 | |
CN111241651A (zh) | 新型水驱图版的制作方法 | |
CN110206530B (zh) | 试油作业中数据处理计量方法及系统 | |
CN117520996B (zh) | 一种用于涵管的实时压力测试方法 | |
CN109472083B (zh) | 基于统计过程控制的土木工程结构损伤检测方法 | |
CN114722328B (zh) | 一种水库常年回水区淤积泥沙粒径的计算方法 | |
CN114528700B (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 |