CN105372701A - Method for detecting soil liquefaction degree of port by using seismic waves - Google Patents

Method for detecting soil liquefaction degree of port by using seismic waves Download PDF

Info

Publication number
CN105372701A
CN105372701A CN201510918325.2A CN201510918325A CN105372701A CN 105372701 A CN105372701 A CN 105372701A CN 201510918325 A CN201510918325 A CN 201510918325A CN 105372701 A CN105372701 A CN 105372701A
Authority
CN
China
Prior art keywords
seismic
liquefaction
wave
degree
soil property
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.)
Granted
Application number
CN201510918325.2A
Other languages
Chinese (zh)
Other versions
CN105372701B (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.)
China Harbour Engineering Co Ltd
Original Assignee
China Harbour Engineering Co Ltd
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 China Harbour Engineering Co Ltd filed Critical China Harbour Engineering Co Ltd
Priority to CN201510918325.2A priority Critical patent/CN105372701B/en
Publication of CN105372701A publication Critical patent/CN105372701A/en
Application granted granted Critical
Publication of CN105372701B publication Critical patent/CN105372701B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

The invention discloses a method for detecting the soil liquefaction degree of a port by using seismic waves. The method comprises the following steps that: a seismic wave mechanical source and a seismic wave receiver are distributed in an area to be detected of the port; clutter in vibration information can be filtered out; seismic wave spectrum strength and seismic wave longitudinal wave velocity can be obtained through calculation and analysis; basic data are inputted to a data processing device; the soil density of the area to be detected, seismic wave spectrum intensity and the dynamic shear stress ratio of the area to be detected can be calculated according to the seismic wave longitudinal wave velocity, and liquefaction resistance rate is calculated, and a calculation formula of the liquefaction resistance rate is described in the descriptions of the invention, and therefore, the liquefaction degree of the area to be detected can be judged. With the method for detecting the soil liquefaction degree of the port by using the seismic waves of the invention adopted, the soil liquefaction degree of the area to be detected of the port can be obtained fast and conveniently.

Description

