CN110826249A - 一种铀矿综合测井曲线井眼环境校正方法 - Google Patents

一种铀矿综合测井曲线井眼环境校正方法 Download PDF

Info

Publication number
CN110826249A
CN110826249A CN201911128857.0A CN201911128857A CN110826249A CN 110826249 A CN110826249 A CN 110826249A CN 201911128857 A CN201911128857 A CN 201911128857A CN 110826249 A CN110826249 A CN 110826249A
Authority
CN
China
Prior art keywords
correction
curve
logging
shaft
formula
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
CN201911128857.0A
Other languages
English (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.)
Beijing Research Institute of Uranium Geology
Original Assignee
Beijing Research Institute of Uranium Geology
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 Beijing Research Institute of Uranium Geology filed Critical Beijing Research Institute of Uranium Geology
Priority to CN201911128857.0A priority Critical patent/CN110826249A/zh
Publication of CN110826249A publication Critical patent/CN110826249A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK 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
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK 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/005—Testing the nature of borehole walls or the formation by using drilling mud or cutting data

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明属于砂岩铀矿地球物理测井技术领域，具体涉及一种铀矿综合测井曲线井眼环境校正方法。本发明包括如下步骤：步骤1建立井筒模型；步骤2建立自然伽马测井曲线的校正图版和校正公式；步骤3建立补偿密度测井曲线的校正图版和校正公式；步骤4建立声波时差测井曲线的校正图版和校正公式；步骤5建立三侧向电阻率测井曲线的校正图版和校正公式；步骤6编写砂岩铀矿综合测井曲线可视化环境校正软件。本发明能够满足砂岩型铀矿地球物理测井资料精细化处理与解释的需求。

Description

