CN117195511B - Quantitative calculation method for initial crust thickness and expansion coefficient - Google Patents

Quantitative calculation method for initial crust thickness and expansion coefficient Download PDF

Info

Publication number
CN117195511B
CN117195511B CN202311063649.3A CN202311063649A CN117195511B CN 117195511 B CN117195511 B CN 117195511B CN 202311063649 A CN202311063649 A CN 202311063649A CN 117195511 B CN117195511 B CN 117195511B
Authority
CN
China
Prior art keywords
crust
thickness
initial
expansion coefficient
crust thickness
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
Application number
CN202311063649.3A
Other languages
Chinese (zh)
Other versions
CN117195511A (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.)
South China Sea Institute of Oceanology of CAS
Original Assignee
South China Sea Institute of Oceanology of CAS
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 South China Sea Institute of Oceanology of CAS filed Critical South China Sea Institute of Oceanology of CAS
Priority to CN202311063649.3A priority Critical patent/CN117195511B/en
Publication of CN117195511A publication Critical patent/CN117195511A/en
Application granted granted Critical
Publication of CN117195511B publication Critical patent/CN117195511B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

In order to solve the problem that the calculation result is inaccurate based on the assumption that the initial crust thickness is a constant value in the calculation of the crust expansion coefficient of the passive land edge or the fractured basin, a new quantitative calculation method for the initial crust thickness and the expansion coefficient is provided, namely, the initial crust thickness and the expansion coefficient are inverted through a uniform limited stretching simulation technology under the constraint of dual data of the current crust thickness and the subsidence period structure subsidence. The method does not assume that the initial crust thickness is a constant value any more, but provides a novel method for jointly inverting the initial crust thickness and the expansion coefficient based on double constraints of the current crust thickness and structural settlement. The new method for calculating the initial crust thickness and the expansion coefficient provided by the invention has important significance for accurately revealing the crust expansion and thinning history, the passive land edge and basin forming process, the hydrocarbon basin sedimentation history, the heat history and the hydrocarbon generation history.

Description

