CN111859632B - Rock physical model construction method and processing terminal for hydrate reservoir - Google Patents

Rock physical model construction method and processing terminal for hydrate reservoir Download PDF

Info

Publication number
CN111859632B
CN111859632B CN202010625752.2A CN202010625752A CN111859632B CN 111859632 B CN111859632 B CN 111859632B CN 202010625752 A CN202010625752 A CN 202010625752A CN 111859632 B CN111859632 B CN 111859632B
Authority
CN
China
Prior art keywords
hydrate
water
modulus
gas
bulk modulus
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
CN202010625752.2A
Other languages
Chinese (zh)
Other versions
CN111859632A (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.)
Guangzhou Marine Geological Survey
Original Assignee
Guangzhou Marine Geological Survey
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 Guangzhou Marine Geological Survey filed Critical Guangzhou Marine Geological Survey
Priority to CN202010625752.2A priority Critical patent/CN111859632B/en
Publication of CN111859632A publication Critical patent/CN111859632A/en
Application granted granted Critical
Publication of CN111859632B publication Critical patent/CN111859632B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V20/00Geomodelling in general
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention relates to a petrophysical model construction method and a processing terminal of a hydrate reservoir, wherein the method comprises the following steps: step 1: generating an initial self-consistent approximate medium, and calculating the bulk modulus and the shear modulus of the self-consistent approximate medium according to a formula; step 2: generating solid-phase rock, and calculating the bulk modulus and the shear modulus of the solid-phase rock according to a formula; step 3: calculating the bulk modulus and shear modulus of the dry rock; step 4: and adding the biphase fluid to generate a petrophysical model, and calculating attenuation factor parameters of the petrophysical model according to a formula. According to the invention, a self-adaptive theory, an equivalent differential theory, an effective medium theory and a spot saturation model are combined, a petrophysical model aiming at a marine region stratum rich in argillaceous and containing natural gas hydrate is established, and hydrate reservoir attenuation research of a subsequent seismic frequency band can be effectively carried out according to the petrophysical model.

Description

