CN102788845A - Barker coding excitation ultrasonic detection method of concrete structure defect - Google Patents

Barker coding excitation ultrasonic detection method of concrete structure defect Download PDF

Info

Publication number
CN102788845A
CN102788845A CN2012103230769A CN201210323076A CN102788845A CN 102788845 A CN102788845 A CN 102788845A CN 2012103230769 A CN2012103230769 A CN 2012103230769A CN 201210323076 A CN201210323076 A CN 201210323076A CN 102788845 A CN102788845 A CN 102788845A
Authority
CN
China
Prior art keywords
signal
signals
barker code
transducer
big
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN2012103230769A
Other languages
Chinese (zh)
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.)
Yellow River Institute of Hydraulic Research
Original Assignee
Yellow River Institute of Hydraulic Research
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 Yellow River Institute of Hydraulic Research filed Critical Yellow River Institute of Hydraulic Research
Priority to CN2012103230769A priority Critical patent/CN102788845A/en
Publication of CN102788845A publication Critical patent/CN102788845A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

The invention discloses a barker coding excitation ultrasonic detection method of a concrete structure defect. The method comprises 1 using a buck coding excitation energy converter at an emission end, enabling carrier waves to be cosine signals and enabling the energy converter to transmit width time signals; and 2 using a matched filtering method at a receiving end, conducting pulse compression on received reflected signals, and finally obtaining narrow pulse signals with high resolution ratio and big major and minor lobes. Barker codes are adopted at the emission end to modulate coding excitation and send modulated big-time-width signals, corresponding matched filtering processing is carried out on the received signals, and large-amplitude narrow-pulse signals can be obtained. The big-time-width signals are big in energy so that a certain range of reflected signals can be obtained when depth defects are detected, the major and minor lobes after the matched filtering processing are big, echo signals are high in range, defect positions can be accurately judged according to defect reflection time, and detection capability of the defects is improved accordingly.

Description