Quantitative calculation method for initial crust thickness and expansion coefficient
Technical Field
The invention relates to the technical field of crust, in particular to a quantitative calculation method for initial crust thickness and expansion coefficient.
Background
The expansion coefficient is an important parameter for quantifying the expansion and thinning degree of the crust of the passive land margin or the fractured basin, and is also the basic data of sedimentation history, heat history and hydrocarbon generation history calculation in oil gas evaluation. By definition, the stretch factor is the ratio of the thickness of the crust before and after stretching, i.e. the ratio of the initial crust thickness to the present crust thickness. However, since the thickness of the crust before stretching and thinning cannot be obtained by direct measurement, most of practical calculation is performed by assuming that the initial crust thickness is a constant value and then calculating according to definition. However, by analyzing the crust structure in different construction contexts, it is found that the crust thickness has a large variation characteristic in space even in the same construction context, and therefore, calculating the crust expansion coefficient based on the assumed initial crust thickness as a constant value brings a large error to the calculation, makes the result inaccurate, and directly affects the evaluation and understanding of the crust expansion thinning history, the passive land and basin formation process, the hydrocarbon basin sedimentation history, the heat history, and the hydrocarbon generation history.
Disclosure of Invention
In order to solve the problem that the calculation result is inaccurate based on the assumption that the initial crust thickness is a constant value in the calculation of the crust expansion coefficient of the passive land edge or the fractured basin, a novel quantitative calculation method for the initial crust thickness and the expansion coefficient is provided, namely the initial crust thickness and the expansion coefficient are inverted by a uniform finite stretching simulation technology under the constraint of dual data of the current crust thickness and the subsidence period structure settlement
In order to achieve the above purpose, the technical scheme of the invention is as follows:
a method for quantitatively calculating initial crust thickness and expansion coefficient, comprising:
s1, establishing an annual stratum grid and a crust structure model in a time domain;
s2, performing paleo-water depth recovery and ablation range definition and ablation stratum recovery to obtain paleo-water depth and ablation stratum thickness;
s3, carrying out deep conversion on the established chronostratigraphic framework and the crust structure model, and solving the thickness of a deposition layer and the thickness of the current crust;
S4, under the control of stratum division, analysis of degraded stratum thickness, lithology, paleo-water depth and sea level, the obtained thickness of the sedimentary deposit is debulked and stripped back, and structural settlement of the actually measured fracture stage is calculated;
S5, giving an initial value to the initial crust thickness, giving an initial strain rate, carrying out uniform limited stretching simulation, calculating structural settlement of a theoretical cracking period, if the error between the structural settlement of the theoretical structure and the structural settlement of an actual cracking period is larger than a given threshold value, gradually adjusting the initial strain rate according to a given step length until the structural settlement of the theoretical structure is calculated to meet the precision, obtaining a rate of change, then solving an expansion coefficient based on the strain rate and the stretching time, and calculating the thickness of the crust after expansion;
S6, comparing the thickness of the crust after stretching and thinning with the thickness of the crust nowadays, and if the difference is smaller than a given threshold value, considering that the initial crust thickness and the stretching coefficient are reasonable; if the difference is greater than the given threshold, the initial crust thickness is required to be continuously adjusted according to the given step length, and the step S5 is repeatedly executed until the error between the crust thickness after the expansion and thinning and the current crust thickness is within a reasonable range, and finally the initial crust thickness and the expansion coefficient of the research area are obtained.
Further, the chronostratigraphic framework and the crustal structure model are established by the following modes:
And carrying out stratum division and interpretation on the two-dimensional and three-dimensional reflection seismic network data under the constraint of drilling and logging data, identifying stratum denudation and earthquake phases, and carrying out full-area tracking on the substrate and the Moholothurian reflection to establish an annual stratum grid and crust structure model.
Further, ancient water depth restoration and ablation range delineation and ablation formation restoration are performed by:
antique depth recovery and ablation range delineation and ablation formation recovery are performed based on drilling lithology, logging interpretation, paleobiology, paleogeography, seismic facies, sedimentary facies data.
Further, the time-depth conversion is carried out on the established chronostratigraphic grid and crust structure model under the control of the well drilling VSP data, the reflection seismic superposition velocity spectrum and the refraction seismic layer velocity material.
Further, the step S3 further includes: the influence of the value of the evaluation speed on the calculation result of the thickness of the earth crust nowadays.
Further, in the step S4, the method further includes: in the back stripping calculation, the influence of lithology, stratum erosion, paleo-water depth and sea level parameters on the result and the errors caused by the influence are fully considered.
Further, in the step S5, the threshold value is 1km.
Further, in the step S6, the threshold value is 0.1km.
Further, in the step S5, an initial crust thickness is given an initial value of 50km.
Further, in the step S5, the initial strain rate is 10 -30.
Compared with the prior art, the invention has the beneficial effects that:
compared with the prior art, the method does not assume that the initial crust thickness is a constant value, but provides a novel method for jointly inverting the initial crust thickness and the expansion coefficient based on double constraints of the current crust thickness and structural settlement. The new method for calculating the initial crust thickness and the expansion coefficient is not only beneficial to analyzing the historic geography and historic construction background before the land edge and the basin fracture, but also has important significance for accurately revealing the crust expansion and thinning history, the passive land edge and basin formation process, the hydrocarbon basin sedimentation history, the heat history and the hydrocarbon generation history.
Drawings
FIG. 1 is a flow chart of the method for quantitatively calculating initial crust thickness and expansion coefficient provided in example 1 of the present invention;
FIG. 2 is a graph of the thickness of the deposit and the thickness of the crystalline crust today;
FIG. 3 is a structural subsidence of a measured open period of the basin;
FIG. 4 is an initial crust thickness of the investigation region;
fig. 5 shows the expansion coefficients of the study area.
Detailed Description
Examples:
the technical scheme of the invention is further described below with reference to the accompanying drawings and examples.
The passive land edge of the north and south China sea is affected by the diving action before stretching, the volcanic arc structure is widely developed, the space of the crust structure is not uniform, and the specific process of applying the method is illustrated by taking the northeast basin of the land edge of the north and south China sea as an embodiment:
s1, carrying out stratum division and explanation on two-dimensional and three-dimensional reflection seismic network data of the Qionsoutheast basin under the constraint of a large amount of drilling and logging data, identifying stratum ablation phases, tracking the whole area of the substrate and Moholothurian reflection, and establishing basin chronostratigraphic grid and crust structure models in a time domain.
S2, carrying out paleo-water depth recovery, formation ablation range definition and ablation formation thickness recovery according to the drilling lithology, logging interpretation, paleobion, paleogeographic, seismic facies, sedimentary and other data of the Qiongtong basin.
S3, carrying out deep conversion on the Qiondong basin stratum grid and the crust structure under the control of the data such as the well drilling VSP data, the reflection seismic superposition velocity spectrum, the refraction seismic layer velocity and the like, solving the thickness of a deposited layer and the thickness of the current crystalline crust (see figure 2), and evaluating the influence of the velocity value on the calculation result of the current crust thickness.
S4, under the control of parameters such as stratum division, analysis of degraded stratum thickness, lithology, ancient water depth, sea level and the like, the compaction and stripping back are carried out on the deposition layer of the southeast basin, and the structural settlement of the actually measured fracture phase of the basin is calculated (see figure 3).
S5, assigning 50km to the initial crust thickness of the southeast agar basin, giving an initial strain rate of 10 -30, carrying out uniform limited stretching simulation, calculating structural settlement of a theoretical open period, and if the error between the structural settlement of the theoretical structure and the structural settlement of the actual peeling back structure is larger than a given threshold value of 1m, gradually adjusting the initial strain rate according to a given step length until the structural settlement of the theoretical structure meets the accuracy, and obtaining a value of a speed change; the expansion coefficient is then solved based on the strain rate and the stretch time, and the thickness of the crust after expansion is calculated.
S6, comparing the thickness of the crust after stretching and thinning with the thickness of the crust nowadays, and if the difference value is smaller than a given threshold value of 0.1km, considering that the initial crust thickness and the stretching coefficient are reasonable; if the difference is greater than the given threshold value of 0.1km, it is necessary to continuously adjust the initial crust thickness by a given step and repeatedly perform step S5 until the error between the stretched thinned crust thickness and the current crust thickness is within a reasonable range, and finally obtain the initial crust thickness of the investigation region (fig. 4) and the stretch coefficient (fig. 5).
It can be seen that the method of the present invention no longer assumes a constant initial crust thickness, but provides a new method for jointly inverting the initial crust thickness and the expansion coefficient based on the double constraints of the current crust thickness and the structural settlement. The new method for calculating the initial crust thickness and the expansion coefficient is not only beneficial to analyzing the historic geography and historic construction background before the land edge and the basin fracture, but also has important significance for accurately revealing the crust expansion and thinning history, the passive land edge and basin formation process, the hydrocarbon basin sedimentation history, the heat history and the hydrocarbon generation history.
The above embodiments are only for illustrating the technical concept and features of the present invention, and are intended to enable those skilled in the art to understand the content of the present invention and implement the same, and are not intended to limit the scope of the present invention. All equivalent changes or modifications made in accordance with the essence of the present invention are intended to be included within the scope of the present invention.