Rock physical model construction method and processing terminal for hydrate reservoir
Technical Field
The invention relates to the technical field of natural gas hydrate exploration, in particular to a rock physical model construction method and a processing terminal of a hydrate reservoir.
Background
Along with the increase of the hydrate content, the increase of the longitudinal wave speed and the transverse wave speed is often caused, and the change of attenuation is also caused, so that the speed dispersion and attenuation characteristic analysis of the hydrate-containing sediment layer is helpful for estimating the hydrate content. Whereas the petrophysical model based on hydrate sedimentary formations (i.e. reservoirs) can analyze the velocity dispersion and attenuation characteristics of the hydrate, the existing petrophysical model for researching sea-area mud-rich high pore hydrate reservoirs is missing or under-estimated and inaccurate in the seismic frequency band (low frequency, less than 500 Hz). The existing research on the petrophysical model includes Biot theory (Chinese is also called than-Australian consolidation theory), jet flow mechanism and BISQ theory combining the two, but the above theories cannot well describe the seismic attenuation characteristics in the seismic frequency band (low frequency, less than 500 Hz), namely cannot well describe the attenuation characteristics of the hydrate reservoir. Therefore, it is desirable to construct a petrophysical model of a hydrate reservoir so that the attenuation characteristics of the hydrate can be analyzed in the seismic frequency band.
Disclosure of Invention
Aiming at the defects of the prior art, one of the purposes of the invention is to provide a method for constructing a petrophysical model of a hydrate reservoir, which can solve the problem of constructing the petrophysical model of the hydrate reservoir;
it is a further object of the present invention to provide a processing terminal that is capable of solving the problem of building a petrophysical model of a hydrate reservoir.
The technical scheme for realizing one of the purposes of the invention is as follows: a method for constructing a petrophysical model of a hydrate reservoir, comprising the steps of:
step 1: generating a self-consistent approximate medium comprising mudstone and hydrate, and calculating the bulk modulus K of the self-consistent approximate medium according to formula (1) SC And shear modulus G SC
Wherein K is ma Represents the bulk modulus of mudstone, G ma Represents the shear modulus of mudstone, S in Representing hydratesContent of K in Represents the bulk modulus of the hydrate, G in Shear modulus of hydrate, P, Q and beta ma Are intermediate variables, and alpha represents the aperture aspect ratio of the self-consistent approximate medium;
step 2: quartz is added on the basis of the self-consistent approximate medium to generate solid-phase rock, and the bulk modulus K of the solid-phase rock is calculated according to the formula (2) dem And shear modulus G dem
Wherein S represents quartz content, K sat And G sat Respectively representing the bulk modulus and the shear modulus, P, of quartz 1 、Q 1 And beta m Are all intermediate variables, alpha 1 Represents the pore aspect ratio of solid-phase rock;
step 3: calculating the bulk modulus K of the dry rock according to the formula (3) dry And shear modulus G dry
Wherein ρ is b Representing the bulk density ρ w Represents the water density, g represents the gravitational acceleration, D represents the sea bottom depth of the hydrate reservoir, n represents the average connection coefficient of the particles,indicating porosity->Represents critical porosity, v represents poisson's ratio, P 2 Representing the effective pressure, Z represents an intermediate variable;
step 4: adding a fluid on the basis of the step 3, wherein the fluid comprises water and gas to obtain a petrophysical model, and calculating according to formulas (4) and (5) to obtain an attenuation factor Q p
Wherein i=1, 2, d 1 And d 2 Representing reservoir thickness of water and gas, respectively, P 3 And P 4 All are the intermediate variables of the two-way valve,andcomplex moduli of propagation of P-waves of water and gas in orthogonal directions, respectively, I 1 And I 2 Respectively represent the impedance of water and gas and the impedance of slow longitudinal wave, r 1 And r 2 The ratio of the fast longitudinal wave fluid tension to the total normal pressure of water and gas respectively,
kappa denotes permeability, omega denotes angular frequency, phi denotes sediment porosity, M i As intermediate variable, alpha 2 As intermediate variable, k i As an intermediate variable, the number of the variables,represents the bulk modulus, eta of water 1 Represents the viscosity of water, which is constant; />Represents the bulk modulus, eta of the gas 2 The viscosity of the gas is expressed as a constant,
obtaining attenuation factor Q p Thereby obtaining the petrophysical model.
The second technical scheme for realizing the purpose of the invention is as follows: a processing terminal comprising, a memory for storing program instructions;
and a processor for executing the program instructions to perform the steps of the petrophysical model building method of the hydrate reservoir.
The beneficial effects of the invention are as follows: according to the invention, a self-adaptive theory, an equivalent differential theory, an effective medium theory and a spot saturation model are combined, a petrophysical model aiming at a marine region stratum rich in argillaceous and containing natural gas hydrate is established, and hydrate reservoir attenuation research of a subsequent seismic frequency band can be effectively carried out according to the petrophysical model.
Drawings
FIG. 1 is a flow chart of a preferred embodiment of the present invention;
FIG. 2 is a schematic of a hydrate reservoir fully saturated with water in the sediment pores;
FIG. 3 is a schematic of a hydrate reservoir with incompletely saturated water in the sediment pores;
FIG. 4 is a schematic diagram comparing the present invention and other models with real data;
fig. 5 is a schematic structural diagram of a processing terminal.
Detailed description of the preferred embodiments
The invention will be further described with reference to the accompanying drawings and detailed description:
prior to describing the implementation of model construction, the relevant background art of geologic models is described. Currently, subterranean formations are generally considered to be composed of solid-phase rock, including natural gas hydrates (hereinafter referred to as hydrates), quartz, and most mudstones, pores, and fluids in the pores, including water and gas. The invention also provides a model construction method based on the theory.
As shown in fig. 1-4, a method for constructing a petrophysical model of a hydrate reservoir includes the following steps:
step 1: generating an initial self-consistent approximate medium, wherein the self-consistent approximate medium comprises mudstone and hydrate, and calculating the bulk modulus K of the self-consistent approximate medium according to the following formula SC And shear modulus G SC
Wherein K is ma Represents the bulk modulus of mudstone, G ma Represents the shear modulus of mudstone, S in Represents the hydrate content, K in Represents the bulk modulus of the hydrate, G in Shear modulus of hydrate, P, Q and beta ma Are intermediate variables, and α represents the aperture aspect ratio of the self-consistent approximation medium. The above formula is based on self-consistent approximation theory.
Step 2: quartz is added on the basis of the self-consistent approximate medium in the step 1 to generate solid-phase rock, and the bulk modulus K of the solid-phase rock is calculated according to the following formula dem And shear modulus G dem
Wherein S represents quartz content, K sat And G sat Respectively representing the bulk modulus and the shear modulus, P, of quartz 1 、Q 1 And beta m Are all intermediate variables, alpha 1 Representing the pore aspect ratio of the solid-phase rock. The above formula is based on the theory of equivalent difference.
Step 3: the bulk modulus K of the dry rock is calculated according to the following formula dry And shear modulus G dry
Wherein ρ is b Representing the bulk density of the dry rock ρ w Represents water density, g represents gravitational acceleration, D represents the sea depth of the hydrate reservoir, i.e. the sea depth of the research target object, the hydrate reservoirs in different areas are researched, the sea depths are different, n represents average connection coefficient of particles, and is constant,representing the porosity of the dry rock, < >>Represents critical porosity, v represents poisson's ratio, P 2 Representing the effective pressure, Z represents an intermediate variable;
step 4: generating a petrophysical model, and adding fluid on the basis of the step 3, wherein the fluid comprises water and gas, and the added fluid is a two-phase fluid, namely the water and the gas are contained in pores, so as to obtain the petrophysical model. Calculating the attenuation factor Q of the final petrophysical model according to the plaque saturation model p Thereby determining the final petrophysical model. Attenuation factor Q p Can be calculated by the following formula:
wherein d 1 And d 2 Representing reservoir thickness of water and gas, respectively, P 3 And P 4 All are the intermediate variables of the two-way valve,and->Complex moduli of propagation of P-waves of water and gas in orthogonal directions, respectively, I 1 And I 2 Respectively represent the impedance of water and gas and the impedance of slow longitudinal wave, r 1 And r 2 The ratio of the fast longitudinal wave fluid tension to the total normal pressure of water and gas respectively,
in the above formula, the relevant parameters of the two-layer medium (water and gas) are calculated by the following formula:
wherein κ represents permeability, ω represents angular frequency, φ represents sediment porosity, M i As intermediate variable, alpha 2 As intermediate variable, k i As an intermediate variable, the number of the variables,represents the bulk modulus, eta of water 1 Represents the viscosity of water, which is constant; />Represents the bulk modulus, eta of the gas 2 The viscosity of the gas is shown as a constant.
Obtaining attenuation factor Q p Thereby constructing a petrophysical model. The seismic frequency band of the petrophysical model obtained by the method is below 500Hz, and belongs to a low-frequency petrophysical model.
The rock physical model obtained by the method is consistent with the actual situation and can be applied to analysis and research such as hydrate reservoir attenuation and the like.
As shown in fig. 2 and 3, fig. 2 is a schematic of a hydrate reservoir fully saturated with water in the sediment pores, and fig. 3 is a schematic of a hydrate reservoir not fully saturated with water (i.e., water and gas) in the sediment pores. In fig. 2, it is assumed that the sediment has a muddy content c=60%, a sediment porosity Φ=50%, a pore aspect ratio α=0.05, and a water content S w =100%, natural gas hydrate content S h Linearly increasing from 10% to 80%. As the frequency f (logarithmic on the abscissa in the figure) increases gradually from 1Hz to 1000Hz, the longitudinal wave attenuation I/Q (on the ordinate in the figure) increases non-linearly to the peak and then decreases gradually, and as the hydrate content increases, the attenuation also increases and the attenuation peak moves gradually to a lower frequency, and when the hydrate content reaches a certain level, the attenuation slowly tends to decrease gradually, i.e., flattens out from the view of the graph.
In fig. 3, it is assumed that the sediment has a muddy content c=60%, a sediment porosity Φ=50%, a pore aspect ratio α=0.05, and a water content S w =90%, gas content S g =10% natural gas hydrate content S h From 0 to 60% (i.e., represented by sh=0.0, 0.2, 0.4, 0.6 in fig. 3). As the frequency f (logarithmic in the graph) increases gradually from 1Hz to 1000Hz, the longitudinal wave attenuation I/Q (ordinate in the graph) increases nonlinearly to the peak valueGradually decreasing, the peak occurs at a corresponding frequency higher than that of a saturated deposit, i.e., the peak of fig. 3 occurs at a higher frequency than that of the peak of fig. 2. And as the hydrate content increases, the attenuation also increases, and the attenuation peak gradually moves toward lower frequency, when the hydrate content reaches a certain level, the attenuation slowly moves toward lower frequency at a significantly faster speed, similar to fig. 2, but compared to fig. 2.
As shown in FIG. 4, another practical example application is shown in Guerin & Goldberg 2005 published paper (Guerin G and Goldberg D2005. Modeling of acoustic wave dissipation in gas hydrate-modeling segments. Geochemistry geosystems 6 (7): 1-16), which is referred to as a reference paper, where measured data sources are compared with Mallik 5L-38 gas hydrate test wells in the large West. Compared with the model proposed by the reference paper (shown as a dotted line, guerin & Goldberg, 2005), the method provided by the invention (shown as a solid line, this student) and the measured data (shown as dots, malik 5L-38 (Guerin & Goldberg, 2005)) are closer, that is, the petrophysical model obtained by the method is closer to the real condition, and the abscissa GH Saturation in FIG. 4 represents the hydrate Saturation.
As shown in fig. 5, the present invention also relates to a processing terminal 100 of an entity apparatus implementing the above method, which includes,
a memory 101 for storing program instructions;
a processor 102 for executing the program instructions to perform steps in the petrophysical model construction method of the hydrate reservoir.
The embodiment disclosed in the present specification is merely an illustration of one-sided features of the present invention, and the protection scope of the present invention is not limited to this embodiment, and any other functionally equivalent embodiment falls within the protection scope of the present invention. Various other corresponding changes and modifications will occur to those skilled in the art from the foregoing description and the accompanying drawings, and all such changes and modifications are intended to be included within the scope of the present invention as defined in the appended claims.