The Barker code excitation supersonic detection method of xoncrete structure defective
Technical field:
The present invention relates to a kind of xoncrete structure defect inspection method, relate in particular to a kind of Barker code excitation supersonic detection method that encourages the xoncrete structure defective of ultrasonic signal based on Barker code.
Background technology:
The means of at present ultrasound detection being carried out in the size and the position of concrete inner structure defective can be divided into two types substantially: the one, and pulse echo method at concrete surface excitation ultrasound signal, receives echoed signal in the same side and carries out the identification of defective; Another kind of is the sound wave through transmission technique, transmits in concrete one side, receive at opposite side, when receiving the sound of signal and changes in amplitude judge whether inner structure exists defective.Ultrasonic detection equipment is simpler, operate also very convenient, so be widely used in the xoncrete structure defects detection.
At the scene in the testing process, ground unrest disturbs greatlyyer sometimes, makes echo and noise signal aliasing, and perhaps echo is submerged in the noise signal fully, and the pulse excitation that places an order of this situation is difficult to judge flaw echoes, makes troubles to defects detection.On the other hand, for the detection of thicker concrete component, because sonic propagation length is bigger apart from the amplitude decay, be difficult for determining echoed signal in the reception signal, this also is one of problem of puzzlement engineering detecting personnel.
Summary of the invention:
, the echoed signal time low to xoncrete structure defects detection signal to noise ratio (S/N ratio) is difficult for interpretation; And it is long apart from concrete component detection problem; The purpose of this invention is to provide a kind of wide pulse signal through the emission macro-energy; In the receiving end matched filtering, obtain the echoed signal that main lobe is narrower, amplitude is higher, thus the Barker code excitation supersonic detection method of the xoncrete structure defective of the detectability of raising defective.
Technical scheme of the present invention realizes in the following manner:
A kind of Barker code excitation supersonic detection method of xoncrete structure defective: this method comprises:
1. adopt the code-excited transducer of Barker code at transmitting terminal, wherein carrier wave is the cosine impulse signal, makes the signal of transducer length launch time broad; 2. take the method for matched filtering at receiving end, the reflected signal that receives is carried out pulse compression, finally obtain the narrow pulse signal that resolution is higher, principal subsidiary lobe is bigger;
The practical implementation step:
Adopting the longest Barker code is N=13, {+1 ,+1 ,+1 ,+1 ,+1 ,-1 ,-1 ,+1 ,+1 ,-1 ,+1 ,-1 ,+1 },
1. confirm coding parameter according to defective approximate depth ; Barker code length
Figure 2012103230769100002DEST_PATH_IMAGE004
; Carrier frequency (transducer dominant frequency)
Figure 2012103230769100002DEST_PATH_IMAGE006
, unit chip carrier cycle number is n; The concrete velocity of wave is
Figure 2012103230769100002DEST_PATH_IMAGE008
,
Figure 2012103230769100002DEST_PATH_IMAGE010
; The duration of transmitted wave is less than the reflection interval that echo is propagated;
2. the Barker code coded signal is calculated in programming:
Figure 2012103230769100002DEST_PATH_IMAGE012
; Barker code is 1 o'clock
Figure 2012103230769100002DEST_PATH_IMAGE014
, and Barker code is-1 o'clock
Figure 2012103230769100002DEST_PATH_IMAGE016
.
3. coded signal
Figure 2012103230769100002DEST_PATH_IMAGE018
lead-in signal generator is used for excitation; Transmitting transducer directly docks with receiving transducer; Recorder signal ;
Figure 825787DEST_PATH_IMAGE020
is as the matched filtering function; Signal length is
Figure 2012103230769100002DEST_PATH_IMAGE022
, and SF remains unchanged in all testing processes;
4. to seized concrete component; Adopt reflectometry or transmission beam method to detect;
Figure 680610DEST_PATH_IMAGE006
is consistent with the transducer dominant frequency for transmission carrier frequency; Coded signal
Figure 460347DEST_PATH_IMAGE018
is through signal generator stimulated emission transducer; The recorder signal is , and signal length is
Figure 2012103230769100002DEST_PATH_IMAGE026
;
5. the echoed signal that receives
Figure 793240DEST_PATH_IMAGE024
is calculated with
Figure 548225DEST_PATH_IMAGE020
, obtained signal
Figure 2012103230769100002DEST_PATH_IMAGE028
;
Figure 2012103230769100002DEST_PATH_IMAGE032
; As if
Figure 2012103230769100002DEST_PATH_IMAGE034
, get
Figure 2012103230769100002DEST_PATH_IMAGE036
in the computation process. individual point of
Figure 2012103230769100002DEST_PATH_IMAGE038
representative
Figure 934469DEST_PATH_IMAGE020
waveform,
Figure 335156DEST_PATH_IMAGE038
is increased to
Figure 2012103230769100002DEST_PATH_IMAGE040
from 1.
N is a n point of waveform;
6 According to
Figure 941718DEST_PATH_IMAGE042
the peak time
Figure 2012103230769100002DEST_PATH_IMAGE044
, and amplitude determined position of the defect within the concrete.
Can also adopt short Barker code is N=2, N=3, N=4, N=5, N=7 or N=11.
Good effect of the present invention is:
Barker code is a kind of of pseudorandom two-phase sign indicating number, and the Barker code pumping signal has desirable ambiguity function, has good distance resolution.Adopt the excitation of Barker code modulating-coding at transmitting terminal, send bandwidth signals when ovennodulation big,, can obtain large amplitude, narrow pulse signal receiving signal through the corresponding matched Filtering Processing.Bandwidth signals is bigger because of energy when big; So when detecting dark defective, still can obtain the reflected signal of certain amplitude, bigger through the principal subsidiary lobe of matched filter processing, the echoed signal that amplitude is higher; Can accurately judge defective locations according to the defect reflection time, thereby improve the detectability of defective.
Description of drawings:
Fig. 1 is the principle schematic that concrete component of the present invention detects;
Peak value contrasted synoptic diagram after Fig. 2 received signal peak and the matched filtering of Barker code pumping signal for monopulse of the present invention encourages;
Fig. 3 is a concrete structure member crevices illustraton of model of the present invention;
Fig. 4 is that Barker code of the present invention detects result.
Embodiment:
Below in conjunction with accompanying drawing the present invention is done further explain:
Embodiment 1,
(1) detection in hole
(20cm * 20cm * 20cm), as shown in Figure 1, it is the horizontal cylindrical hole of 2cm that there is a diameter in test block central authorities to select a rock test block; Take both sides that the single-shot list of wearing is received detection mode; Measuring point is 10 and evenly distributes that wherein the 6th measuring point is just in time relative with the hole, and the transducer dominant frequency is 500kHz.At first detect with the monopulse energisation mode.Use carrier frequency to detect as 13 Barker code excitations of 500kHz then, main lobe peak value after the matched filtering and rect.p. detect the contrast of main lobe peak value and see Fig. 2.
Detect step:
1. Barker code length
Figure 2012103230769100002DEST_PATH_IMAGE046
; Carrier frequency (transducer dominant frequency)
Figure 2012103230769100002DEST_PATH_IMAGE048
500kHz, unit chip carrier cycle number is n=1; The concrete velocity of wave is
Figure 2012103230769100002DEST_PATH_IMAGE050
,
Figure 810448DEST_PATH_IMAGE010
.This formula shows: the duration of transmitted wave is less than the reflection interval that echo is propagated;
2. Barker code coding letter
Figure 2012103230769100002DEST_PATH_IMAGE052
is calculated in programming
Figure 502460DEST_PATH_IMAGE012
; Barker code is 1 o'clock
Figure 239472DEST_PATH_IMAGE014
, and Barker code is-1 o'clock
Figure 700540DEST_PATH_IMAGE016
.
3. with coded signal
Figure 130385DEST_PATH_IMAGE018
The lead-in signal generator is used for excitation, and transmitting transducer directly docks with receiving transducer, the recorder signal
Figure 309693DEST_PATH_IMAGE020
, As the matched filtering function, signal length does
Figure 228288DEST_PATH_IMAGE022
=3000, SF is 1 * 10 in all testing processes -8S remains unchanged;
4. to seized concrete component; Adopt reflectometry or transmission beam method to detect;
Figure 501137DEST_PATH_IMAGE006
is consistent with the transducer dominant frequency for transmission carrier frequency; Coded signal
Figure 230059DEST_PATH_IMAGE018
is through signal generator stimulated emission transducer; The recorder signal is
Figure 237767DEST_PATH_IMAGE024
, and signal length is
Figure 735745DEST_PATH_IMAGE026
=40000.
5. the echoed signal that receives
Figure 507392DEST_PATH_IMAGE024
is calculated with
Figure 395713DEST_PATH_IMAGE020
, obtained signal .
Figure DEST_PATH_IMAGE030A
Figure 506069DEST_PATH_IMAGE032
; As if , get
Figure 824235DEST_PATH_IMAGE036
in the computation process.
6 According to
Figure 244852DEST_PATH_IMAGE042
The peak time
Figure 655105DEST_PATH_IMAGE044
, and amplitude determined position of the defect within the concrete.
Can find out that by Fig. 2 the peak value of signal is minimum after the 6th the measuring point matched filtering of Barker code excitation, can judge the position in hole thus.And the 6th measuring point peak value of monopulse mode of excitation is not minimum.Analyze the former because defective build-up effect to the long pulse pumping signal, the excitation signal pulse is long more, and the energy loss that runs into behind the hole is big more, show the height of peak value after the matched filtering, and the hole is less to the energy affect of single pulse signal.
Embodiment 2,
Crack Detection:
(Fig. 3 is seen at the dark 12cm of sillar place in the crack for 24cm * 24cm * 24cm), the wide about 1.5mm in crack to select sillar that horizontal fracture is arranged in the middle of one.The detection mode of taking the single-shot list to receive, transmitting probe and receiving transducer are at the sillar top; Use carrier frequency to detect as 13 Barker code excitations of 500kHz, Fig. 4 is the testing result of Barker code excitation, and visible main lobe peak value is narrower, and peak amplitude is higher, also can obviously find out secondary and three echoed signals.Can find out that by Fig. 4 there is reflected signal at 55 μ s places, the position that can judge the crack by the concrete velocity of wave is about 12cm.