Claims (7)

1. A method for quantitatively calculating an initial crust thickness and expansion coefficient, comprising:
s1, establishing an annual stratum grid and a crust structure model in a time domain;
s2, performing paleo-water depth recovery and ablation range definition and ablation stratum recovery to obtain paleo-water depth and ablation stratum thickness;
s3, carrying out deep conversion on the established chronostratigraphic framework and the crust structure model, and solving the thickness of a deposition layer and the thickness of the current crust;
S4, under the control of stratum division, analysis of degraded stratum thickness, lithology, paleo-water depth and sea level, the obtained thickness of the sedimentary deposit is debulked and stripped back, and structural settlement of the actually measured fracture stage is calculated;
S5, giving an initial value to the initial crust thickness, giving an initial strain rate, carrying out uniform limited stretching simulation, calculating structural settlement of a theoretical cracking period, if the error between the structural settlement of the theoretical structure and the structural settlement of an actual cracking period is larger than a given threshold value, gradually adjusting the initial strain rate according to a given step length until the structural settlement of the theoretical structure is calculated to meet the precision, obtaining a rate of change, then solving an expansion coefficient based on the strain rate and the stretching time, and calculating the thickness of the crust after expansion;
S6, comparing the thickness of the crust after stretching and thinning with the thickness of the crust nowadays, and if the difference is smaller than a given threshold value, considering that the initial crust thickness and the stretching coefficient are reasonable; if the difference value is larger than a given threshold value, the initial crust thickness is required to be continuously adjusted according to a given step length, and the step S5 is repeatedly executed until the error between the crust thickness after the expansion and thinning and the current crust thickness is within a reasonable range, and finally the initial crust thickness and the expansion coefficient of the research area are obtained;
The chronostratigraphic framework and the crust structure model are established in the following manner:
Performing formation division and interpretation on two-dimensional and three-dimensional reflection seismographic network data under the constraint of drilling and logging data
Identifying a layer ablation seismic phase, tracking the whole area of the substrate and the Moholo surface reflection, and establishing an annual stratum grid and crust structure model;
The paleo-water depth recovery and denudation range delineation and denudation formation recovery were carried out by:
antique depth recovery and ablation range delineation and ablation formation recovery are performed based on drilling lithology, logging interpretation, paleobiology, paleogeography, seismic facies, sedimentary facies data.
2. The method of claim 1, wherein the established chronostratigraphic grid and crust structural model is converted in depth under control of borehole VSP data, reflected seismic superposition velocity spectrum, refracted seismic layer velocity profile.
3. The method for quantitatively calculating the initial crust thickness and expansion coefficient according to claim 1, further comprising, in said step S3: the influence of the value of the evaluation speed on the calculation result of the thickness of the earth crust nowadays.
4. The method for quantitatively calculating the initial crust thickness and expansion coefficient according to claim 1, further comprising, in said step S4: in the back stripping calculation, the influence of lithology, stratum erosion, paleo-water depth and sea level parameters on the result and the errors caused by the influence are fully considered.
5. The method for quantitative calculation of initial crust thickness and expansion coefficient according to claim 1, wherein in said step S5, said threshold value is 1km.
6. The method for quantitative calculation of initial crust thickness and expansion coefficient according to claim 1, wherein in said step S6, said threshold value is 0.1km.
7. The method for quantitatively calculating the initial crust thickness and the expansion coefficient according to claim 1, wherein in the step S5, the initial crust thickness is given an initial value of 50km.
CN202311063649.3A 2023-08-23 2023-08-23 Quantitative calculation method for initial crust thickness and expansion coefficient Active CN117195511B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311063649.3A CN117195511B (en) 2023-08-23 2023-08-23 Quantitative calculation method for initial crust thickness and expansion coefficient

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311063649.3A CN117195511B (en) 2023-08-23 2023-08-23 Quantitative calculation method for initial crust thickness and expansion coefficient