Claims (2)

1. The rock physical model construction method of the hydrate reservoir is characterized by comprising the following steps of:
step 1: generating a self-consistent approximate medium comprising mudstone and hydrate, and calculating the bulk modulus K of the self-consistent approximate medium according to formula (1) SC And shear modulus G SC
Wherein K is ma Represents the bulk modulus of mudstone, G ma Represents the shear modulus of mudstone, S in Represents the hydrate content, K in Represents the bulk modulus of the hydrate, G in Shear modulus of hydrate, P, Q and beta ma Are intermediate variables, and alpha represents the aperture aspect ratio of the self-consistent approximate medium;
step 2: quartz is added on the basis of the self-consistent approximate medium to generate solid-phase rock, and the bulk modulus K of the solid-phase rock is calculated according to the formula (2) dem And shear modulus G dem
Wherein S represents quartz content, K sat And G sat Respectively representing the bulk modulus and the shear modulus, P, of quartz 1 、Q 1 And beta m Are all intermediate variables, alpha 1 Represents the pore aspect ratio of solid-phase rock;
step 3: calculating the bulk modulus K of the dry rock according to the formula (3) dry And shear modulus G dry
Wherein ρ is b Representing the bulk density ρ w Represents the water density, g represents the gravitational acceleration, D represents the sea bottom depth of the hydrate reservoir, n represents the average connection coefficient of the particles,indicating porosity->Represents critical porosity, v represents poisson's ratio, P 2 Representing the effective pressure, Z represents an intermediate variable;
step 4: adding a fluid on the basis of the step 3, wherein the fluid comprises water and gas to obtain a petrophysical model, and calculating according to formulas (4) and (5) to obtain an attenuation factor Q p
Wherein i=1, 2, d 1 And d 2 Representing reservoir thickness of water and gas, respectively, P 3 And P 4 All are the intermediate variables of the two-way valve,and->Complex moduli of propagation of P-waves of water and gas in orthogonal directions, respectively, I 1 And I 2 Respectively represent the impedance of water and gas and the impedance of slow longitudinal wave, r 1 And r 2 The ratio of the fast longitudinal wave fluid tension to the total normal pressure of water and gas respectively,
kappa denotes permeability, omega denotes angular frequency, phi denotes sediment porosity, M i As intermediate variable, alpha 2 As an intermediate variable, the number of the variables,k i as an intermediate variable, the number of the variables,represents the bulk modulus, eta of water 1 Represents the viscosity of water, which is constant; />Represents the bulk modulus, eta of the gas 2 The viscosity of the gas is expressed as a constant,
obtaining attenuation factor Q p Thereby obtaining the petrophysical model.
2. A processing terminal, characterized in that it comprises,
a memory for storing program instructions;
a processor for executing the program instructions to perform the steps of the petrophysical model construction method of a hydrate reservoir of claim 1.
CN202010625752.2A 2020-07-01 2020-07-01 Rock physical model construction method and processing terminal for hydrate reservoir Active CN111859632B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010625752.2A CN111859632B (en) 2020-07-01 2020-07-01 Rock physical model construction method and processing terminal for hydrate reservoir

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010625752.2A CN111859632B (en) 2020-07-01 2020-07-01 Rock physical model construction method and processing terminal for hydrate reservoir