The method of a kind of seismic event detection port area soil property Degree of Liquefaction
Technical field
The present invention relates to a kind of geological exploration method.More particularly, the present invention relates to a kind of method utilizing seismic event to detect port area soil property Degree of Liquefaction.
Background technology
Along with the continuous expansion of market abroad, port engineering construction area worldwide progressively expands.For the port structure being positioned at distribution of strong earthquakes band, the research of its seismic design progressively becomes the Focal point and difficult point in port engineering design aspect.In run in early stage is in the past implemented, seismic design is the important bottleneck of domestic design unit going global.For an overseas port engineering project, often need the seismic fortification criterion and the specification that relate to the every countries such as the U.S., Europe and engineering location or area, carry out the aseismic analysis of engineering structure by referring to all standard criterion.In addition, when carrying out seismic design, except meeting designing requirement that specification specifies, the on-site background of the side's of design incorporation engineering, construction requirement and earthquake history etc. are needed to adopt more detailed and complicated Seismic Design Method toward contact, when saturated loose sand (comprising some silt) is subject to vibrating, if pore water has little time to discharge, pore water pressure constantly increases, and effective stress reduces gradually.When effective stress is reduced to zero, the state that the soil body can produce sharply changes, and becomes liquid the same, completely loses shearing strength and load-bearing capacity, cause ground and building destruction, the liquefaction phenomenon that Here it is often says.Liquefaction can occur at any position in saturated sand layers, can on ground, also can at a certain depth of subsurface, depend on state and the Vibration Condition of sand, the layer of sand did not sometimes liquefy originally, but because layer of sand liquefies below, hydrostatic excess pressure dissipates with the upwards flowing of water, if at this moment hydraulic gradient is excessive, upward water flow may make the seepage erosion of layer of sand generation above and unstability, even if or do not reach unstability, but the load-bearing capacity of soil layer above also can be made to decline.The construction of anti-Degree of Liquefaction to harbour of harbour soil property plays vital reference role.
Summary of the invention
An object of the present invention is to solve at least the problems referred to above, and the advantage will illustrated at least is below provided.
A further object of the invention is to provide a kind of method utilizing seismic event to detect port area soil property Degree of Liquefaction, and it can be quick, knows region to be measured, harbour soil property Degree of Liquefaction easily.
In order to realize according to these objects of the present invention and other advantage, provide a kind of method utilizing seismic event to detect port area soil property Degree of Liquefaction, it comprises the following steps:
Step one, lays seismic event mechanical source and seismic receiver in region to be measured, harbour, and on described seismic receiver access data harvester and data processing equipment;
Step 2, described data collector gathers the ground of seismic receiver collection and the vibration information of underground, the clutter in vibration information described in described data processing equipment filtering, and computational analysis obtains seismic event wave spectrum intensity and seismic event longitudinal wave velocity;
Step 3, inputs described data processing equipment by basic data; Described basic data comprises: focal point location expression, soil property sampling description, stacking fold, seismic wave frequency spectrum analysis, wave field separation, one or more in rock & soil mechanical parameter;
Step 4, utilize seismic event longitudinal wave velocity to calculate region to be measured soil property density, formula is as follows:
wherein VP is longitudinal wave velocity, and K is Lame's constant, and ∧ is shearing elasticity modulus;
Utilize seimic wave velocity spectral intensity value to calculate seismic wave reaction spectral intensity, computing formula is as follows:
wherein, Sv is velocity spectrum intensity level, and T is the cycle;
Calculate the cyclic shear stress ratio in region to be measured, computing formula is as follows:
Wherein,
α={0.0052(α max/SI)-0.0163}10 -2
b={0.091O(α max/SI)+O.0787}10 -2
Step 5, calculate liquefaction opposing rate, its computing formula is as follows:
F L = C W R L / L , R L = 0.0882 ρ N / 1.7 - - - ( I I I )
Wherein, C wfor ground motion characteristic correction factor; R lfor triaxial strength ratio, N is the N value that standard penetration test (SPT) measures; σ vfor whole top pressure; σ v' be the effective pressure in top; The degree of depth of X distance ground surface; Thus the Degree of Liquefaction in region to be measured can be judged.
Preferably, described seismic event mechanical source is mechanical type impact origin or fluid pressure type impact origin.
Preferably, in described step one, lay multiple seismic receiver successively in region to be detected after, lay a mechanical source again; Distance between seismic receiver is not more than 1/4 of the seismic event wavelength that described mechanical source produces.
Preferably, described seismic receiver comprises for being embedded in underground and ground detection differently three-component seismometer string and hemisphere seismic receiver in the unidirectional string of seismic wave, land three-component seismometer string, hole.
Preferably, described data processing equipment utilizes the data of seismic receiver described in the process of SIMQKE_GR software, obtains seismic wave reaction spectrum file and saves as .srf formatted file.
Preferably, described multiple seismic receiver is that vertical direction connects bunchiness, and the string of described seismic receiver composition, circumferentially distributes for the center of circle with described seismic event mechanical source.
Preferably, the clutter in vibration information described in described filtering specifically refers to ground roll, direct wave, side echo, shear wave and scattering wave.
Preferably, the degree of depth of unidirectional wave detector embedding described in described seismic receiver in underground is 20-30m.
The present invention at least comprises following beneficial effect: method placement of mechanical focus and the seismic receiver in region to be measured, harbour utilizing seismic event to detect port area soil property Degree of Liquefaction of the present invention, corresponding reflection wave can be obtained quickly and easily, after utilizing SIMQKE_GR software to carry out processing the related data about reflection wave, the data such as more accurate seismic wave reaction spectrum can be obtained, for calculating below provides foundation.Applicant utilizes the correlation formula in step 5 of the present invention the anti-liquefied fraction calculating region to be measured, harbour or even pore water pressure ratio etc. and the harbour soil property of simple and fast more to liquefy relevant parameter, thus judge harbour soil property Degree of Liquefaction quickly and easily according to correlation parameter, for port area antidetonation construction provides foundation.
Part is embodied by explanation below by other advantage of the present invention, target and feature, part also will by research and practice of the present invention by those skilled in the art is understood.
Accompanying drawing explanation
The distribution schematic diagram that Fig. 1 is mechanical source described in the method utilizing seismic event to detect port area soil property Degree of Liquefaction of the present invention and seismic receiver;
Fig. 2 is the schematic diagram of seismic wave reaction spectrum described in the method utilizing seismic event to detect port area soil property Degree of Liquefaction of the present invention;
Fig. 3 is the curve map of different monitoring points cyclic shear stress ratio described in the method utilizing seismic event to detect port area soil property Degree of Liquefaction of the present invention.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in further detail, can implement according to this with reference to instructions word to make those skilled in the art.
Should be appreciated that used hereinly such as " to have ", other element one or more do not allotted in " comprising " and " comprising " term or the existence of its combination or interpolation.
As shown in Fig. 1 and 2 and 3, the invention provides the method utilizing seismic event to detect port area soil property Degree of Liquefaction, it comprises the following steps:
Step one, lays seismic event mechanical source 1 and seismic receiver 2 in region to be measured, harbour, and on described seismic receiver 2 access data harvester and data processing equipment;
Step 2, described data collector gathers the ground of seismic receiver 2 collection and the vibration information of underground, the clutter in vibration information described in described data processing equipment filtering, and computational analysis obtains seismic event wave spectrum intensity and seismic event longitudinal wave velocity;
Step 3, inputs described data processing equipment by basic data; Described basic data comprises: focal point location expression, soil property sampling description, stacking fold, seismic wave frequency spectrum analysis, wave field separation, one or more in rock & soil mechanical parameter; Such as, after utilizing SIMQKE_GR software process process related data, the seismic wave reaction produced as shown in Figure 2 is set a song to music line chart.
Step 4, utilize seismic event longitudinal wave velocity to calculate region to be measured soil property density, formula is as follows:
wherein VP is longitudinal wave velocity, and K is Lame's constant, and ∧ is shearing elasticity modulus;
Utilize seimic wave velocity spectral intensity value to calculate seismic wave reaction spectral intensity, computing formula is as follows:
wherein, Sv is velocity spectrum intensity level, and T is the cycle;
Calculate the cyclic shear stress ratio in region to be measured, computing formula is as follows:
Wherein,
α={0.0052(α max/SI)-0.0163}10 -2
b={0.0910(α max/SI)+0.0787}10 -2
Step 5, calculate liquefaction opposing rate, its computing formula is as follows:
F L = C W R L / L , R L = 0.0882 ρ N / 1.7 - - - ( I I I )
Wherein, C wfor ground motion characteristic correction factor; R lfor triaxial strength ratio, N is the N value that standard penetration test (SPT) measures; σ vfor whole top pressure; σ v' be the effective pressure in top; The degree of depth of X distance ground surface; Thus the Degree of Liquefaction in region to be measured can be judged.As shown in Figure 3, also there is very big difference in its cyclic shear stress ratio of monitoring point that the cyclic shear stress ratio in region to be measured is different.The harbour difference of identical soil property granularity or soil property density, its anti-Degree of Liquefaction is also different greatly, so, need to calculate the parameter values such as anti-Degree of Liquefaction further and carry out more accurately to react harbour soil property situation, for port construction provides abundanter foundation.
Wherein in an embodiment, as shown in Figure 1, described seismic event mechanical source 1 is mechanical type impact origin or fluid pressure type impact origin.
Wherein in an embodiment, as shown in Figure 1, in described step one, lay multiple seismic receiver 2 successively in region to be detected after, lay a mechanical source 1 again; Distance between seismic receiver 2 is not more than 1/4 of the seismic event wavelength that described mechanical source 1 produces.
Wherein in an embodiment, as shown in Figure 1, described seismic receiver 2 comprises for being embedded in underground and ground detection differently three-component seismometer string and hemisphere seismic receiver in the unidirectional string of seismic wave, land three-component seismometer string, hole.
Wherein in an embodiment, as shown in Figure 1, described data processing equipment utilizes the data of seismic receiver described in the process of SIMQKE_GR software, obtains seismic wave reaction spectrum file and saves as .srf formatted file.
Wherein in an embodiment, as shown in Figure 1, described multiple seismic receiver 2 connects bunchiness in vertical direction, and the string of described seismic receiver composition, circumferentially distributes for the center of circle with described seismic event mechanical source 1.
Wherein in an embodiment, the clutter in vibration information described in described filtering specifically refers to ground roll, direct wave, side echo, shear wave and scattering wave.
Wherein in an embodiment, as shown in Figure 1, described in described seismic receiver 2, the degree of depth of unidirectional wave detector embedding in underground is 20-30m.
Although embodiment of the present invention are open as above, but it is not restricted to listed in instructions and embodiment utilization, it can be applied to various applicable the field of the invention completely, for those skilled in the art, can easily realize other amendment, therefore do not deviating under the universal that claim and equivalency range limit, the present invention is not limited to specific details and illustrates here and the legend described.