一种铀矿综合测井曲线井眼环境校正方法
技术领域
本发明属于砂岩铀矿地球物理测井技术领域，具体涉及一种铀矿综合测井曲线井眼环境校正方法。
背景技术
野外采集地球物理测井数据时，钻孔内的环境如井径扩大、泥浆密度与矿化度变化、泥饼和泥浆侵入带厚度、地层温度与压力，以及测井探管外径等非地层因素会对地球物理测井曲线产生重要的影响。
国内外石油和煤田等测井公司均制作了与他们所生产仪器相配套的测井曲线环境校正解释图版，也研发了相配套的环境校正软件。而国内地矿相关单位却没有开展针对小口径钻孔测井曲线的环境校正工作。并且测井曲线的环境影响校正图版及校正软件是针对各自测井系统的特征在不同环境条件下实验模拟的，具有一定的局限性和适用范围。
基于以上问题，有必要根据HD4002系统的特征，建立一套适合该系统铀矿综合测井曲线的井眼环境校正方法，并研发出相应的处理软件，以此消除测井过程中的环境影响，得到反映原状地层的真实测井曲线。
发明内容
本发明的目的是建立一套适合HD4002系统铀矿综合测井曲线的井眼环境校正方法，满足砂岩型铀矿地球物理测井资料精细化处理与解释的需求。
本发明采用的技术方案：
一种铀矿综合测井曲线井眼环境校正方法，包括如下步骤：步骤1建立井筒模型；步骤2建立自然伽马测井曲线的校正图版和校正公式；步骤3建立补偿密度测井曲线的校正图版和校正公式；步骤4建立声波时差测井曲线的校正图版和校正公式；步骤5建立三侧向电阻率测井曲线的校正图版和校正公式；步骤6编写砂岩铀矿综合测井曲线可视化环境校正软件。
所述步骤1中，井筒模型包括三大类：一类是井筒直径大于测井探头直径1mm；二类是井筒直径均匀变化，并在井筒内放置水；三类是井筒直径均匀变化，并在井筒内放置不同密度和矿化度的泥浆液体。
所述步骤2的具体步骤为，将采集自然伽马测井曲线的探头放置在三种类型的井筒模型中，分别采集自然伽马测井数据，剔除异常值，取平均值，建立自然伽马校正图版，并根据自然伽马校正图版中曲线的拟合得到的校正公式。
所述校正公式为
GRc＝A·GR·exp[0.0358·ρm·(d-di)-0.6021] (1)
式中，GR和GRc分别表示环境校正前和校正后的伽马测井值，单位CPS；d和di分别表示井筒直径和仪器外径，单位为mm；ρm为钻孔内泥浆密度，单位为g/cm3；A是复合参数，它与仪器的外径、扩径率和地区地质背景有关。
所述步骤3的具体步骤为，将采集补偿密度测井曲线的探头放置在三种类型的井筒模型中，分别采集补偿密度测井数据，剔除异常值，取平均值，建立补偿密度校正图版，并根据补偿密度校正图版中曲线的拟合得到的校正公式。
所述校正公式为
DENc＝DEN+A·(d-BIT)·(DENt-ρm) (2)
式中，DEN和DENc分别表示环境校正前和校正后的密度测井值，单位为g/cm3；BIT表示钻头直径，单位mm；d为井筒直径，单位mm；DENt为纯砂岩地层密度，单位为g/cm3；ρm为钻孔内泥浆密度，单位为g/cm3；A为受地区限制的校正比例系数。
所述步骤4的具体步骤为，将采集声波时差测井曲线的探头放置在一类和二类井筒模型中，采集声波时差测井数据，剔除异常值，取平均值，建立声波时差校正图版，并根据声波时差校正图版中曲线的拟合得到的校正公式。
所述校正公式为
ACc＝AC-A·((d-BIT)/12)·BIT (3)
式中，AC和ACc分别表示环境校正前和校正后的密度测井值，单位为μs/m；BIT表示钻头直径，单位mm；d为井筒直径，单位mm；A为受地区限制的校正比例系数。
所述步骤5的具体步骤为，将采集三侧向电阻率测井曲线的探头放置在三种类型的井筒模型中，采集三侧向电阻率测井数据，剔除异常值，取平均值，建立三侧向电阻率校正图版，并根据三侧向电阻率校正图版中曲线的拟合得到的校正公式。
所述校正公式为
Rt_corre/Rt＝1.22-0.039·lg(Rt/Rm+1.983) (4)
式中，Rt和Rt_corre分别为校正前、后的三侧向电阻率和浅侧向电阻率，Rm为泥浆电阻率。
本发明的有益效果：
砂泥岩地层井筒内扩径和泥浆对测井资料解译地层的影响是巨大的，需要进行环境校正。研究测井资料的环境校正技术，求取能够反应真实地层电阻率、密度、声波等信息，校正后的测井曲线能显著提高测井资料解译的精度和合理性，对其它地区开展测井环境校正具有一定的参考和借鉴意义。
附图说明
图1是一类井筒模型；
图2为二类井筒模型；
图3为三类井筒模型；
图4是HD4002系统砂岩铀矿综合测井曲线井眼环境影响校正流程图；
图5是HD4002系统自然伽马测井曲线井径和泥浆校正图版；
图6是HD4002系统补偿密度测井曲线井径和泥浆校正图版；
图7是HD4002系统声波测井曲线井径校正图版；
图8是HD4002系统三侧向电阻率测井曲线井径和泥浆校正图版。
具体实施方式
下面结合附图和具体实施例对本发明提供的一种铀矿综合测井曲线井眼环境校正方法作进一步详细描述。
如图4所示，本发明提供的一种铀矿综合测井曲线井眼环境校正方法，包括如下步骤：
步骤(1)建立如图1-图3所示的井筒模型，其中一类井筒模型的直径设置有41mm，46mm和61mm(采集自然伽马测井曲线的探头直径为40mm，采集声波时差测井曲线的探头直径为45mm，采集三侧向电阻率和补偿密度测井曲线的探头直径为60mm)。
二类井筒模型的直径设置有80、120、160、200、240、280和320mm。
三类井筒模型的直径设置有80、120、160、200、240、280和320mm，泥浆密度设置有1.0、1.2、1.4、1.6、1.8和2.0g/cm3，泥浆电阻率设置有1.0、1.2、1.4、1.6、1.8和2.0Ω·m；
步骤(2)建立自然伽马测井曲线的校正图版和校正公式。将采集自然伽马测井曲线的探头放置在三种类型的井筒模型中，分别采集不受井径扩大和泥浆影响的自然伽马测井曲线、仅受井径扩大影响的自然伽马测井曲线、同时受井径扩大和泥浆影响的自然伽马测井曲线。
在每种类型井筒中采集多组自然伽马测井数据，剔除异常值，取平均值，建立自然伽马校正图版(图5所示)，横坐标设置为井筒直径和仪器直径的差值，纵坐标为伽马校正值和测量值的比值。
通过对图3校正图版中曲线的拟合得到的校正公式如下：
GRc＝A·GR·exp[0.0358·ρm·(d-di)-0.6021](1)
式中，GR和GRc分别表示环境校正前和校正后的伽马测井值，单位CPS；d和di分别表示井筒直径和仪器外径，单位为mm；ρm为钻孔内泥浆密度，单位为g/cm3；A是复合参数，它与仪器的外径、扩径率和地区地质背景有关。
步骤(3)建立补偿密度测井曲线的校正图版和校正公式。将采集补偿密度测井曲线的探头放置在三种类型的井筒模型中，设置7种砂岩地层的密度，分别为1.8、2.0、2.2、2.4、2.6、2.8和3.0g/cm3。分别采集不受井径扩大和泥浆影响的补偿密度测井曲线、仅受井径扩大影响的补偿密度测井曲线、同时受井径扩大和泥浆影响的补偿密度测井曲线。
在每种类型井筒中采集多组补偿密度测井数据，剔除异常值，取平均值，建立补偿密度校正图版(图6所示)，横坐标为钻孔直径，纵坐标为密度曲线的校正量。
通过对图4校正图版中曲线的拟合得到的校正公式如下：
DENc＝DEN+A·(d-BIT)·(DENt-ρm) (2)
式中，DEN和DENc分别表示环境校正前和校正后的密度测井值，单位为g/cm3；BIT表示钻头直径，单位mm；d为井筒直径，单位mm；DENt为纯砂岩地层密度，单位为g/cm3；ρm为钻孔内泥浆密度，单位为g/cm3；A为受地区限制的校正比例系数，在0.008左右进行调整。
步骤(4)建立声波时差测井曲线的校正图版和校正公式。HD4002系统补偿声波测井探头对井眼影响有较强的补偿作用。因此，声波时差测井相对于密度测井，其受井眼影响较小，但当扩径严重(扩径率大于30％)或者井壁很不规则时，声波时差明显增大。受泥浆密度影响不大。
将采集声波时差测井曲线的探头放置在一类和二类井筒模型中，采集不受井径扩大影响和受井径扩大影响的声波时差测井曲线。在每种类型井筒中采集多组声波时差测井数据，剔除异常值，取平均值，建立声波时差校正图版(图7所示)，横坐标为钻孔直径，纵坐标为声波时差曲线的校正量。
通过对图5校正图版中曲线的拟合得到的校正公式如下：
ACc＝AC-A·((d-BIT)/12)·BIT (3)
式中，AC和ACc分别表示环境校正前和校正后的密度测井值，单位为μs/m；BIT表示钻头直径，单位mm；d为井筒直径，单位mm；A为受地区限制的校正比例系数，在0.05左右进行调整。
步骤(5)建立三侧向电阻率测井曲线的校正图版和校正公式。将采集三侧向电阻率测井曲线的探头放置在三种类型的井筒模型中，采集不受井径扩大和泥浆影响的三侧向电阻率测井曲线、仅受井径扩大影响的三侧向电阻率测井曲线、同时受井径扩大和泥浆影响的三侧向电阻率测井曲线。
在每种类型井筒中采集多组三侧向电阻率测井数据，剔除异常值，取平均值，建立三侧向电阻率校正图版(图8所示)，横坐标为地层与泥浆电阻率的比值，纵坐标为三侧向电阻率的校正量。
通过对图6校正图版中曲线的拟合得到的校正公式如下：
Rt_corre/Rt＝1.22-0.039·lg(Rt/Rm+1.983) (4)
式中，Rt和Rt_corre分别为校正前、后的三侧向电阻率和浅侧向电阻率，Rm为泥浆电阻率。
步骤(6)编写HD4002系统砂岩铀矿综合测井曲线可视化环境校正软件。将预处理后的测井曲线输入到环境校正软件中，实现测井曲线的自动环境校正。