Publications (2)

Publication Number Publication Date
CN111859632A CN111859632A (en) 2020-10-30
CN111859632B true CN111859632B (en) 2023-09-22

Family

ID=72989018

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010625752.2A Active CN111859632B (en) 2020-07-01 2020-07-01 Rock physical model construction method and processing terminal for hydrate reservoir

Country Status (1)

Country Link
CN (1) CN111859632B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112946783B (en) * 2021-01-29 2024-02-06 中国石油大学(北京) Hydrate saturation determination method, device and equipment
CN113552630B (en) * 2021-08-13 2022-03-04 广州海洋地质调查局 Unconsolidated formation permeability prediction method based on elastic impedance and processing terminal
CN113466963B (en) * 2021-08-13 2022-02-01 广州海洋地质调查局 Fracture filling II type hydrate saturation estimation method and processing terminal
CN113791457B (en) * 2021-09-08 2022-11-15 中国海洋石油集团有限公司 Method and device for calculating rock skeleton modulus of natural gas hydrate reservoir

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106842364A (en) * 2015-12-04 2017-06-13 中国石油化工股份有限公司 Unconformity biogenic limestone seismic reservoir petrophysical model method for building up
CN109582923A (en) * 2018-11-06 2019-04-05 广州海洋地质调查局 A kind of method and processing terminal based on two-parameter calculated hydration object saturation degree
CN110909487A (en) * 2019-12-06 2020-03-24 西南石油大学 Method for establishing fine anisotropic shale rock physical model

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10534871B2 (en) * 2011-03-09 2020-01-14 Schlumberger Technology Corporation Method and systems for reservoir modeling, evaluation and simulation
US9482771B2 (en) * 2012-11-05 2016-11-01 Fugro Marine Geoservices, Inc. Method of indicating the presence of gas hydrate and shallow gas in deepwater environment

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106842364A (en) * 2015-12-04 2017-06-13 中国石油化工股份有限公司 Unconformity biogenic limestone seismic reservoir petrophysical model method for building up
CN109582923A (en) * 2018-11-06 2019-04-05 广州海洋地质调查局 A kind of method and processing terminal based on two-parameter calculated hydration object saturation degree
CN110909487A (en) * 2019-12-06 2020-03-24 西南石油大学 Method for establishing fine anisotropic shale rock physical model