Claims (8)

1. utilize seismic event to detect a method for port area soil property Degree of Liquefaction, it is characterized in that, comprise the following steps:
Step one, lays seismic event mechanical source and seismic receiver in region to be measured, harbour, and on described seismic receiver access data harvester and data processing equipment;
Step 2, described data collector gathers the ground of seismic receiver collection and the vibration information of underground, the clutter in vibration information described in described data processing equipment filtering, and computational analysis obtains seismic event wave spectrum intensity and seismic event longitudinal wave velocity;
Step 3, inputs described data processing equipment by basic data; Described basic data comprises: focal point location expression, soil property sampling description, stacking fold, seismic wave frequency spectrum analysis, wave field separation, one or more in rock & soil mechanical parameter;
Step 4, utilize seismic event longitudinal wave velocity to calculate region to be measured soil property density, formula is as follows:
wherein V pfor longitudinal wave velocity, K is Lame's constant, and ∧ is shearing elasticity modulus;
Utilize seimic wave velocity spectral intensity value to calculate seismic wave reaction spectral intensity, computing formula is as follows:
wherein, Sv is velocity spectrum intensity level, and T is the cycle;
Calculate the cyclic shear stress ratio in region to be measured, computing formula is as follows:
L=(b-ax)SIσ v/σ′ v(II)
Wherein,
a={0.0052(α max/SI)-0.0163}10 -2
b={0.0910(α max/SI)+0.0787}10 -2
Step 5, calculate liquefaction opposing rate, its computing formula is as follows:
F L = C W R L / L , R L = 0.0882 ρ N / 1.7 - - - ( I I I )
Wherein, C wfor ground motion characteristic correction factor; R lfor triaxial strength ratio, N is the N value that standard penetration test (SPT) measures; σ vfor whole top pressure; σ v 'for the effective pressure in top; The degree of depth of X distance ground surface; Thus the Degree of Liquefaction in region to be measured can be judged.
2. utilize seismic event to detect the method for port area soil property Degree of Liquefaction as claimed in claim 1, it is characterized in that, described seismic event mechanical source is mechanical type impact origin or fluid pressure type impact origin.
3. utilize seismic event to detect the method for port area soil property Degree of Liquefaction as claimed in claim 2, it is characterized in that, in described step one, lay multiple seismic receiver successively in region to be detected after, lay a mechanical source again; Distance between seismic receiver is not more than 1/4 of the seismic event wavelength that described mechanical source produces.
4. utilize seismic event to detect the method for port area soil property Degree of Liquefaction as claimed in claim 1, it is characterized in that, described seismic receiver comprises for being embedded in underground and ground detection differently three-component seismometer string and hemisphere seismic receiver in the unidirectional string of seismic wave, land three-component seismometer string, hole.
5. utilize seismic event to detect the method for port area soil property Degree of Liquefaction as claimed in claim 1, it is characterized in that, described data processing equipment utilizes the data of seismic receiver described in the process of SIMQKE_GR software, obtains seismic wave reaction spectrum file and saves as .srf formatted file.
6. utilize seismic event to detect the method for port area soil property Degree of Liquefaction as claimed in claim 4, it is characterized in that, described multiple seismic receiver is that vertical direction connects bunchiness, and the string of described seismic receiver composition, circumferentially distributes for the center of circle with described seismic event mechanical source.
7. utilize seismic event to detect the method for port area soil property Degree of Liquefaction as claimed in claim 1, it is characterized in that, the clutter in vibration information described in described filtering specifically refers to ground roll, direct wave, side echo, shear wave and scattering wave.
8. utilize seismic event to detect the method for port area soil property Degree of Liquefaction as claimed in claim 1, it is characterized in that, the degree of depth of unidirectional wave detector embedding described in described seismic receiver in underground is 20-30m.
CN201510918325.2A 2015-12-10 2015-12-10 A kind of method that seismic wave detects port area soil property Degree of Liquefaction Active CN105372701B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510918325.2A CN105372701B (en) 2015-12-10 2015-12-10 A kind of method that seismic wave detects port area soil property Degree of Liquefaction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510918325.2A CN105372701B (en) 2015-12-10 2015-12-10 A kind of method that seismic wave detects port area soil property Degree of Liquefaction