Claims (10)

1.一种铀矿综合测井曲线井眼环境校正方法，其特征在于：包括如下步骤：步骤(1)建立井筒模型；步骤(2)建立自然伽马测井曲线的校正图版和校正公式；步骤(3)建立补偿密度测井曲线的校正图版和校正公式；步骤(4)建立声波时差测井曲线的校正图版和校正公式；步骤(5)建立三侧向电阻率测井曲线的校正图版和校正公式；步骤(6)编写砂岩铀矿综合测井曲线可视化环境校正软件。
2.根据权利要求1所述一种铀矿综合测井曲线井眼环境校正方法，其特征在于：所述步骤(1)中，井筒模型包括三大类：一类是井筒直径大于测井探头直径1mm；二类是井筒直径均匀变化，并在井筒内放置水；三类是井筒直径均匀变化，并在井筒内放置不同密度和矿化度的泥浆液体。
3.根据权利要求2所述一种铀矿综合测井曲线井眼环境校正方法，其特征在于：所述步骤(2)的具体步骤为，将采集自然伽马测井曲线的探头放置在三种类型的井筒模型中，分别采集自然伽马测井数据，剔除异常值，取平均值，建立自然伽马校正图版，并根据自然伽马校正图版中曲线的拟合得到的校正公式。
4.根据权利要求3所述一种铀矿综合测井曲线井眼环境校正方法，其特征在于：所述校正公式为
GRc＝A·GR·exp[0.0358·ρm·(d-di)-0.6021] (1)
式中，GR和GRc分别表示环境校正前和校正后的伽马测井值，单位CPS；d和di分别表示井筒直径和仪器外径，单位为mm；ρm为钻孔内泥浆密度，单位为g/cm3；A是复合参数，它与仪器的外径、扩径率和地区地质背景有关。
5.根据权利要求3所述一种铀矿综合测井曲线井眼环境校正方法，其特征在于：所述步骤(3)的具体步骤为，将采集补偿密度测井曲线的探头放置在三种类型的井筒模型中，分别采集补偿密度测井数据，剔除异常值，取平均值，建立补偿密度校正图版，并根据补偿密度校正图版中曲线的拟合得到的校正公式。
6.根据权利要求5所述一种铀矿综合测井曲线井眼环境校正方法，其特征在于：所述校正公式为
DENc＝DEN+A·(d-BIT)·(DENt-ρm) (2)
式中，DEN和DENc分别表示环境校正前和校正后的密度测井值，单位为g/cm3；BIT表示钻头直径，单位mm；d为井筒直径，单位mm；DENt为纯砂岩地层密度，单位为g/cm3；ρm为钻孔内泥浆密度，单位为g/cm3；A为受地区限制的校正比例系数。
7.根据权利要求5所述一种铀矿综合测井曲线井眼环境校正方法，其特征在于：所述步骤(4)的具体步骤为，将采集声波时差测井曲线的探头放置在一类和二类井筒模型中，采集声波时差测井数据，剔除异常值，取平均值，建立声波时差校正图版，并根据声波时差校正图版中曲线的拟合得到的校正公式。
8.根据权利要求7所述一种铀矿综合测井曲线井眼环境校正方法，其特征在于：所述校正公式为
ACc＝AC-A·((d-BIT)/12)·BIT (3)
式中，AC和ACc分别表示环境校正前和校正后的密度测井值，单位为μs/m；BIT表示钻头直径，单位mm；d为井筒直径，单位mm；A为受地区限制的校正比例系数。
9.根据权利要求7所述一种铀矿综合测井曲线井眼环境校正方法，其特征在于：所述步骤(5)的具体步骤为，将采集三侧向电阻率测井曲线的探头放置在三种类型的井筒模型中，采集三侧向电阻率测井数据，剔除异常值，取平均值，建立三侧向电阻率校正图版，并根据三侧向电阻率校正图版中曲线的拟合得到的校正公式。
10.根据权利要求9所述一种铀矿综合测井曲线井眼环境校正方法，其特征在于：所述校正公式为
Rt_corre/Rt＝1.22-0.039·lg(Rt/Rm+1.983) (4)
式中，Rt和Rt_corre分别为校正前、后的三侧向电阻率和浅侧向电阻率，Rm为泥浆电阻率。
CN201911128857.0A 2019-11-18 2019-11-18 一种铀矿综合测井曲线井眼环境校正方法 Pending CN110826249A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911128857.0A CN110826249A (zh) 2019-11-18 2019-11-18 一种铀矿综合测井曲线井眼环境校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911128857.0A CN110826249A (zh) 2019-11-18 2019-11-18 一种铀矿综合测井曲线井眼环境校正方法

