CN111783308A - Method for accurately predicting tunnel surrounding rock displacement - Google Patents
Method for accurately predicting tunnel surrounding rock displacement Download PDFInfo
- Publication number
- CN111783308A CN111783308A CN202010651430.5A CN202010651430A CN111783308A CN 111783308 A CN111783308 A CN 111783308A CN 202010651430 A CN202010651430 A CN 202010651430A CN 111783308 A CN111783308 A CN 111783308A
- Authority
- CN
- China
- Prior art keywords
- model
- regression
- models
- displacement
- tunnel
- 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
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
The invention discloses a method for accurately predicting tunnel surrounding rock displacement, which comprises the steps of collecting displacement data of surface subsidence, vault subsidence, periphery convergence and the like in time after a tunnel is excavated, recording time, and drawing a time-displacement scatter diagram; performing regression analysis simultaneously according to an exponential model, a logarithmic model and a hyperbolic model given in technical specification for monitoring and measuring the railway tunnel (Q/CR 9218-; according to the probability theory and the theory related to mathematical statistics, the three models are respectively subjected to the significance test of a regression model and a regression coefficient; after the regression model and the significance test of the regression coefficient are passed, the model with the minimum residual square sum is used as the regression model; the invention can avoid adopting one type of model of exponential model, logarithmic model and hyperbolic model to carry out regression, can avoid adopting an unobvious model to carry out regression, and can use the model with the best regression effect as the regression model, wherein the best regression effect means that the regression value is closest to the measured value.
Description
Technical Field
The invention belongs to the technical field of tunnel engineering, and particularly relates to a method for accurately predicting tunnel surrounding rock displacement.
Background
For tunnel engineering, displacement monitoring and measuring items such as surface subsidence, vault subsidence and peripheral convergence belong to indispensable measuring items and also belong to important components of new Austrian construction. In the tunnel construction process, time-displacement scatter diagrams are generally drawn according to time for the measured displacements, approved mathematical models are selected for fitting analysis, the mathematical models can be selected according to index models, logarithm models and hyperbolic models given by technical regulations for railway tunnel monitoring and measuring (Q/CR9218-2015), then the maximum or final displacement is subjected to prediction analysis and is compared and analyzed with control reference values given by technical regulations for tunnel monitoring and measuring, and the working states of tunnel surrounding rocks and lining supporting structures are comprehensively analyzed and researched in combination with tunnel construction conditions. If the fitted accumulated displacement curve tends to a fixed value along with time, the tunnel surrounding rock is in a stable state, and the lining supporting structure is reliable; otherwise, indicating that the tunnel surrounding rock and the lining supporting structure are in an unstable state, and taking corresponding remedial measures in time. Therefore, the accurate prediction of the tunnel surrounding rock displacement is of great importance to the safety evaluation of tunnel construction.
In the regression fitting process of tunnel surrounding rock displacement, the following three problems are found by reviewing the existing literature, wherein 1, part of the literature only adopts one type of model in the three types of models to carry out regression; 2. part of the literature was not tested for significance during the regression fitting process; 3. some documents use the above three types of models to perform regression, but select the final model based on the correlation coefficient.
Therefore, it is urgently needed to develop a method for accurately predicting the displacement of the surrounding rock of the tunnel to solve the above problems.
Disclosure of Invention
To solve the problems set forth in the background art described above. The invention provides a method for accurately predicting tunnel surrounding rock displacement.
In order to achieve the purpose, the invention provides the following technical scheme:
a method for accurately predicting tunnel surrounding rock displacement comprises the following steps:
s1, after the tunnel is excavated, collecting displacement data such as surface subsidence, vault subsidence and peripheral convergence in time, recording time, and drawing a time-displacement scatter diagram;
s2, carrying out regression analysis simultaneously according to an exponential model, a logarithmic model and a hyperbolic model given in technical specification Q/CR9218-2015 for railway tunnel monitoring and measuring;
s3, according to the probability theory and the mathematical statistics correlation theory, respectively carrying out significance test on the regression model and the regression coefficient on the three models, and judging whether the three models pass or not, if the three models pass, the step S4 is carried out, and if the three models do not pass, the step S5 is carried out;
s4, after the significance test of the regression model and the regression coefficient is carried out, the model with the minimum residual square sum is used as the regression model;
and S5, kicking off the regression model after not passing the significance test of the regression model and the regression coefficient.
In steps S2, S3, S4, S5: the three analysis processes of adopting three types of displacement duration models to carry out regression simultaneously, carrying out significance test on the regression model and regression coefficients and adopting a model with minimum residual square sum as the regression model are all absent.
Compared with the prior art, the invention has the beneficial effects that:
after the method is adopted, step S2 can avoid adopting one type of model of exponential model, logarithmic model and hyperbolic model for regression, namely can avoid adopting a model with poor regression effect for regression; step S3 may avoid using insignificant models for regression; in step S4, the model with the best regression effect, i.e., the regression value closest to the measured value, may be used as the regression model. The correlation coefficient represents the closeness of the linear correlation between the dependent variable and the independent variable, and the closer the correlation coefficient is to 1, the closer the linear correlation relationship between the dependent variable and the independent variable is. It can be seen that the correlation coefficient cannot represent the degree of deviation between the measured value and the regression value, and for selecting the regression model closest to the measured value, the selection can be performed according to the sum of squares of the residuals, the smaller the sum of squares of the residuals is, the closer the regression value is to the measured value is, and the better the regression effect is.
Drawings
FIG. 1 is a flow chart of a method for accurately predicting tunnel surrounding rock displacement according to the present invention.
Fig. 2 is a time-displacement scatter diagram of the peripheral convergence of the lower dam tunnel DK238+620 according to the present invention.
FIG. 3 is a block diagram of the nonlinear curve fitting module of origin software according to the present invention.
FIG. 4 is a diagram of the new functions of origin software according to the present invention.
FIG. 5 is a block diagram of a dialog for entering a functional form into origin software in accordance with the present invention.
FIG. 6 is a block diagram of the dialog of origin software for performing regression analysis in accordance with the present invention.
FIG. 7 is a diagram of an exponential model obtained by regression in accordance with the present invention.
FIG. 8 is a graph of the regression model of the present invention.
FIG. 9 is a hyperbolic model diagram obtained by regression in accordance with the present invention.
FIG. 10 is a comparison of exponential, logarithmic and hyperbolic models of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention provides the following technical scheme:
as shown in fig. 1, a method for accurately predicting displacement of surrounding rock of a tunnel includes:
s1, after the tunnel is excavated, collecting displacement data such as surface subsidence, vault subsidence and peripheral convergence in time, recording time, and drawing a time-displacement scatter diagram;
s2, carrying out regression analysis simultaneously according to an exponential model, a logarithmic model and a hyperbolic model given in technical specification Q/CR9218-2015 for railway tunnel monitoring and measuring;
the exponential model is:
U=Ae(-B/t);
the logarithmic model is:
U=Alg(1+t)+B;
the hyperbolic model is:
in the above formula: u is a deformation value; A. b is a regression coefficient; t is the observation time (d) of the measurement point.
S3, according to the probability theory and the mathematical statistics correlation theory, respectively carrying out significance test on the regression model and the regression coefficient on the three models, and judging whether the three models pass or not, if the three models pass, the step S4 is carried out, and if the three models do not pass, the step S5 is carried out;
s4, after the significance test of the regression model and the regression coefficient is carried out, the model with the minimum residual square sum is used as the regression model;
and S5, kicking off the regression model after not passing the significance test of the regression model and the regression coefficient.
In steps S2, S3, S4, S5: the three analysis processes of adopting three types of displacement duration models to carry out regression simultaneously, carrying out significance test on the regression model and regression coefficients and adopting a model with minimum residual square sum as the regression model are all absent.
Each of the above steps is explained below with reference to specific examples.
Step 1: after the tunnel is excavated, displacement data such as surface subsidence, vault subsidence and peripheral convergence are collected in time, time is recorded, and a time-displacement scatter diagram is drawn. The method specifically comprises the following steps:
according to the monitoring and measuring frequency given by technical specification of railway tunnel monitoring and measuring (Q/CR9218-2015), displacement data such as ground surface subsidence, vault subsidence and peripheral convergence are collected, instruments such as a high-precision total station can be adopted as a collecting tool, and a time-displacement scatter diagram is drawn. The invention takes the peripheral convergence data of the mileage DK238+620 of a lower dam tunnel of a newly-built Symphytus-Hengchu railway (Chuanchuan segment) as an example, the data is shown in a table 1, and a drawn time-displacement scatter diagram is shown in a table 2.
Table 1 perimeter convergence of lower dam tunnel DK238+620
Step 2: and (4) carrying out regression analysis on the displacement duration models of the exponential model, the logarithmic model and the hyperbolic model simultaneously. The method specifically comprises the following steps:
the nonlinear curve fitting module of origin software (see fig. 3) is adopted, then a new function (see fig. 4) is selected, the function forms of an exponential model, a logarithmic model and a hyperbolic model are respectively input (see fig. 5), finally, the regression analysis can be completed by clicking fitting (see fig. 6), and the exponential model (see fig. 7), the logarithmic model (see fig. 8) and the hyperbolic model (see fig. 9) obtained by regression are obtained.
And step 3: according to the probability theory and the theory related to mathematical statistics, the three models are respectively subjected to the significance test of the regression model and the regression coefficient, and whether the three models pass or not is judged, specifically:
origin software gives the results of the regression model and also the results of the significance test. The results of the significance tests for the exponential model regression model and the regression coefficients are shown in tables 2 and 3, respectively. The results of the significance tests on the logistic model regression model and the regression coefficients are shown in tables 4 and 5, respectively. The results of the significance tests on the hyperbolic model regression model and the regression coefficients are shown in tables 6 and 7, respectively.
TABLE 2 results of significance testing of regression models of exponential models
TABLE 3 significance test results of regression coefficients of exponential models
TABLE 4 results of significance testing of logistic model regression models
TABLE 5 significance test results for regression coefficients of log models
TABLE 6 hyperbolic model regression model significance test results
TABLE 7 hyperbolic model regression coefficient significance test results
As can be seen from tables 2 to 7, the values in the last column are all less than 0.05, i.e., the p-values for the exponential model, the logarithmic model, the hyperbolic model regression model, and the regression coefficients are all less than 0.05. Therefore, both the regression model and the significance test of the regression coefficients were passed.
And 4, step 4: after passing the significance test of the regression model and the regression coefficient, the model with the minimum sum of squared residuals was used as the regression model. The method specifically comprises the following steps:
as can be seen from table 2, when the measured values in table 1 were regressed by the exponential model, the sum of squares of the residuals was 251.86; as can be seen from table 4, when the logarithmic model is used to regress the measured values in table 1, the sum of squares of the residuals is 876.22; as can be seen from table 6, when the hyperbolic model is used to regress the actual measurement values in table 1, the sum of squares of the residuals is 593.76. Namely, the sum of the squares of the residuals is minimum when an exponential model is adopted for regression, namely, the regression value at the moment is closest to the measured value, namely, the displacement predicted based on the model is most accurate. The comparison results of the above three types of models are shown in FIG. 10.
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.
Claims (1)
1. A method for accurately predicting tunnel surrounding rock displacement is characterized by comprising the following steps:
s1, after the tunnel is excavated, collecting displacement data such as surface subsidence, vault subsidence and peripheral convergence in time, recording time, and drawing a time-displacement scatter diagram;
s2, performing regression analysis simultaneously according to an exponential model, a logarithmic model and a hyperbolic model given in the technical specification for monitoring and measuring railway tunnels (Q/CR 9218-2015);
s3, according to the probability theory and the mathematical statistics correlation theory, respectively carrying out significance test on the regression model and the regression coefficient on the three models, and judging whether the three models pass or not, if the three models pass, the step S4 is carried out, and if the three models do not pass, the step S5 is carried out;
s4, after the significance test of the regression model and the regression coefficient is carried out, the model with the minimum residual square sum is used as the regression model;
and S5, kicking off the regression model after not passing the significance test of the regression model and the regression coefficient.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010651430.5A CN111783308A (en) | 2020-07-08 | 2020-07-08 | Method for accurately predicting tunnel surrounding rock displacement |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010651430.5A CN111783308A (en) | 2020-07-08 | 2020-07-08 | Method for accurately predicting tunnel surrounding rock displacement |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111783308A true CN111783308A (en) | 2020-10-16 |
Family
ID=72758170
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010651430.5A Pending CN111783308A (en) | 2020-07-08 | 2020-07-08 | Method for accurately predicting tunnel surrounding rock displacement |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111783308A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115331394A (en) * | 2022-08-30 | 2022-11-11 | 重庆地质矿产研究院 | Method for reducing failure rate of geological disaster early warning system based on key parameter predicted value |
CN115824813A (en) * | 2023-02-23 | 2023-03-21 | 叙镇铁路有限责任公司 | Test device and method for testing surrounding rock plastic zone range caused by tunnel excavation |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05187997A (en) * | 1992-01-10 | 1993-07-27 | Japan Tobacco Inc | Water content measurement device by infrared rays |
US20050234688A1 (en) * | 2004-04-16 | 2005-10-20 | Pinto Stephen K | Predictive model generation |
WO2011089258A2 (en) * | 2010-01-25 | 2011-07-28 | Imec | A variability-aware reliability simulation method of electronic systems |
CN103020736A (en) * | 2012-12-12 | 2013-04-03 | 西安电子科技大学 | Information mining and progress forecasting method based on heterogeneous system integration |
CN103646180A (en) * | 2013-12-19 | 2014-03-19 | 山东大学 | Method for forecasting acute toxicity of organic compounds by building quantitative structure-activity relationship model with quantum chemistry method |
CN105608061A (en) * | 2015-12-23 | 2016-05-25 | 郝士华 | Tunnel monitoring displacement measurement end-section data regression analysis method |
CN106156892A (en) * | 2016-07-13 | 2016-11-23 | 国网福建省电力有限公司 | A kind of method for building up of grid line loss rate forecast model |
CN106600055A (en) * | 2016-12-12 | 2017-04-26 | 华北电力大学(保定) | Wind speed prediction method the basis of self excitation threshold autoregression model |
CN107014683A (en) * | 2017-04-12 | 2017-08-04 | 山东东山王楼煤矿有限公司 | A kind of acquisition methods of earth's surface Rock Displacement Movement parameter |
CN109635375A (en) * | 2018-11-23 | 2019-04-16 | 中国航空工业集团公司沈阳飞机设计研究所 | Aircraft off-the-shelf hardware weight estimation method based on multiple linear regression |
CN109635344A (en) * | 2018-11-14 | 2019-04-16 | 中国航空工业集团公司沈阳飞机设计研究所 | Effectiveness Evaluation Model preparation method and device based on l-G simulation test |
CN109753684A (en) * | 2018-11-29 | 2019-05-14 | 国网江苏省电力有限公司盐城供电分公司 | One kind being used for the modified multiple linear regression modeling method of substation's energy consumption benchmark |
CN109766517A (en) * | 2018-11-29 | 2019-05-17 | 国网江苏省电力有限公司盐城供电分公司 | A kind of energy consumption benchmark modification method for substation energy efficiency assessment |
CN110322371A (en) * | 2019-07-08 | 2019-10-11 | 长沙理工大学 | The area Gao Suntai multiplexing electric abnormality user based on multiple linear regression analysis detects localization method |
CN111123188A (en) * | 2019-12-20 | 2020-05-08 | 国网山东省电力公司电力科学研究院 | Electric energy meter comprehensive verification method and system based on improved least square method |
-
2020
- 2020-07-08 CN CN202010651430.5A patent/CN111783308A/en active Pending
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05187997A (en) * | 1992-01-10 | 1993-07-27 | Japan Tobacco Inc | Water content measurement device by infrared rays |
US20050234688A1 (en) * | 2004-04-16 | 2005-10-20 | Pinto Stephen K | Predictive model generation |
WO2011089258A2 (en) * | 2010-01-25 | 2011-07-28 | Imec | A variability-aware reliability simulation method of electronic systems |
CN103020736A (en) * | 2012-12-12 | 2013-04-03 | 西安电子科技大学 | Information mining and progress forecasting method based on heterogeneous system integration |
CN103646180A (en) * | 2013-12-19 | 2014-03-19 | 山东大学 | Method for forecasting acute toxicity of organic compounds by building quantitative structure-activity relationship model with quantum chemistry method |
CN105608061A (en) * | 2015-12-23 | 2016-05-25 | 郝士华 | Tunnel monitoring displacement measurement end-section data regression analysis method |
CN106156892A (en) * | 2016-07-13 | 2016-11-23 | 国网福建省电力有限公司 | A kind of method for building up of grid line loss rate forecast model |
CN106600055A (en) * | 2016-12-12 | 2017-04-26 | 华北电力大学(保定) | Wind speed prediction method the basis of self excitation threshold autoregression model |
CN107014683A (en) * | 2017-04-12 | 2017-08-04 | 山东东山王楼煤矿有限公司 | A kind of acquisition methods of earth's surface Rock Displacement Movement parameter |
CN109635344A (en) * | 2018-11-14 | 2019-04-16 | 中国航空工业集团公司沈阳飞机设计研究所 | Effectiveness Evaluation Model preparation method and device based on l-G simulation test |
CN109635375A (en) * | 2018-11-23 | 2019-04-16 | 中国航空工业集团公司沈阳飞机设计研究所 | Aircraft off-the-shelf hardware weight estimation method based on multiple linear regression |
CN109753684A (en) * | 2018-11-29 | 2019-05-14 | 国网江苏省电力有限公司盐城供电分公司 | One kind being used for the modified multiple linear regression modeling method of substation's energy consumption benchmark |
CN109766517A (en) * | 2018-11-29 | 2019-05-17 | 国网江苏省电力有限公司盐城供电分公司 | A kind of energy consumption benchmark modification method for substation energy efficiency assessment |
CN110322371A (en) * | 2019-07-08 | 2019-10-11 | 长沙理工大学 | The area Gao Suntai multiplexing electric abnormality user based on multiple linear regression analysis detects localization method |
CN111123188A (en) * | 2019-12-20 | 2020-05-08 | 国网山东省电力公司电力科学研究院 | Electric energy meter comprehensive verification method and system based on improved least square method |
Non-Patent Citations (3)
Title |
---|
葛东旭: "《数据挖掘原理与应用》", 30 April 2020, pages: 245 - 249 * |
郝士华 等: "大断面隧道台阶法开挖净空收敛数据分段回归分析", 《公路交通科技(应用技术版)》 * |
郝士华 等: "大断面隧道台阶法开挖净空收敛数据分段回归分析", 《公路交通科技(应用技术版)》, 31 December 2015 (2015-12-31), pages 1 - 6 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115331394A (en) * | 2022-08-30 | 2022-11-11 | 重庆地质矿产研究院 | Method for reducing failure rate of geological disaster early warning system based on key parameter predicted value |
CN115331394B (en) * | 2022-08-30 | 2023-09-01 | 重庆地质矿产研究院 | Method for reducing fault rate of geological disaster early warning system based on key parameter predicted value |
CN115824813A (en) * | 2023-02-23 | 2023-03-21 | 叙镇铁路有限责任公司 | Test device and method for testing surrounding rock plastic zone range caused by tunnel excavation |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111783308A (en) | Method for accurately predicting tunnel surrounding rock displacement | |
CN102629300A (en) | Step stress accelerated degradation data assessment method based on gray prediction models | |
CN109271319B (en) | Software fault prediction method based on panel data analysis | |
CN111367959B (en) | Zero-time-lag nonlinear expansion Granger causal analysis method | |
CN114417697A (en) | Neural network-based TBM hob abrasion real-time prediction method and system | |
CN105239608A (en) | Landslide displacement prediction method based on wavelet transform-rough set-support vector regression (WT-RS-SVR) combination | |
CN111104989A (en) | Gas pipeline corrosion prediction method and device | |
CN111143768A (en) | Air quality prediction algorithm based on ARIMA-SVM combined model | |
CN118246128B (en) | Reinforced concrete frame structure quality detection method, medium and system | |
CN111143934A (en) | Structural deformation prediction method based on time convolution network | |
Wu et al. | Remaining useful life estimation based on a nonlinear Wiener process model with CSN random effects | |
CN112149296A (en) | Method for judging stability type of hydrological time sequence | |
CN115451803A (en) | Water gate real-time monitoring and early warning method based on fusion of GNSS and measuring robot | |
CN117233347B (en) | Carbon steel spheroidization grade measuring method, system and equipment | |
CN114492074A (en) | Probabilistic damage tolerance assessment analysis method | |
CN110726542A (en) | Analysis method for fatigue life of spring | |
CN116954624B (en) | Compiling method based on software development kit, software development system and server | |
CN102621018A (en) | Online automatic measuring method for rubber hardness | |
CN110895624B (en) | Method for testing consistency of degraded data of accelerated storage and natural storage based on maximum entropy spectrum estimation | |
CN114382490A (en) | Shield tunneling machine cutter head wear assessment and prediction method | |
RU2502974C1 (en) | Method to determine residual resource of technical facilities | |
CN114995338A (en) | Industrial process micro-fault detection method based on normative variable analysis and JS divergence fusion | |
Uzielli et al. | Uncertainty-based characterisation of Troll marine clay | |
CN110889083B (en) | Degraded data consistency checking method based on window spectrum estimation | |
CN115858403B (en) | False alarm rate prediction method of electronic system |
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: 20201016 |