Publications (2)

Publication Number Publication Date
CN105372701A true CN105372701A (en) 2016-03-02
CN105372701B CN105372701B (en) 2017-09-12

Family

ID=55375049

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510918325.2A Active CN105372701B (en) 2015-12-10 2015-12-10 A kind of method that seismic wave detects port area soil property Degree of Liquefaction

Country Status (1)

Country Link
CN (1) CN105372701B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106324682A (en) * 2016-09-19 2017-01-11 中国石油化工股份有限公司 Surface structure investigation method applied to permafrost regions
CN106768847A (en) * 2017-02-14 2017-05-31 华侨大学 The analogy method and device of a kind of Potential of Seabed Under Wave Loading liquefaction water-jet depth

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6308135B1 (en) * 1999-10-07 2001-10-23 Golder Sierra Llc Soil liquefaction prevention by electro-osmosis during an earthquake event
CN104989411A (en) * 2015-06-16 2015-10-21 河海大学 Liquefaction discriminating device and method of pipe piece surrounding saturated sand caused by shield construction

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6308135B1 (en) * 1999-10-07 2001-10-23 Golder Sierra Llc Soil liquefaction prevention by electro-osmosis during an earthquake event
CN104989411A (en) * 2015-06-16 2015-10-21 河海大学 Liquefaction discriminating device and method of pipe piece surrounding saturated sand caused by shield construction

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YOUD T L: "Liquefaction resistance of soils: summery report from the NCEER and 1998 NCEER/NSF workshops on evaluation of liquefaction resistance of soils", 《JOURNAL OF GEOTECHNICAL AND GEOENVIRONMENTAL ENGINEERING》 *
张宗亮: "《200m级以上高心墙堆石坝关键技术研究及工程应用》", 31 December 2011 *
索丽生: "《水工设计手册》", 31 August 1984 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106324682A (en) * 2016-09-19 2017-01-11 中国石油化工股份有限公司 Surface structure investigation method applied to permafrost regions
CN106324682B (en) * 2016-09-19 2019-02-26 中国石油化工股份有限公司 A kind of surface structure survey method applied to frozen soil layer area
CN106768847A (en) * 2017-02-14 2017-05-31 华侨大学 The analogy method and device of a kind of Potential of Seabed Under Wave Loading liquefaction water-jet depth
CN106768847B (en) * 2017-02-14 2018-11-27 华侨大学 A kind of analogy method and device of Potential of Seabed Under Wave Loading liquefaction water-jet depth