Publications (2)

Publication Number Publication Date
CN117195511A CN117195511A (en) 2023-12-08
CN117195511B true CN117195511B (en) 2024-04-30

Family

ID=88998984

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311063649.3A Active CN117195511B (en) 2023-08-23 2023-08-23 Quantitative calculation method for initial crust thickness and expansion coefficient

Country Status (1)

Country Link
CN (1) CN117195511B (en)

Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104200039A (en) * 2014-09-17 2014-12-10 中国石油大学(华东) Quantitative forecasting method of tectonic fissure occurrence
CN104459795A (en) * 2014-12-08 2015-03-25 中国科学院南海海洋研究所 Depth-varying-to-density earth crust extension coefficient thermal calibration gravity anomaly retrieval method
CN106886047A (en) * 2017-02-28 2017-06-23 中国地质大学(北京) A kind of method of receiver function and gravity Inversion CRUSTAL THICKNESS and ripple ratio
CN107015290A (en) * 2017-03-13 2017-08-04 西北大学 A kind of method that reworked garden basin primary deposit looks are recovered
CN109181643A (en) * 2017-03-03 2019-01-11 侯英翔 Metallic ore, nonmetallic ore and coal mine dig lane, when opencast mining, dust-removing method
CN109799533A (en) * 2018-12-28 2019-05-24 中国石油化工股份有限公司 A kind of method for predicting reservoir based on bidirectional circulating neural network
CN110517794A (en) * 2019-08-23 2019-11-29 长安大学 A method of that establishes shale gas reservoir buries-thermal evolution history figure
CN110688728A (en) * 2019-08-20 2020-01-14 中国石油大学(华东) Method for quantitatively analyzing sedimentation characteristics of one-dimensional sediments in time domain and water environment
CN111009179A (en) * 2019-10-23 2020-04-14 上海同继地质科技有限公司 Method and device for determining denudation thickness
CN111475920A (en) * 2020-03-13 2020-07-31 中海石油深海开发有限公司 Method and system for acquiring ancient water depth of deep basin, electronic equipment and storage medium
CN111771144A (en) * 2018-02-28 2020-10-13 沙特阿拉伯石油公司 Locating new hydrocarbon fields and predicting reservoir performance from hydrocarbon migration
AU2020102157A4 (en) * 2020-07-31 2020-10-15 China University Of Geosciences, Beijing Forward Modelling Method For The Subsidence Process Of A Foreland Basin Under A Point Load Condition With Inversion Constraints
CN113970796A (en) * 2020-07-23 2022-01-25 中国石油化工股份有限公司 Method for accurately recovering ancient water depth of sedimentary basin
CN115183739A (en) * 2022-07-13 2022-10-14 中国科学院南海海洋研究所 Method for calculating basin structure settlement based on fault activity weighted extension strain
CN115267909A (en) * 2022-08-09 2022-11-01 中国科学院南海海洋研究所 Surface structure settlement calculation method and device
CN115373024A (en) * 2022-08-09 2022-11-22 中国科学院南海海洋研究所 Method and device for inverting passive land edge crustal structure based on stratum recording settlement
CN115639597A (en) * 2022-10-25 2023-01-24 中国石油大学(华东) Carrying system fine depicting method based on seismic attributes

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104200039A (en) * 2014-09-17 2014-12-10 中国石油大学(华东) Quantitative forecasting method of tectonic fissure occurrence
CN104459795A (en) * 2014-12-08 2015-03-25 中国科学院南海海洋研究所 Depth-varying-to-density earth crust extension coefficient thermal calibration gravity anomaly retrieval method
CN106886047A (en) * 2017-02-28 2017-06-23 中国地质大学(北京) A kind of method of receiver function and gravity Inversion CRUSTAL THICKNESS and ripple ratio
CN109181643A (en) * 2017-03-03 2019-01-11 侯英翔 Metallic ore, nonmetallic ore and coal mine dig lane, when opencast mining, dust-removing method
CN107015290A (en) * 2017-03-13 2017-08-04 西北大学 A kind of method that reworked garden basin primary deposit looks are recovered
CN111771144A (en) * 2018-02-28 2020-10-13 沙特阿拉伯石油公司 Locating new hydrocarbon fields and predicting reservoir performance from hydrocarbon migration
CN109799533A (en) * 2018-12-28 2019-05-24 中国石油化工股份有限公司 A kind of method for predicting reservoir based on bidirectional circulating neural network
CN110688728A (en) * 2019-08-20 2020-01-14 中国石油大学(华东) Method for quantitatively analyzing sedimentation characteristics of one-dimensional sediments in time domain and water environment
CN110517794A (en) * 2019-08-23 2019-11-29 长安大学 A method of that establishes shale gas reservoir buries-thermal evolution history figure
CN111009179A (en) * 2019-10-23 2020-04-14 上海同继地质科技有限公司 Method and device for determining denudation thickness
CN111475920A (en) * 2020-03-13 2020-07-31 中海石油深海开发有限公司 Method and system for acquiring ancient water depth of deep basin, electronic equipment and storage medium
CN113970796A (en) * 2020-07-23 2022-01-25 中国石油化工股份有限公司 Method for accurately recovering ancient water depth of sedimentary basin
AU2020102157A4 (en) * 2020-07-31 2020-10-15 China University Of Geosciences, Beijing Forward Modelling Method For The Subsidence Process Of A Foreland Basin Under A Point Load Condition With Inversion Constraints
CN111931367A (en) * 2020-07-31 2020-11-13 中国地质大学(北京) Forward modeling method for foreland basin settlement process under inversion constraint point load condition
CN115183739A (en) * 2022-07-13 2022-10-14 中国科学院南海海洋研究所 Method for calculating basin structure settlement based on fault activity weighted extension strain
CN115267909A (en) * 2022-08-09 2022-11-01 中国科学院南海海洋研究所 Surface structure settlement calculation method and device
CN115373024A (en) * 2022-08-09 2022-11-22 中国科学院南海海洋研究所 Method and device for inverting passive land edge crustal structure based on stratum recording settlement
CN115639597A (en) * 2022-10-25 2023-01-24 中国石油大学(华东) Carrying system fine depicting method based on seismic attributes

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
下地壳流对南海北缘白云凹陷地壳伸展变形的影响;王晓芳等;《海洋地质与第四纪地质》;20071231;第37卷(第6期);全文 *
南海及邻域基础地质问题与中国大洋钻探选址;尉建功等;《地质学报》;20220831;第96卷(第8期);全文 *
被动大陆边缘张-破裂过程与岩浆活动:南海的归属;孙珍等;《地球科学》;20210331;第46卷(第3期);全文 *