Publications (1)

Publication Number Publication Date
CN110826249A true CN110826249A (zh) 2020-02-21

Family

ID=69556393

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911128857.0A Pending CN110826249A (zh) 2019-11-18 2019-11-18 一种铀矿综合测井曲线井眼环境校正方法

Country Status (1)

Country Link
CN (1) CN110826249A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112228049A (zh) * 2020-09-30 2021-01-15 核工业二〇八大队 地浸工艺钻孔的测量方法
CN113356834A (zh) * 2021-07-14 2021-09-07 中国石油天然气股份有限公司 井眼声波时差的校正方法、校正装置和电子设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5021653A (en) * 1990-02-07 1991-06-04 Schlumberger Technology Corporation Geochemical logging apparatus and method for determining concentrations of formation elements next to a borehole
US20030024305A1 (en) * 2001-07-31 2003-02-06 Ireland Peter J. Method and apparatus for measurement alignment
CN105114063A (zh) * 2015-08-04 2015-12-02 中国石油大学(华东) 一种斜井双侧向测井曲线校正方法
CN107092031A (zh) * 2017-06-01 2017-08-25 中国石油天然气集团公司 一种声波测井数据校正方法及装置
CN109407173A (zh) * 2018-09-29 2019-03-01 核工业北京地质研究院 基于常规测井曲线的岩性精细化和自动化识别方法
CN110145301A (zh) * 2019-04-28 2019-08-20 中国石油天然气集团有限公司 一种随钻密度成像测井仪环境校正方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5021653A (en) * 1990-02-07 1991-06-04 Schlumberger Technology Corporation Geochemical logging apparatus and method for determining concentrations of formation elements next to a borehole
US20030024305A1 (en) * 2001-07-31 2003-02-06 Ireland Peter J. Method and apparatus for measurement alignment
CN105114063A (zh) * 2015-08-04 2015-12-02 中国石油大学(华东) 一种斜井双侧向测井曲线校正方法
CN107092031A (zh) * 2017-06-01 2017-08-25 中国石油天然气集团公司 一种声波测井数据校正方法及装置
CN109407173A (zh) * 2018-09-29 2019-03-01 核工业北京地质研究院 基于常规测井曲线的岩性精细化和自动化识别方法
CN110145301A (zh) * 2019-04-28 2019-08-20 中国石油天然气集团有限公司 一种随钻密度成像测井仪环境校正方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
田志杰 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112228049A (zh) * 2020-09-30 2021-01-15 核工业二〇八大队 地浸工艺钻孔的测量方法
CN113356834A (zh) * 2021-07-14 2021-09-07 中国石油天然气股份有限公司 井眼声波时差的校正方法、校正装置和电子设备
CN113356834B (zh) * 2021-07-14 2023-04-18 中国石油天然气股份有限公司 井眼声波时差的校正方法、校正装置和电子设备