Also Published As

Publication number Publication date
CN111859632A (en) 2020-10-30

Similar Documents

Publication Publication Date Title
CN111859632B (en) Rock physical model construction method and processing terminal for hydrate reservoir
CN105468886B (en) The method for calculating strata pressure based on petrophysical parameters
CN103645509B (en) The inverting of compact reservoir pore components and S-Wave Velocity Predicted Method
CN103643949B (en) A kind of reservoir contains quantitative forecast method and the device of oil gas
EP2616848A1 (en) Method of predicting the pressure sensitivity of seismic velocity within reservoir rocks
CN113376709B (en) Method for predicting reservoir natural gas hydrate saturation by using logging data
CN107831540A (en) The direct new method for extracting of reservoir physical parameter
CN117055114B (en) Quantitative analysis method for free gas saturation of reservoir sediment
CN114114409B (en) Petrophysical modeling method for sea natural gas hydrate, electronic equipment and medium
Aljazaeri et al. Modeling of thermal and burial histories for selected deep formations in the Middle-Jurassic to Lower Cretaceous in Siba Gas field, Southern Iraq
CN115857006B (en) Submarine acoustic and physical parameter detection method, medium and system
Chouet et al. Seismic source dynamics of gas‐piston activity at Kı̄lauea Volcano, Hawai ‘i
CN110967742A (en) Porosity inversion method and system
CN113534263B (en) Oil-gas saturation prediction method independent of logging information
Liu et al. Rock-physics models of hydrate-bearing sediments in permafrost, Qilian Mountains, China
CN115201915A (en) Ocean wide azimuth data deep dessert prediction method, storage medium and equipment
CN111650644B (en) Quantitative prediction method for salt cemented sandstone
CN110297280B (en) Prediction method for carbonate rock overpressure spatial distribution characteristics
CN110879424A (en) Fluid property identification method and device
CN113552630B (en) Unconsolidated formation permeability prediction method based on elastic impedance and processing terminal
CN118094919B (en) Method, device, equipment and storage medium for estimating gas saturation of sand body
CN107764697A (en) Gas potential detection method based on the progressive equation non-linear inversion of pore media
CN115079261B (en) Compact sandstone gas reservoir evaluation method based on multiparameter quantitative interpretation template
CN114076983B (en) Earthquake prediction method and device based on tight sandstone oil and gas reservoir effective reservoir
CN109425904B (en) Method for acquiring pore pressure of carbonate rock stratum

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