CN106646602A - 一种基于多种震源模型的震后快速生成震动图的方法 - Google Patents
一种基于多种震源模型的震后快速生成震动图的方法 Download PDFInfo
- Publication number
- CN106646602A CN106646602A CN201611256348.2A CN201611256348A CN106646602A CN 106646602 A CN106646602 A CN 106646602A CN 201611256348 A CN201611256348 A CN 201611256348A CN 106646602 A CN106646602 A CN 106646602A
- Authority
- CN
- China
- Prior art keywords
- earthquake
- seismic
- region
- analysis
- model
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000004458 analytical method Methods 0.000 claims abstract description 21
- 238000001228 spectrum Methods 0.000 claims description 9
- 230000001133 acceleration Effects 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 238000012544 monitoring process Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000012732 spatial analysis Methods 0.000 claims description 3
- 240000007594 Oryza sativa Species 0.000 claims description 2
- 235000007164 Oryza sativa Nutrition 0.000 claims description 2
- 235000009566 rice Nutrition 0.000 claims description 2
- 238000007405 data analysis Methods 0.000 abstract description 3
- 230000010355 oscillation Effects 0.000 abstract 4
- 239000011435 rock Substances 0.000 abstract 1
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 241001416181 Axis axis Species 0.000 description 1
- 235000011464 Pachycereus pringlei Nutrition 0.000 description 1
- 240000006939 Pachycereus weberi Species 0.000 description 1
- 235000011466 Pachycereus weberi Nutrition 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/614—Synthetically generated data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
- G01V2210/665—Subsurface modeling using geostatistical modeling
- G01V2210/6652—Kriging
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于多种震源模型的震后快速生成震动图的方法,涉及电通信技术领域。所述方法包括:从测震台获取震区的地震信息,检索并获得震区的地质构造,将地震信息与地质构造叠加分析后,获得震区的震源模型;根据震源模型智能匹配地震动衰减关系,根据地震动衰减关系,计算获得虚拟测震台站点的工作区范围内的基岩层的地震动参数;将虚拟测震台站点插值到相邻站点,间距为500‑1000米;对插值后得到的栅格图进行配色处理,得到最终的地震动烈度图。本发明提高对地震发生后的数据分析和处理效率,提高震后预测准确度。
Description
技术领域
本发明涉及电通信技术领域,尤其涉及一种基于多种震源模型的震后快速生成震动图的方法。
背景技术
现有震后生成震动图是基于烈度衰减方程以及中国东西向点椭圆衰减公式,绘制烈度等值线实现,为震后评估提供数据支撑。
但是现有方法存在以下不足:震源模型技术单一,因仅考虑点椭圆模型,对发生大震时需要的线源以及面源模型支持较弱,同时,现有方法也未考虑实际地质构造信息对烈度线的影响,无法比对不同的衰减关系对模型结果的影响程度。
发明内容
本发明的目的在于提供一种基于多种震源模型的震后快速生成震动图的方法,从而解决现有技术中存在的前述问题。
为了实现上述目的,本发明所述基于多种震源模型的震后快速生成震动图的方法,其特征在于,所述方法包括:
S1,从测震台获取震区的地震信息,检索并获得震区的地质构造,将地震信息与地质构造叠加分析后,获得震区的震源模型;
S2,根据震源模型智能匹配地震动衰减关系,根据地震动衰减关系,计算获得虚拟测震台站点的工作区范围内的基岩层的地震动参数;
S3,将虚拟测震台站点插值到相邻站点,间距为500-1000米;
S4,对插值后得到的栅格图进行配色处理,得到最终的地震动烈度图。
优选地,步骤S1中,从测震台获取震区的地震信息,具体按照下述步骤实现:实时监控国家数字测震台网的数据备份中心服务器,如果监测到新增地震数据,则直接获取并分析新增的地震数据,得到震中位置,然后根据震中位置检索并获得震区的地址构造;如果没有监测到新增地震数据,则继续监控。
更优选地,所述分析新增地震数据包括分析震中经纬度数据和分析震级数据。
优选地,步骤S1中,将地震信息与地质构造叠加分析后,获得震区的震源模型,具体按照下述步骤实现:获取震中点位和地址构造面状图层,将震中点位与地质构造面状图层进行空间叠加分析后,确认评估当前地震使用的模型,所述模型包括点源模型和点椭圆模型。
更优选地,如果使用的模型为点椭圆模型,则求震区长短轴的信息,具体为:对地震信息中的震中信息和震区的断裂带面状数据进行空间分析,若震中点位未落在面状要素里,则判定此次地震为点源模型;若震中点位落在面状要素里,则得到当前断裂带的走向,所述走向方向为点椭圆的长轴方向。
优选地,步骤S2中所述地震动衰减关系根据衰减公式获得,所述衰减公式为公式(1):
lgY=a+b*M-c*lg(R+d*exp(e*M))+sigma (1);
根据公式(1)计算获得虚拟测震台站点的工作区范围内网格点位的峰值和加速度的值,公式(1)中,Y为加速度、速度或者谱值,M为地震震级,R为震中距;系数a、b、c、常数sigma根据地震等烈度线资料回归统计求得。
优选地,步骤S3中,通过克里金或者反距离权重插值法,将网格点位插值到相邻站点间距为800米的更密的点位。
本发明的有益效果是:
本发明所述方法提高对地震发生后的数据分析和处理效率,提高震后预测准确度,在考虑烈度衰减的同时结合地质构造信息,获取震源机制和震源破裂过程详细结果,然后自行匹配适当的衰减关系,从而计算出峰值加速度等值线、峰值速度等值线,Sa(0.3s)谱值等值线、Sa(1s)谱值计等值线、Sa(3s)谱值等值线以及地震动强度预测图,快速准确估算受灾面积、极震区分布,为后继的灾害评估提供依据。
附图说明
图1是基于多种震源模型的震后快速生成震动图的方法的流程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本发明,并不用于限定本发明。
实施例
本实施例所述基于多种震源模型的震后快速生成震动图的方法,所述方法包括:
S1,从测震台获取震区的地震信息,检索并获得震区的地质构造,将地震信息与地质构造叠加分析后,获得震区的震源模型;
S2,根据震源模型智能匹配地震动衰减关系,根据地震动衰减关系,计算获得虚拟测震台站点的工作区范围内的基岩层的地震动参数;
S3,将虚拟测震台站点插值到相邻站点,通过克里金或者反距离权重插值法,将网格点位插值到相邻站点间距为800米的更密的点位;
S4,对插值后得到的栅格图进行配色处理,得到最终的地震动烈度图。
更详细的解释说明为:
(一)步骤S1中,从测震台获取震区的地震信息,具体按照下述步骤实现:实时监控国家数字测震台网的数据备份中心服务器,如果监测到新增地震数据,则直接获取并分析新增的地震数据,得到震中位置,然后根据震中位置检索并获得震区的地址构造;如果没有监测到新增地震数据,则继续监控。
所述分析新增地震数据包括分析震中经纬度数据和分析震级数据。
(二)步骤S1中,将地震信息与地质构造叠加分析后,获得震区的震源模型,具体按照下述步骤实现:获取震中点位和地址构造面状图层,将震中点位与地质构造面状图层进行空间叠加分析后,确认评估当前地震使用的模型,所述模型包括点源模型和点椭圆模型。
如果使用的模型为点椭圆模型,则求震区长短轴的信息,具体为:对地震信息中的震中信息和震区的断裂带面状数据进行空间分析,若震中点位未落在面状要素里,则判定此次地震为点源模型;若震中点位落在面状要素里,则得到当前断裂带的走向,所述走向方向为点椭圆的长轴方向。
(三)步骤S2中所述地震动衰减关系根据衰减公式获得,所述衰减公式为公式(1):
lgY=a+b*M-c*lg(R+d*exp(e*M))+sigma (1);
根据公式(1)计算获得虚拟测震台站点的工作区范围内网格点位的峰值和加速度的值,公式(1)中,Y为加速度、速度或者谱值,M为地震震级,R为震中距;系数a、b、c、常数sigma根据地震等烈度线资料回归统计求得。
在步骤S2中,对新疆、川藏、华北、中强区域分别建立衰减关系库,根据震中位置自动判定此次模型需要采用的衰减关系,衰减关系的存储如表1,包含的字段分别为区域名称area、地震动参数名称psa、峰值加速度Tpga、长短轴axis、a1、b1是震级小于6.5级的系数,a2、b2是震级大于6.5级的系数。根据获取到的系数代入到公式(1)计算出点位的峰值加速度。
表格1为对新疆、川藏、华北、中强区域分别建立衰减关系表
通过采用本发明公开的上述技术方案,得到了如下有益的效果:
本发明所述方法提高对地震发生后的数据分析和处理效率,提高震后预测准确度,在考虑烈度衰减的同时结合地质构造信息,获取震源机制和震源破裂过程详细结果,然后自行匹配适当的衰减关系,从而计算出峰值加速度等值线、峰值速度等值线,Sa(0.3s)谱值等值线、Sa(1s)谱值计等值线、Sa(3s)谱值等值线以及地震动强度预测图,快速准确估算受灾面积、极震区分布,为后继的灾害评估提供依据。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视本发明的保护范围。
Claims (7)
1.一种基于多种震源模型的震后快速生成震动图的方法,其特征在于,所述方法包括:
S1,从测震台获取震区的地震信息,检索并获得震区的地质构造,将地震信息与地质构造叠加分析后,获得震区的震源模型;
S2,根据震源模型智能匹配地震动衰减关系,根据地震动衰减关系,计算获得虚拟测震台站点的工作区范围内的基岩层的地震动参数;
S3,将虚拟测震台站点插值到相邻站点,间距为500-1000米;
S4,对插值后得到的栅格图进行配色处理,得到最终的地震动烈度图。
2.根据权利要求1所述方法,其特征在于,步骤S1中,从测震台获取震区的地震信息,具体按照下述步骤实现:实时监控国家数字测震台网的数据备份中心服务器,如果监测到新增地震数据,则直接获取并分析新增的地震数据,得到震中位置,然后根据震中位置检索并获得震区的地址构造;如果没有监测到新增地震数据,则继续监控。
3.根据权利要求2所述方法,其特征在于,所述分析新增地震数据包括分析震中经纬度数据和分析震级数据。
4.根据权利要求1所述方法,其特征在于,步骤S1中,将地震信息与地质构造叠加分析后,获得震区的震源模型,具体按照下述步骤实现:
获取震中点位和地址构造面状图层,将震中点位与地质构造面状图层进行空间叠加分析后,确认评估当前地震使用的模型,所述模型包括点源模型和点椭圆模型。
5.根据权利要求4所述方法,其特征在于,如果使用的模型为点椭圆模型,则求震区长短轴的信息,具体为:
对地震信息中的震中信息和震区的断裂带面状数据进行空间分析,若震中点位未落在面状要素里,则判定此次地震为点源模型;若震中点位落在面状要素里,则得到当前断裂带的走向,所述走向方向为点椭圆的长轴方向。
6.根据权利要求1所述方法,其特征在于,步骤S2中所述地震动衰减关系根据衰减公式获得,所述衰减公式为公式(1):
lgY=a+b*M-c*lg(R+d*exp(e*M))+sigma (1);
根据公式(1)计算获得虚拟测震台站点的工作区范围内网格点位的峰值和加速度的值,公式(1)中,Y为加速度、速度或者谱值,M为地震震级,R为震中距;系数a、b、c、常数sigma根据地震等烈度线资料回归统计求得。
7.根据权利要求1所述方法,其特征在于,步骤S3中,通过克里金或者反距离权重插值法,将网格点位插值到相邻站点间距为800米的更密的点位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611256348.2A CN106646602A (zh) | 2016-12-30 | 2016-12-30 | 一种基于多种震源模型的震后快速生成震动图的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611256348.2A CN106646602A (zh) | 2016-12-30 | 2016-12-30 | 一种基于多种震源模型的震后快速生成震动图的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106646602A true CN106646602A (zh) | 2017-05-10 |
Family
ID=58837281
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611256348.2A Pending CN106646602A (zh) | 2016-12-30 | 2016-12-30 | 一种基于多种震源模型的震后快速生成震动图的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106646602A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109375252A (zh) * | 2018-12-13 | 2019-02-22 | 中国地震局地球物理研究所 | 考虑不同发震构造最大可信地震的地震动参数评价方法 |
CN110244355A (zh) * | 2019-07-25 | 2019-09-17 | 西南交通大学 | 一种基于震源断层模型的脉冲地震动模拟方法 |
CN110688773A (zh) * | 2019-10-14 | 2020-01-14 | 中国电建集团成都勘测设计研究院有限公司 | 一种流域堵塞快速定位的系统及方法 |
CN111060962A (zh) * | 2019-12-31 | 2020-04-24 | 武汉地震科学仪器研究院有限公司 | 一种基于地震仪实测及数值模拟近实时烈度图生成方法 |
CN117058844A (zh) * | 2023-08-03 | 2023-11-14 | 天津君秒安减灾科技有限公司 | 一种地震预警报警方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4766514B2 (ja) * | 2006-02-16 | 2011-09-07 | 株式会社竹中工務店 | 震災リスクの評価方法及び震災リスク評価用の断層モデル形成プログラム |
CN104462814A (zh) * | 2014-12-03 | 2015-03-25 | 沈阳工业大学 | 一种近地表地震动模拟的网格分级方法 |
-
2016
- 2016-12-30 CN CN201611256348.2A patent/CN106646602A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4766514B2 (ja) * | 2006-02-16 | 2011-09-07 | 株式会社竹中工務店 | 震災リスクの評価方法及び震災リスク評価用の断層モデル形成プログラム |
CN104462814A (zh) * | 2014-12-03 | 2015-03-25 | 沈阳工业大学 | 一种近地表地震动模拟的网格分级方法 |
Non-Patent Citations (4)
Title |
---|
陈鲲等: "2014年2月12日新疆于田7.3级地震震动图", 《地质地震》 * |
陈鲲等: "基于地震记录的震动图校正方法研究", 《应用基础与工程科学学报》 * |
陈鲲等: "考虑场地效应的ShakeMap系统研究", 《中国地震》 * |
陈鲲等: "考虑震源破裂过程的青海玉树地震震动图研究", 《中国地震》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109375252A (zh) * | 2018-12-13 | 2019-02-22 | 中国地震局地球物理研究所 | 考虑不同发震构造最大可信地震的地震动参数评价方法 |
CN109375252B (zh) * | 2018-12-13 | 2020-05-05 | 中国地震局地球物理研究所 | 考虑不同发震构造最大可信地震的地震动参数评价方法 |
CN110244355A (zh) * | 2019-07-25 | 2019-09-17 | 西南交通大学 | 一种基于震源断层模型的脉冲地震动模拟方法 |
CN110688773A (zh) * | 2019-10-14 | 2020-01-14 | 中国电建集团成都勘测设计研究院有限公司 | 一种流域堵塞快速定位的系统及方法 |
CN111060962A (zh) * | 2019-12-31 | 2020-04-24 | 武汉地震科学仪器研究院有限公司 | 一种基于地震仪实测及数值模拟近实时烈度图生成方法 |
CN117058844A (zh) * | 2023-08-03 | 2023-11-14 | 天津君秒安减灾科技有限公司 | 一种地震预警报警方法 |
CN117058844B (zh) * | 2023-08-03 | 2024-04-19 | 天津君秒安减灾科技有限公司 | 一种地震预警报警方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106646602A (zh) | 一种基于多种震源模型的震后快速生成震动图的方法 | |
Mai et al. | Hybrid broadband ground-motion simulations: Combining long-period deterministic synthetics with high-frequency multiple S-to-S backscattering | |
Ornthammarath et al. | Probabilistic seismic hazard assessment for Thailand | |
CN105787281B (zh) | 一种海浪有效波高场的融合反演方法和装置 | |
CN105807316B (zh) | 基于振幅叠加的地面观测微地震速度模型校正方法 | |
Courboulex et al. | High‐frequency directivity effect for an M w 4.1 earthquake, widely felt by the population in southeastern France | |
AU2019216706A1 (en) | Iterative ray-tracing for autoscaling of oblique ionograms | |
CN105093319B (zh) | 基于三维地震数据的地面微地震静校正方法 | |
CN102279410A (zh) | 矿山地下开采活动实时监测系统及其方法 | |
Bakun et al. | Epistemic uncertainty in the location and magnitude of earthquakes in Italy from macroseismic data | |
Ugurhan et al. | Strong‐ground‐motion simulation of the 6 April 2009 L’Aquila, Italy, earthquake | |
Palano et al. | Strain and stress fields in the Southern Apennines (Italy) constrained by geodetic, seismological and borehole data | |
Zafarani et al. | Stochastic simulation of strong‐motion records from the 2012 Ahar–Varzaghan Dual Earthquakes, Northwest of Iran | |
Cardona et al. | Crustal seismicity associated to rpid surface uplift at Laguna del Maule Volcanic Complex, Southern Volcanic Zone of the Andes | |
CN104459784A (zh) | 基于单台、双台和双事件数据二维Lg波Q值层析成像方法 | |
CN104122586B (zh) | 一种地震勘探无线节点数据监控方法 | |
Manea et al. | Estimating geophysical bedrock depth using single station analysis and geophysical data in the extra-Carpathian area of Romania | |
Quintero et al. | Near-regional CMT and multiple-point source solution of the September 5, 2012, Nicoya, Costa Rica Mw 7.6 (GCMT) earthquake | |
Hsieh et al. | Efficient inversions for earthquake slip distributions in 3D structures | |
Yang et al. | Towards fast focal mechanism inversion of shallow crustal earthquakes in the Chinese mainland | |
CN107918147A (zh) | 绕射波成像方法和装置 | |
Iyengar et al. | Strong ground motion estimation during the Kutch, India earthquake | |
CN108132140A (zh) | 一种星载光学遥感器全视场柱辐射率确定方法 | |
CN105929445A (zh) | 微破裂向量扫描方法 | |
Konovalov et al. | Testing of a New GMPE Model in Probabilistic Seismic Hazard Analysis for the Sakhalin Region |
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 |
Application publication date: 20170510 |
|
RJ01 | Rejection of invention patent application after publication |