CN103278855B - Method for eliminating influence of roadways and terrains on apparent resistivity in direct-current exploration - Google Patents
Method for eliminating influence of roadways and terrains on apparent resistivity in direct-current exploration Download PDFInfo
- Publication number
- CN103278855B CN103278855B CN201310173286.9A CN201310173286A CN103278855B CN 103278855 B CN103278855 B CN 103278855B CN 201310173286 A CN201310173286 A CN 201310173286A CN 103278855 B CN103278855 B CN 103278855B
- Authority
- CN
- China
- Prior art keywords
- prime
- formula
- resistivity
- current
- apparent resistivity
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000004088 simulation Methods 0.000 claims abstract description 6
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 4
- 238000012876 topography Methods 0.000 claims description 15
- 238000005259 measurement Methods 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 5
- 239000004744 fabric Substances 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 3
- 238000011835 investigation Methods 0.000 claims description 3
- 239000000523 sample Substances 0.000 claims description 3
- 235000013350 formula milk Nutrition 0.000 description 26
- 230000008034 disappearance Effects 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000004615 ingredient Substances 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010291 electrical method Methods 0.000 description 1
- 230000005294 ferromagnetic effect Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method for eliminating the influence of roadways and terrains on the apparent resistivity in direct-current exploration. The method comprises the following steps of establishing a model according to a roadway shape, the roadway cutting depth, topographic relief elevation and the like; taking the minimum value of the formation resistivity of an exploration area as the resistivity of a uniform half-space containing the roadways and the terrains; putting a power supply electrode at a corresponding position in the model according to the type of the adopted device for direct-current exploration and the electrode distribution mode; performing one-time three-dimensional numerical simulation every time the power supply electrode moves, obtaining potential response of a background field, and interpolating to obtain potential difference between receiving electrodes as potential difference of the background field; deducing an apparent resistivity formula from a potential formula of a point current source in the horizontal uniform self-space; and subtracting the potential difference of the background field from the measured potential difference to obtain potential different without the influence of the roadways and the terrains, and substituting the potential difference into an apparent resistivity formula of a current source in an underground self-space to obtain the relative apparent resistivity. The method can be widely applied to direct-current resistivity exploration projects.
Description
Technical field
The invention belongs to electricity and electromagnetic prospecting field, be specifically related to eliminate the technology that in DC resistivity exploration, tunnel and landform affect apparent resistivity.
Background technology
Apparent resistivity is one of important method of resistivity prospecting data interpretation.In underground and direct current of ground electrical survey (-ing), the apparent resistivity by tunnel and the influence of topography distorts, and often mixes with the exception of geologic objective body, causes interpretation errors even to be judged by accident.Tunnel repair distorts to the one of constant current field normal distribution, and with device form, cloth pole mode, topographic relief degree and tunnel geometric configuration, the many factors such as stratum electric conductivity are relevant.This impact is nonlinear with the change of transmitting electrode, attempts to find factor of influence
[1]it is difficult for realizing correcting.The influence of topography is also distort to the one of constant current field.Common ratiometric correction method
[2], because the homogeneous half space resistivity value of simulating the response of pure landform is not easily determined, when underground electrical structure becomes complexity, the precision of this correcting method reduces, and even draws error result.If rely on the three-dimensional numerical value inverting comprising tunnel and landform completely
[3], ingredient necessary during this Data Inversion of disappearance apparent resistivity is explained.And the ill-posedness of Geophysical Inverse Problem, the problems such as the instability in 3-d inversion and multi-solution, calculated amount are huge need to solve.For this reason, the present invention is according to the sumproperties of current field in linear medium, utilize three-dimensional numerical value just drilling the impact of peeling off tunnel, landform, from the homogeneous half space current potential of point electric current source in underground, provide underground and direct current of ground explores the blanket apparent resistivity formula of various device.
There is following defect in prior art:
1., in the DC electrical method detection on underground and ground, the apparent resistivity by tunnel and the influence of topography distorts, and often mixes with the exception of geologic objective body, causes interpretation errors even to be judged by accident;
2. tunnel repair and device form, cloth pole mode, topographic relief degree and tunnel geometric configuration, the many factors such as stratum electric conductivity are relevant, and non-linearly change with polar distance of supplying power, and it is difficult for attempting to find factor of influence to realize correcting;
3. the method for the common ratiometric correction influence of topography, because the homogeneous half space resistivity value of simulating the response of pure landform is not easily determined, when underground electrical structure becomes complexity, the precision of this correcting method reduces, and even draws error result;
If 4. rely on the three-dimensional numerical value inverting comprising tunnel and landform completely, by ingredient necessary in this data interpretation of disappearance apparent resistivity.Moreover 3-d inversion self also also exists the problem such as instability and height nonuniqueness;
5. original underground and direct current of ground apparent resistivity, respectively from homogeneous half space current source current potential formula on the surface, and in the even total space, the current potential formula of current source is derived, can not be general.For overhead nearer subsurface investigation, the apparent resistivity that the current potential formula of original current source from the even total space is derived has larger error.
Documents
[1] Yue Jianhua, Li Zhidan. the tunnel repair in mine direct current method exploration. coal journal, 1999,24 (1): 7-10
[2] Holcomble HT, Jiracek G R. Three dimensional terrain correction in resistivitysurveys. Geophysics, 1984, 49(4): 439~452
[3] Wu little Ping. resistivity inversion under non-equilibrium drilling condition. Chinese Journal of Geophysics, 2005,48 (4): 932-936.
Summary of the invention
The object of the present invention is to provide a kind of method eliminated tunnel and landform and direct-current exploration apparent resistivity is affected, on the impact of apparent resistivity, improve the Explanation Accuracy of D.C. method to eliminate tunnel and landform in DC resistivity exploration.
In order to solve above technical matters, the present invention by the following technical solutions.
Eliminate the method that tunnel and landform affect direct-current exploration apparent resistivity, it is characterized in that the additivity according to constant current field in linear medium, the tunnel in homogeneous half space and landform are separated from ground electricity structure, specifically comprises the following steps:
Step one, according to roadway shape and the Modling model such as cutting depth, topographic relief elevation;
Step 2, using the resistivity of the minimum value of exploration area formation resistivity as the homogeneous half space containing tunnel and band landform;
Step 3, by transmitting electrode
with
the type of device adopted according to direct-current exploration and cloth pole mode, place correspondence position in a model;
Step 4, transmitting electrode is often mobile once carries out a Three-dimensional simulation, obtain the potential response of ambient field, then interpolation obtains receiving electrode
with
between potential difference (PD)
the potential difference (PD) of field as a setting;
Step 5, is arranged in horizontal homogeneous semispace current potential formula from point electric current source, derivation apparent resistivity formula: transmitting electrode in ground
with
form current return, have any point in ground
current potential formula
(1)
In formula
be ground resistivity,
supply current,
with
power supply respectively
and virtual source
arrive
distance,
with
power supply respectively
and virtual source
arrive
distance, any point in ground
current potential be
(2)
Formula (1) and formula (2) subtract each other, and obtain receiving electrode
with
between potential difference (PD)
(3)
Formula (3) is derived
as apparent resistivity, computing formula is as follows
(4)
When in formula (4)
for actual measurement potential difference (PD), the earth be non-horizontal homogeneous half space time, the ground resistivity in formula
be apparent resistivity, above formula is the formula calculating apparent resistivity, wherein
(5)
For electrode coefficient;
Step 6, by described actual measurement potential difference (PD)
deduct the potential difference (PD) of field as a setting
, obtain the potential difference (PD) without tunnel and the influence of topography
, utilize following formula
(6)
Calculate restivity
, be namely a cancellation the apparent resistivity of tunnel and the influence of topography.
Described restivity is added described homogeneous half space resistivity has been eliminated the absolute apparent resistivity of tunnel and the influence of topography.
The described step 5 apparent resistivity that current source current potential formula is derived from homogeneous half space, can be applicable to subsurface investigation and ground location simultaneously.
Described method can be applicable to whole electrical resistivity survey methods on underground and ground, comprise section and depth measurement, forward probe, roof and floor and lateral wall detection etc., and one pole-dipole, dipole-monopole, one pole-one pole, dipole-dipole etc. four kinds basic device type and combination in any.
The present invention has beneficial effect.The rock forming mineral of the overwhelming majority except ferromagnetic mineral, in the field strength range of D.C. method, is linear medium.Utilize the sumproperties of constant current field, eliminate tunnel and the influence of topography by partition method.This separation null method, solves the difficulty correcting and determine factor of influence in tunnel repair, does not have common ratioing technigue topographic correction homogeneous half space resistivity to determine the improper error brought.In the apparent resistivity of inversion interpretation important component part, applying three-dimensional numerical value is just drilled, and avoids and relies on the problems such as 3-d inversion instability, multi-solution and the calculated amount brought be huge completely.The apparent resistivity that current source current potential formula is derived from homogeneous half space, can switching naturally in underground and ground location, not only need not transformation for mula, and ensure that the precision of the apparent resistivity when tunnel is near ground.Whole electrical resistivity survey methods on the present invention compatibly lower and ground, as section and depth measurement, forward probe, roof and floor and lateral wall detection etc., and one pole-dipole, dipole-monopole, one pole-one pole, dipole-dipole etc. four kinds basic device type and combination in any.
Accompanying drawing explanation
Fig. 1 is underground point current field distribution plan, wherein
be point electric current source in ground,
the virtual source of symmetric position,
source and the virtual source distance to ground,
be any point in ground, the dotted line of band arrow is current line.
Fig. 2 is the rugged topography homogeneous half space computation model schematic diagram containing tunnel, wherein
be the earth with Air Interface,
be tunnel and Air Interface,
underground cutoff boundary,
emitting electrode, wherein
be power supply positive electrode,
power cathode,
receiving electrode,
source point to frontier point radius vector,
the outer normal vector on border, arrow indication device
moving direction.
Embodiment
Below in conjunction with accompanying drawing, technical scheme of the present invention is described in further detail.
For symmetrical quadrupole section detection in tunnel, mountain region.
According to the step one in technical scheme, according to roadway shape and the Modling model such as cutting depth, topographic relief elevation, as shown in Figure 2;
According to the step 2 in technical scheme, the minimum value of Qu Ce district formation resistivity
as the homogeneous half space resistivity containing tunnel and band landform;
According to the step 3 in technical scheme, by the transmitting electrode of symmetrical quadrupole device
with
place correspondence position in a model, as shown in Figure 2;
According to the step 4 in technical scheme, transmitting electrode often moves once, carries out a Three-dimensional simulation.If adopt the field domain Meta algorithm such as finite element, there is following boundary condition:
At the interface of earth's surface and air
on
(7)
At the interface of tunnel and air
on
(8)
At the cutoff boundary of underground infinite space
on, the first kind, Equations of The Second Kind, third boundary condition can be set.Wherein third boundary condition is general form
(9)
Above in formula,
current potential,
the distance that initial point is shown up a little,
the radial vector of source point to border
with outer normal vector
between angle.If tunnel is buried comparatively dark, ground-air border can be moved down, become border in ground, adopt the boundary condition of formula (9).Boundary condition formula (7) and (8) are accurate, and the boundary condition (9) of underground boundary is approximate.If possible, coboundary extends to ground as far as possible, obtains the Three-dimensional simulation result of degree of precision.
After utilizing Three-dimensional simulation to obtain the potential response of ambient field, then interpolation obtains receiving electrode
with
between potential difference (PD)
the potential difference (PD) of field as a setting;
According to formula (5) the calculation element coefficient of step 5 in technical scheme, the peak of the liftoff shape of roadway
if,
, so
(10)
Formula (10) is substituted into formula (5), and electrode coefficient is
(11)
According to the step 6 in technical scheme, potential difference (PD) will be surveyed
deduct the potential difference (PD) of field as a setting
, obtain the potential difference (PD) without tunnel and the influence of topography
, utilize formula (6)
(12)
Obtain the restivity without tunnel and the influence of topography;
Described in technical scheme, by restivity
add upper half-space resistivity
, the absolute apparent resistivity of be eliminated tunnel and the influence of topography.
Claims (4)
1. eliminate the method that tunnel and landform affect direct-current exploration apparent resistivity, it is characterized in that the additivity according to constant current field in linear medium, the tunnel in homogeneous half space and landform are separated from ground electricity structure, specifically comprises the following steps:
Step one, according to roadway shape and cutting depth, topographic relief elevation Modling model;
Step 2, using the resistivity of the minimum value of exploration area formation resistivity as the homogeneous half space containing tunnel and band landform;
Step 3, type of device electrodes of A and B adopted according to direct-current exploration and cloth pole mode, place correspondence position in a model;
Step 4, transmitting electrode is often mobile once carries out a Three-dimensional simulation, obtain the potential response of ambient field, then interpolation obtains the potential difference (PD) Δ U between receiving electrode M and N "
mNthe potential difference (PD) of field as a setting;
Step 5, is arranged in horizontal homogeneous semispace current potential formula from point electric current source, derivation apparent resistivity formula: in ground, electrodes of A and B form current return, have the current potential formula of any point M in ground
In formula, ρ is ground resistivity, I is supply current, AM and A ' M is the distance of power supply A and virtual source A ' to M respectively, BM and B ' M is the distance of power supply B and virtual source B ' to M respectively, and in ground, the current potential of any point N is
Formula (1) and formula (2) subtract each other, and obtain the potential difference (PD) Δ U between receiving electrode M and N
mN
The ρ that formula (3) is derived is as apparent resistivity, and computing formula is as follows
As Δ U in formula (4)
mNfor actual measurement potential difference (PD), the earth be non-horizontal homogeneous half space time, the ground resistivity ρ in formula is apparent resistivity, above formula be calculating apparent resistivity formula, wherein
For electrode coefficient;
Step 6, by described actual measurement potential difference (PD) Δ U
mNdeduct the potential difference (PD) Δ U of field as a setting "
mN, obtain the potential difference (PD) Δ U ' without tunnel and the influence of topography
mN, utilize following formula
Calculate restivity ρ ' for eliminating the apparent resistivity in tunnel and the influence of topography.
2. the as claimed in claim 1 method eliminated tunnel and landform and affect direct-current exploration apparent resistivity, is characterized in that, described restivity is added described homogeneous half space resistivity has been eliminated the absolute apparent resistivity of tunnel and the influence of topography.
3. the as claimed in claim 1 method eliminated tunnel and landform and affect direct-current exploration apparent resistivity, it is characterized in that, the apparent resistivity of described step 5 current source current potential formula derivation from homogeneous half space, can be applicable to subsurface investigation and ground location simultaneously.
4. the as claimed in claim 1 method eliminated tunnel and landform and direct-current exploration apparent resistivity is affected, it is characterized in that, described method can be applicable to whole electrical resistivity survey on underground and ground, comprises section and depth measurement, forward probe, roof and floor and lateral wall detection.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310173286.9A CN103278855B (en) | 2013-05-13 | 2013-05-13 | Method for eliminating influence of roadways and terrains on apparent resistivity in direct-current exploration |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310173286.9A CN103278855B (en) | 2013-05-13 | 2013-05-13 | Method for eliminating influence of roadways and terrains on apparent resistivity in direct-current exploration |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103278855A CN103278855A (en) | 2013-09-04 |
CN103278855B true CN103278855B (en) | 2015-06-24 |
Family
ID=49061426
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310173286.9A Expired - Fee Related CN103278855B (en) | 2013-05-13 | 2013-05-13 | Method for eliminating influence of roadways and terrains on apparent resistivity in direct-current exploration |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103278855B (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105116453B (en) * | 2015-08-14 | 2017-10-17 | 中国石油天然气股份有限公司 | Transient electromagnetic exploration method and device for natural gas hydrate in permafrost zone |
CN107272068A (en) * | 2017-06-09 | 2017-10-20 | 北京市地质调查研究院 | It is a kind of to utilize the method without anti-interference across the hole resistivity method Underground space of cable |
CN108535783B (en) * | 2018-04-16 | 2019-07-05 | 中煤科工集团西安研究院有限公司 | A kind of pair of resistivity section carries out the method and device of layer bit correction |
CN109026021B (en) * | 2018-08-13 | 2024-02-27 | 中国电子科技集团公司第二十二研究所 | Focusing system and method of shield tunneling machine |
CN110297016A (en) * | 2019-06-19 | 2019-10-01 | 安徽国科骄辉科技有限公司 | The calculation method and imaging device of the apparent resistivity of cylinder model |
CN110703346B (en) * | 2019-11-14 | 2021-06-11 | 张道清 | Same-frequency and same-source ratio measurement method for natural electric field |
CN111273355B (en) * | 2020-03-16 | 2021-08-13 | 中国科学技术大学 | Advanced detection method and advanced detection system for roadway |
CN113885086B (en) * | 2021-08-05 | 2023-06-13 | 中煤科工集团西安研究院有限公司 | Underground direct-current equatorial direction dipole dynamic source abnormity self-explicit advanced detection method |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2402217A1 (en) * | 1977-09-01 | 1979-03-30 | Geophysique Cie Gle | ELECTRICAL PROSPECTION OF THE BASEMENT, WITH LINEIC ELECTRODES |
US5809458A (en) * | 1996-09-05 | 1998-09-15 | Western Atlas International, Inc. | Method of simulating the response of a through-casing electrical resistivity well logging instrument and its application to determining resistivity of earth formations |
US6344746B1 (en) * | 1999-12-03 | 2002-02-05 | Baker Hughes Incorporated | Method for processing the lapse measurements |
CN101556340A (en) * | 2008-04-10 | 2009-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | Three-dimensional small surface element magnetotelluric continuous array data acquisition method |
CN102621585A (en) * | 2012-03-23 | 2012-08-01 | 中煤科工集团西安研究院 | Method for rapidly and accurately correcting frame influence for transient electromagnetic data |
-
2013
- 2013-05-13 CN CN201310173286.9A patent/CN103278855B/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2402217A1 (en) * | 1977-09-01 | 1979-03-30 | Geophysique Cie Gle | ELECTRICAL PROSPECTION OF THE BASEMENT, WITH LINEIC ELECTRODES |
US5809458A (en) * | 1996-09-05 | 1998-09-15 | Western Atlas International, Inc. | Method of simulating the response of a through-casing electrical resistivity well logging instrument and its application to determining resistivity of earth formations |
US6344746B1 (en) * | 1999-12-03 | 2002-02-05 | Baker Hughes Incorporated | Method for processing the lapse measurements |
CN101556340A (en) * | 2008-04-10 | 2009-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | Three-dimensional small surface element magnetotelluric continuous array data acquisition method |
CN102621585A (en) * | 2012-03-23 | 2012-08-01 | 中煤科工集团西安研究院 | Method for rapidly and accurately correcting frame influence for transient electromagnetic data |
Non-Patent Citations (2)
Title |
---|
山区电磁测深的地形影响分析;刘宏等;《石油物探》;20050830;第44卷(第04期);404-407 * |
高密度电阻率法地形影响校正;简兴祥等;《物探化探计算技术》;20080715;第30卷(第04期);303-305 * |
Also Published As
Publication number | Publication date |
---|---|
CN103278855A (en) | 2013-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103278855B (en) | Method for eliminating influence of roadways and terrains on apparent resistivity in direct-current exploration | |
Martínez-Moreno et al. | Combined microgravity, electrical resistivity tomography and induced polarization to detect deeply buried caves: Algaidilla cave (Southern Spain) | |
CA2921817C (en) | Mapping resistivity distribution within the earth | |
McGrath et al. | Integrated high-resolution geophysical investigations as potential tools for water resource investigations in karst terrain | |
CN103064124B (en) | A kind of ratio approach correcting electromagnetic survey Considering Terrain Effect | |
Boucher et al. | Using 2D inversion of magnetic resonance soundings to locate a water-filled karst conduit | |
CN103809204B (en) | A kind of collecting method of field audio magnetotelluric method | |
CN107085240B (en) | Slope magnetofluid detection system and method | |
Cardimona | zywvutsrqponmlihgfed | |
Ramazi et al. | Contribution of geophysical inversion theory and geostatistical simulation to determine geoelectrical anomalies | |
Berge | Electrical resistivity tomography investigations on a paleoseismological trenching study | |
CN108873073B (en) | A kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency | |
Dennis et al. | Transient electromagnetic surveys for the measurement of near-surface electrical anisotropy | |
Hennig et al. | The effect of dike geometry on different resistivity configurations | |
Aning et al. | Electrical resistivity as a geophysical mapping tool; a case study of the new art department, Knust-Ghana | |
CA2834518A1 (en) | Method and system for surveying a distribution of chargeability in a volume of earth | |
Neyamadpour et al. | Use of four-electrode arrays in three-dimensional electrical resistivity imaging survey | |
Oskooi et al. | Resolution of airborne VLF data | |
Smith et al. | A comparison of airborne electromagnetic data with ground resistivity data over the Midwest deposit in the Athabasca basin | |
Mustopa et al. | Subsurface structure of Tangkuban Parahu area derived from CSAMT and gravity investigation | |
Amini et al. | Anomaly enhancement in 2D electrical resistivity imaging method using a residual resistivity technique | |
CN104133254A (en) | Electric and seismic comprehensive detection instrument for underground coal mine | |
Qi et al. | An array source TEM apparent resistivity conversion and differential conductivity imaging method and its potential applications | |
Pethö | Geological applications of the VLF method | |
Chen et al. | Three-dimensional observation system of direct current method based on overlapping electrodes and its application |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CB03 | Change of inventor or designer information |
Inventor after: Yan Shu Inventor after: Wu Qiao Inventor after: Zheng Ji Inventor after: Xue Guoqiang Inventor after: Qiu Weizhong Inventor after: Di Qingyun Inventor after: Li Gang Inventor before: Yan Shu Inventor before: Wu Qiao Inventor before: Zheng Ji |
|
COR | Change of bibliographic data | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150624 |
|
CF01 | Termination of patent right due to non-payment of annual fee |