Also Published As

Publication number Publication date
CN117195511A (en) 2023-12-08

Similar Documents

Publication Publication Date Title
Shen et al. Crustal and uppermost mantle structure beneath the United States
Pfiffner et al. The Swiss Alps and their peripheral foreland basin: Stratigraphic response to deep crustal processes
US8868348B2 (en) Well constrained horizontal variable H-V curve constructing method for seismic wave velocity field construction
CN105204069B (en) A kind of Eroded Thickness restoration methods
CN105386756B (en) A method of brittle formation porosity is calculated using dependent variable
CN113031068B (en) Reflection coefficient accurate base tracking prestack seismic inversion method
CN107966732A (en) The seismic properties change rate acquiring method being oriented to based on space structure
CN102967881A (en) Geological small-layer data deep-time conversion method and device for seismic data interpretation
CN113740901A (en) Land seismic data full-waveform inversion method and apparatus based on complex undulating surface
CN107942379A (en) A kind of method for improving complex fault block rate pattern precision
CN111025388B (en) Multi-wave combined prestack waveform inversion method
CN111175825A (en) Depth domain speed modeling method
Xia et al. Correspondences among lacustrine fluctuations, climate changes and the Milankovitch cycles in the Paleogene through tracking onlap points and correlating palaeontology in Liaozhong Depression, Bohai Bay Basin, NE China
Ding et al. Seismic physical modeling study for a shale oil reservoir with complex lithofacies and meter-scale structures
CN117195511B (en) Quantitative calculation method for initial crust thickness and expansion coefficient
CN111413737A (en) Method for explaining seismic horizon to deposition unit in close well pattern area
Prather et al. Stratigraphic analysis of XES02: Implications for the sequence stratigraphic paradigm
CN115183739A (en) Method for calculating basin structure settlement based on fault activity weighted extension strain
CN114563816B (en) Method and device for establishing earthquake interpretation velocity model in oil and gas reservoir evaluation stage
Xu et al. An Improved Fracture Seismic Method for identifying the drilling targets of medium-deep geothermal resources: A case study on heishan geothermal area
CN117094152B (en) Basin sedimentation history simulation method and device for coupling sedimentation
CN111596348A (en) Stratum velocity correction method and device based on acoustic time difference logging data
WO2020070534A1 (en) A method for computing a thermal parameter of a passive margin area and associated system
CN114706141B (en) Underground meandering stream zone quantitative characterization method based on sedimentation analogy
Qin et al. A planar anomaly trend correction time-depth conversion method based on well-seismic constrained tomography velocity inversion

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