Claims (2)

1. the Barker code of an xoncrete structure defective encourages supersonic detection method: this method comprises:
1. adopt the code-excited transducer of Barker code at transmitting terminal, wherein carrier wave is a cosine signal, makes transducer launch signal when wide; 2. take the method for matched filtering at receiving end, the reflected signal that receives is carried out pulse compression, finally obtain the narrow pulse signal that resolution is higher, principal subsidiary lobe is bigger;
The practical implementation step:
Adopting the longest Barker code is N=13, {+1 ,+1 ,+1 ,+1 ,+1 ,-1 ,-1 ,+1 ,+1 ,-1 ,+1 ,-1 ,+1 },
(1). confirms coding parameter according to the defective approximate depth; Barker code length ; Carrier frequency is transducer dominant frequency
Figure 2012103230769100001DEST_PATH_IMAGE006
, and unit chip carrier cycle number is n; The concrete velocity of wave is
Figure 2012103230769100001DEST_PATH_IMAGE008
,
Figure 2012103230769100001DEST_PATH_IMAGE010
; The duration of transmitted wave is less than the reflection interval that echo is propagated;
(2). the Barker code coded signal is calculated in programming:
Figure 2012103230769100001DEST_PATH_IMAGE012
; Barker code is 1 o'clock
Figure 2012103230769100001DEST_PATH_IMAGE014
, and Barker code is-1 o'clock
Figure 2012103230769100001DEST_PATH_IMAGE016
;
(3). with coded signal
Figure 2012103230769100001DEST_PATH_IMAGE018
The lead-in signal generator is used for excitation, and transmitting transducer directly docks with receiving transducer, the recorder signal
Figure 2012103230769100001DEST_PATH_IMAGE020
,
Figure 835327DEST_PATH_IMAGE020
As the matched filtering function, signal length does
Figure 2012103230769100001DEST_PATH_IMAGE022
, SF remains unchanged in all testing processes;
(4). to seized concrete component; Adopt reflectometry or transmission beam method to detect;
Figure 822874DEST_PATH_IMAGE006
is consistent with the transducer dominant frequency for transmission carrier frequency; Coded signal is through signal generator stimulated emission transducer; The recorder signal is
Figure 2012103230769100001DEST_PATH_IMAGE024
, and signal length is
Figure 2012103230769100001DEST_PATH_IMAGE026
;
(5). the echoed signal that receives
Figure 233313DEST_PATH_IMAGE024
is calculated with
Figure 517664DEST_PATH_IMAGE020
, obtained signal
Figure 2012103230769100001DEST_PATH_IMAGE028
;
Figure 2012103230769100001DEST_PATH_IMAGE030
Figure 2012103230769100001DEST_PATH_IMAGE032
; In the computation process if
Figure 2012103230769100001DEST_PATH_IMAGE034
; Get
Figure 2012103230769100001DEST_PATH_IMAGE036
;
Figure 272703DEST_PATH_IMAGE038
individual point of
Figure 2012103230769100001DEST_PATH_IMAGE038
representative
Figure 244705DEST_PATH_IMAGE020
waveform,
Figure 516603DEST_PATH_IMAGE038
is increased to
Figure 2012103230769100001DEST_PATH_IMAGE040
from 1;
N is a n point of
Figure 2012103230769100001DEST_PATH_IMAGE042
waveform;
(6) According to
Figure 768593DEST_PATH_IMAGE042
the peak time and amplitude determine the location of defects inside the concrete.
2. the Barker code excitation supersonic detection method of xoncrete structure defective according to claim 1: it is characterized in that: Barker code is N=2, N=3, N=4, N=5, N=7 or N=11.
CN2012103230769A 2012-09-04 2012-09-04 Barker coding excitation ultrasonic detection method of concrete structure defect Pending CN102788845A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012103230769A CN102788845A (en) 2012-09-04 2012-09-04 Barker coding excitation ultrasonic detection method of concrete structure defect

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012103230769A CN102788845A (en) 2012-09-04 2012-09-04 Barker coding excitation ultrasonic detection method of concrete structure defect