Similar Documents

Publication Publication Date Title
CN106761677B (zh) 页岩气水平井单井产能的测井预测方法
CN107783187B (zh) 一种将测井速度和地震速度结合建立三维速度场的方法
CN103857876A (zh) 用于执行井眼压裂作业的系统和方法
KR101415199B1 (ko) 전기비저항 및 밀도 검층 데이터를 이용한 셰일가스전 스윗스팟 수평시추구간의 슬로우니스값, 영률, 포아송비 및 취성 추정방법
CN107907910B (zh) 一种不同岩性油藏横波测井确定方法
KR101415196B1 (ko) 전기비저항 및 밀도 검층 데이터를 이용한 셰일가스전 스윗스팟 수평시추구간의 파쇄구간 선정방법
CN110826249A (zh) 一种铀矿综合测井曲线井眼环境校正方法
CN106285642B (zh) 一种基于地震资料的地层孔隙压力预测方法
Zhao et al. Key seismic survey technologies for deep complex geological structures: A case study of the northern section of the Longmenshan Fault Fold Belt in the Sichuan Basin
CN116122801A (zh) 一种页岩油水平井体积压裂可压性综合评价方法
KR20180014294A (ko) 물리검층 자료를 이용한 s파 속도 산출방법
CN108678672A (zh) 深水浅部复杂岩性地层导管高效置入及井口稳定预测方法
Gieske et al. An analysis of earth-tide-induced groundwater flow in eastern Botswana
AU2021105676A4 (en) Near-source electromagnetic seismoelectric combined GR inversion method
CN113006779B (zh) 一种砾岩地层岩石可钻性级值的评价方法和校正方法
Ziqing Geosteering methods of a dynamic geological model based on Kalman filter
CN105589098B (zh) 碳酸盐岩去除沉积泥质横向影响的储层反演方法和系统
Batir et al. Preliminary model of fracture and stress state in the Hellisheidi Geothermal Field, Hengill Volcanic System, Iceland
CN111768104B (zh) 单斜砂泥岩互层型坝基岩性变化程度定量评价方法
Li et al. Evaluation of lithology variations in layered red beds with depth: An example of the Yellow River Guxian Dam, NW China
CN109339776B (zh) 一种各向异性地层地应力方位测量方法
CN106990448A (zh) 一种海相碎屑岩水平井隔夹层测井解释方法及装置
CN112483172A (zh) 一种人为再造阻水边界的方法
Zhou Study on reasonable bottom hole flow pressure limits in different development stages of tight reservoirs in Block X
CN118113980A (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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200221