Also Published As

Publication number Publication date
CN105372701B (en) 2017-09-12

Similar Documents

Publication Publication Date Title
Gao et al. Scattering of plane and cylindrical SH waves by a horseshoe shaped cavity
CN102445307A (en) Method for measuring flow rate and flow direction of single-well underground water and leaking point of reservoir, and measuring device thereof
CN104007462B (en) Crack prediction method based on attenuation anisotropy
CN105041331B (en) A kind of tunnel portal section anti-seismic fortified length computational methods
CN107092030A (en) Pass through the method for Rayleigh wave survey soft clay bag
Preisig et al. Modelling discharge rates and ground settlement induced by tunnel excavation
Walter et al. Mapping Rainfall-Triggered Slidequakes and Seismic Landslide-Volume Estimation at Heumoes SlopeAll rights reserved. No part of this periodical may be reproduced or transmitted in any form or by any means, electronic or mechanical, including photocopying, recording, or any information storage and retrieval system, without permission in writing from the publisher
CN112965136A (en) Multi-stage advanced detection method for water-rich karst tunnel
Zhang et al. Evolution trend of the Huangyanwo rock mass under the action of reservoir water fluctuation
Occhiena et al. A microseismic-based procedure for the detection of rock slope instabilities
Pirone et al. Analysis of Safety Factor in unsaturated pyroclastic slope
CN105700017B (en) A kind of method and apparatus of definite petroleum distribution data
CN105626056A (en) Method and device for determining three-dimensional stratum pore pressure of target area
CN105372701A (en) Method for detecting soil liquefaction degree of port by using seismic waves
Shen et al. Seismic response of cracking features in Jubao Mountain during the aftershocks of Jiuzhaigou Ms7. 0 earthquake
CN110377935A (en) Prediction technique of the pile peripheral earth excess pore water pressure along change in depth in pile driving process
CN106192865B (en) A kind of mud-rock flow Blaps femaralis computational methods and its application
Austin et al. Alongshore fluid motions in the swash zone of a sandy and gravel beach
Zaalishvili et al. Expected seismic intensity assessment taking into account local topography site effect
Sungkono et al. Assessment of Sidoarjo mud flow embankment stability using very low frequency electromagnetic method
Su et al. Experimental observation and geophysical modeling of the hydrological effects on Pleistocene glaciation deposits susceptible to landslide
Guo et al. Monitoring performance of slopes via ambient seismic noise recordings: Case study in a colluvium deposit
Raines et al. The application of passive seismic techniques to the detection of buried hollows
Hossain et al. Shear wave velocity as a function of standard penetration number and depth in Dhaka City, Bangladesh
Rodríguez-Castellanos et al. Near shore seismic movements induced by seaquakes using the boundary element method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Sun Huifeng

Inventor after: Liu Shijing

Inventor after: Wang Rong

Inventor after: Chen Jian

Inventor after: Dong Jie

Inventor before: Liu Shijing

Inventor before: Wang Rong

Inventor before: Chen Jian

Inventor before: Dong Jie

GR01 Patent grant
GR01 Patent grant