Publications (1)

Publication Number Publication Date
CN102788845A true CN102788845A (en) 2012-11-21

Family

ID=47154305

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012103230769A Pending CN102788845A (en) 2012-09-04 2012-09-04 Barker coding excitation ultrasonic detection method of concrete structure defect

Country Status (1)

Country Link
CN (1) CN102788845A (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103149274A (en) * 2013-01-28 2013-06-12 中国科学院声学研究所 Defect detecting method of concrete
CN106525969A (en) * 2016-10-27 2017-03-22 中国电建集团贵阳勘测设计研究院有限公司 Device and method for carrying out nondestructive testing on anchor rod by adopting cosine linear scanning signal
CN107281687A (en) * 2017-08-22 2017-10-24 河海大学常州校区 A kind of fire hydrant aqueous condition detection method
CN109781850A (en) * 2019-01-26 2019-05-21 零声科技(苏州)有限公司 A kind of electromagnetic acoustic on-line monitoring system based on impulse compression method
CN110865124A (en) * 2019-11-27 2020-03-06 华东理工大学 Nonlinear ultrasonic guided wave detection system and method based on linear power amplifier
CN112986398A (en) * 2021-03-15 2021-06-18 南昌航空大学 Electromagnetic ultrasonic Lamb wave transducer and online detection system and method
CN112986399A (en) * 2021-03-15 2021-06-18 南昌航空大学 Electromagnetic ultrasonic SH (shear) guided wave transducer and online detection system and method
CN113433226A (en) * 2020-03-23 2021-09-24 中国地震局地球物理研究所 Method and equipment for detecting rock physical property based on digital coding signal

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101224115A (en) * 2008-01-28 2008-07-23 深圳市蓝韵实业有限公司 Electropult of ultrasonic diagnosis equipment
US20080229833A1 (en) * 2003-12-02 2008-09-25 Katsunori Asafusa Ultrasonic Diagnostic Apparatus
CN102393445A (en) * 2011-10-24 2012-03-28 沈阳建筑大学 Pipeline structure damage monitoring method based on piezoelectric ceramic sensor and guide wave analysis
CN102645267A (en) * 2012-05-02 2012-08-22 中国科学院声学研究所 Method and system for detecting elastic wave velocity of rock component

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080229833A1 (en) * 2003-12-02 2008-09-25 Katsunori Asafusa Ultrasonic Diagnostic Apparatus
CN101224115A (en) * 2008-01-28 2008-07-23 深圳市蓝韵实业有限公司 Electropult of ultrasonic diagnosis equipment
CN102393445A (en) * 2011-10-24 2012-03-28 沈阳建筑大学 Pipeline structure damage monitoring method based on piezoelectric ceramic sensor and guide wave analysis
CN102645267A (en) * 2012-05-02 2012-08-22 中国科学院声学研究所 Method and system for detecting elastic wave velocity of rock component

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
W. H. CHEN, ETC.: "Ultrasonic non-destructive testing using Barker code pulse compression techniques", 《ULTRASONICS》, vol. 26, no. 1, 31 January 1988 (1988-01-31), XP025835587, DOI: doi:10.1016/0041-624X(88)90044-3 *
张群逸: "雷达中的相位编码信号与处理", 《火控雷达技术》, vol. 34, 31 December 2005 (2005-12-31), pages 29 - 32 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103149274A (en) * 2013-01-28 2013-06-12 中国科学院声学研究所 Defect detecting method of concrete
CN103149274B (en) * 2013-01-28 2015-10-07 中国科学院声学研究所 A kind of concrete defect detection method
CN106525969A (en) * 2016-10-27 2017-03-22 中国电建集团贵阳勘测设计研究院有限公司 Device and method for carrying out nondestructive testing on anchor rod by adopting cosine linear scanning signal
CN107281687B (en) * 2017-08-22 2020-03-10 河海大学常州校区 Fire hydrant water condition detection method
CN107281687A (en) * 2017-08-22 2017-10-24 河海大学常州校区 A kind of fire hydrant aqueous condition detection method
CN109781850A (en) * 2019-01-26 2019-05-21 零声科技(苏州)有限公司 A kind of electromagnetic acoustic on-line monitoring system based on impulse compression method
CN109781850B (en) * 2019-01-26 2021-06-22 零声科技(苏州)有限公司 Electromagnetic ultrasonic online monitoring system based on pulse compression method
CN110865124A (en) * 2019-11-27 2020-03-06 华东理工大学 Nonlinear ultrasonic guided wave detection system and method based on linear power amplifier
CN110865124B (en) * 2019-11-27 2022-09-27 华东理工大学 Nonlinear ultrasonic guided wave detection system and method based on linear power amplifier
CN113433226A (en) * 2020-03-23 2021-09-24 中国地震局地球物理研究所 Method and equipment for detecting rock physical property based on digital coding signal
CN112986398A (en) * 2021-03-15 2021-06-18 南昌航空大学 Electromagnetic ultrasonic Lamb wave transducer and online detection system and method
CN112986399A (en) * 2021-03-15 2021-06-18 南昌航空大学 Electromagnetic ultrasonic SH (shear) guided wave transducer and online detection system and method
CN112986398B (en) * 2021-03-15 2022-06-28 南昌航空大学 Electromagnetic ultrasonic Lamb wave transducer and online detection system and method
CN112986399B (en) * 2021-03-15 2022-06-28 南昌航空大学 Electromagnetic ultrasonic SH guided wave transducer and online detection system and method

Similar Documents

Publication Publication Date Title
CN102788845A (en) Barker coding excitation ultrasonic detection method of concrete structure defect
CA2949040C (en) An apparatus and method for measuring the pressure inside a pipe or container
CN105066918A (en) Ultrasonic underwater target thickness measuring system and thickness measuring method
EP2799820A1 (en) Liquid surface level measurement device, method, and program
RU2010127782A (en) DEVICE AND METHOD FOR PIPELINE CONTROL USING ULTRASONIC WAVES OF TWO DIFFERENT TYPES
CN102565800B (en) Detection device and detecting method
CN109114436A (en) Line clogging detection system based on low-frequency sound wave
CN109991590B (en) System and method for testing low-frequency emission characteristic of transducer in pressure tank in limited space
CN106290580B (en) Vacuum high-low frequency acoustic measurement device and method
CN105953079A (en) Pipe cleaning ball positioning system for oil and gas pipeline
US10527590B2 (en) Apparatus and method for inspecting a pipeline
CN103149274B (en) A kind of concrete defect detection method
CN102383786A (en) Sector cement bonding and sound wave time difference logging instrument
EA018239B1 (en) Acoustic thickness measurements using gas as a coupling medium
CN104049038A (en) Ultrasonic-acoustic emission detection method for composite material
CN104062663B (en) A kind of multi-beam seabed subbottom profile probing equipment
CN104536003A (en) Ultrasonic distance measuring method and device based on multiple emission frequencies
CN104634873A (en) Ultrasonic detection system and method for damage of steel wires in bridge-cable anchoring area
CN102645267B (en) Method and system for detecting elastic wave velocity of rock component
JP2011047763A (en) Ultrasonic diagnostic device
CN105181047A (en) Coastal water flow rate high-frequency acoustic monitoring system and monitoring method thereof
CN104501909A (en) Liquid level measuring device and measuring method for small measuring range based on supersonic
CA2908682A1 (en) Conical ultrasonic probe
JP6024413B2 (en) measuring device
CN107462588A (en) A kind of trees Inner Defect Testing method using GPR as instrument

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20121121