CN108267722A - Geological radar echo-signal physical property, which deconstructs to know with detection target number reconstruct intelligence, takes method - Google Patents
Geological radar echo-signal physical property, which deconstructs to know with detection target number reconstruct intelligence, takes method Download PDFInfo
- Publication number
- CN108267722A CN108267722A CN201810064500.XA CN201810064500A CN108267722A CN 108267722 A CN108267722 A CN 108267722A CN 201810064500 A CN201810064500 A CN 201810064500A CN 108267722 A CN108267722 A CN 108267722A
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- spectrum energy
- spectrum
- radar
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000000704 physical effect Effects 0.000 title claims abstract description 35
- 238000001514 detection method Methods 0.000 title claims abstract description 34
- 238000001228 spectrum Methods 0.000 claims abstract description 75
- 230000004044 response Effects 0.000 claims abstract description 46
- 238000004458 analytical method Methods 0.000 claims abstract description 18
- 238000005516 engineering process Methods 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims abstract description 13
- 239000002245 particle Substances 0.000 claims description 34
- 230000003595 spectral effect Effects 0.000 claims description 15
- 238000010183 spectrum analysis Methods 0.000 claims description 12
- 238000004321 preservation Methods 0.000 claims description 11
- 230000006870 function Effects 0.000 claims description 9
- 230000000149 penetrating effect Effects 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000003384 imaging method Methods 0.000 claims description 6
- 239000004615 ingredient Substances 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 3
- 210000005252 bulbus oculi Anatomy 0.000 claims description 3
- 239000002131 composite material Substances 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 238000000265 homogenisation Methods 0.000 claims description 3
- 230000008676 import Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000000691 measurement method Methods 0.000 claims description 3
- 239000013028 medium composition Substances 0.000 claims description 3
- 238000013508 migration Methods 0.000 claims description 3
- 230000005012 migration Effects 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 claims description 3
- 239000004744 fabric Substances 0.000 claims description 2
- 230000009471 action Effects 0.000 claims 1
- 238000005096 rolling process Methods 0.000 claims 1
- 230000002547 anomalous effect Effects 0.000 abstract description 3
- 230000006854 communication Effects 0.000 abstract description 3
- 230000000452 restraining effect Effects 0.000 abstract description 3
- 239000000284 extract Substances 0.000 abstract description 2
- 230000008859 change Effects 0.000 description 5
- 230000010287 polarization Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009440 infrastructure construction Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/295—Means for transforming co-ordinates or for evaluating data, e.g. using computers
Abstract
A kind of geological radar echo-signal physical property destructing is known with detection target number reconstruct intelligence takes method, this method has the characteristics that frequency response and detection objective body own frequency response characteristic unique resonant using geological radar reflection wave signal in underground medium communication process, by the spectrum energy contribution rate intensity value for calculating measured target dielectric response frequency, carry out inverting and analysis, efficiently extract identification deep, faint or minor reflective signal, the interference of noise signal is removed, solution translates the spatial distribution of the spectrum energy contribution rate intensity of buried target body or anomalous body.The present invention inherently improves geologic radar detection depth and the accuracy explained, it reduces because of data interpretation personnel difference or the uncertainty of the restraining factors to geological radar result interpretation such as lack experience, solves the multi-solution that geological radar Underground interprets extremely to a certain extent, it can be achieved to identify objective body reconstruct and anomalous body feature intelligentization, so as to provide a kind of Digital Signal Processing and the new method of analysis interpretation for geologic radar detection technology.
Description
Technical field：
The present invention relates to engineering geophysics and geological mapping fields more particularly to a kind of geological radar echo-signal physical property to deconstruct
Method is taken with the intelligent knowledge of detection target number reconstruct.
Background technology：
With the fast development of urban modernization construction, large-scale skyscraper, the underground foundation of ultra-deep and crisscross multiple
The infrastructure constructions such as miscellaneous rail traffic are increasing, and geological prospecting task and the difficult geological faced also increase and dash forward increasingly
Go out, original tradition geology exploratory techniques then show slightly " helpless " in many cases with method.And with quick, lossless, high
The geological radar of the features such as effect has become the effective technology means that current engineering geological investigation field is widely used, but electromagnetism
Communication process of the wave in underground is sufficiently complex and attenuation is serious, and deep formation or small geologic body reflection signal are often very weak, respectively
The interference signal of kind noise and clutter is but very strong, and the useful reflective information for causing deep faint is further under strong jamming background
Reduction even bury or in data processing for reduce noise improve signal-to-noise ratio and excess processes cause echo-signal to damage
It loses, and then the useful information of underground deep small objects body can not be extracted and detected with to meet the big depth of geological radar and high-resolution
The needs of explanation.
It is that geological radar technology encounters " bottleneck " that effective weak reflection echo signal how is extracted in strong jamming background
With prominent problem, correctly identify that various clutters with noise are geology thunders to extract useful weak reflective information from radar return record
Up to the important link and key point of data process and interpretation.
Invention content：
The technical problem to be solved by the present invention is to the reconciliation of geologic radar detection depth can inherently be improved by providing one kind
The accuracy released, reduce because data interpretation personnel difference or lack experience etc. restraining factors to geological radar result interpretation not really
It is qualitative, and solve the multi-solution that interprets extremely of geological radar Underground to a certain extent, it realizes to objective body reconstruct and different
The geological radar echo-signal physical property of normal body characteristics Weigh sensor is deconstructed takes method with the intelligent knowledge of detection target number reconstruct.
The particular technique solution of the present invention is to provide a kind of geological radar echo-signal physical property destructing and detection target
Digital reconstruction intelligence knowledge takes method, and this method includes the following steps：
Step 1. Ground Penetrating Radar Signal acquires and data prediction；
Step 2. is based on the spectrum analysis of Fourier transformation to being included not in the geological radar echo-signal of measured target body
Destructing knowledge is carried out with medium physical property to take；
Step 3. is based on frequency of the spectrum energy analytic approach digitlization different target body in measured target medium institute's preservation space
Spectrum energy contribution rate intensity；
Step 4. reconstructs the space geometry feature of detected target body by two dimensional inversion and dimensional Modeling Technology and is assigned
Deposit the spatial distribution of underground medium；
The intelligent knowledge of step 5. takes.
Preferably, geological radar echo-signal physical property destructing according to the present invention and detection target number reconstruct intelligence
Change to know and take method, wherein, the step 1 specifically includes the following contents：
(1), using conventional geological radar and low-frequency antenna, rational acquisition parameter is set, carries out serial section at the scene
It measures and random point measures two ways and obtains reflection echo signal and initial data respectively, region measurement section and survey are surveyed meeting
As much as possible in design bore position or existing prospecting bore position or the clear and definite position of geologic anomaly in the case of point distribution
It acquires corresponding data more, while the three-dimensional coordinate information of all geological radar measuring points in area is measured using GPS；
(2), the primary reflection wave signal data for acquiring and preserving is converted using conventional Georadar Data processing software
Into the data file of " .DZT ", " * .DT " or " * .DAT " form, and to its data file editor, arrangement, storage；
(3), the methods of being interfered by filtering, zero migration, background denoising, rejecting is to the original Georadar Data of acquisition
Conventional pretreatment is carried out, tentatively to reduce noise jamming, signal-to-noise ratio is improved and obtains effective geological radar reflection echo signal.
Preferably, geological radar echo-signal physical property destructing according to the present invention and detection target number reconstruct intelligence
Change to know and take method, wherein, the step 2 specifically includes the following contents：
(1), in view of the additivity of Fourier transformation and its reciprocity with inverse Fourier transform, pass through inverse Fourier transform
The composite signal of an aperiodic interior different sine waves, conversion are decomposed during can the original radar reflection echo of collection in worksite be recorded
For the spectrum information of the sine wave signal of multiple periodically continued variations in a time-domain, thus obtained intensional data is pre-
Handling result carries out physical property destructing and identification；
(2), the sampled signal after being decomposed to above-mentioned conversion carries out frequency spectrum discerning, and according to different target medium itself, institute is peculiar
Spectral response feature and wave phase Frequency Response determine detection objective body signal, and determine its frequency spectrum characteristic of field；
(3), after known Ground Penetrating Radar carries out gridding subdivision, and to radar reflection echo signal in its cell
Carry out frequency spectrum analysis and Analyzing the amplitude spectrum；
(4), underground different target body is gone out to electromagnetism according to geology radar instruments device configuration and center of antenna frequency demarcating
The frequency response parameter of wave polarization；
(5), the amplitude spectrum and frequency bandwidth for determining detected target body standard reflection echo are analyzed, is finally separating out target
The corresponding basic frequency response characteristic peak value of body.
Preferably, geological radar echo-signal physical property destructing according to the present invention and detection target number reconstruct intelligence
Change to know and take method, wherein, the step 3 specifically includes the following contents：
(1), it is determined according to the unique distinctive dominant frequency response characteristic of the different target body isolated in step 2 itself and analysis
Frequency response range, Fourier integral calculating is carried out in the window particle of corresponding time, using mathematics definite integral to separation
The frequency spectral density function of different target body reflection echo afterwards integrates in the range of its frequency response bound and obtains different mesh
The respective response spectrum energy intensity value of standard type is calculated by formula (1), is situated between to frequency spectral density function in entire detected target
Matter window particle reflection echo range of frequency bandwidths inner product gets the space medium of detected target body window particle institute preservation
The spectrum energy intensity total value of full range domain response possessed by ingredient is calculated by formula (2),
In formula (1) and formula (2)：WnBe different detected target bodies respectively responds frequency in spectrum analysis window particle
Spectral energy intensity value, WTFor full range domain response frequency possessed by the space medium ingredient of detected target body window particle institute preservation
Spectral energy intensity total value, n are the composition medium composition of different detected target body window particles, and ω is angular frequency, and (a, b) is visits
Survey medium frequency response bound；
(2), integral and calculating program is write derived from formula (3) using C++, MATLAB computer programming language, by when window
Roll mode batch processing simultaneously calculates the window particle different target medium on radar cross-section on each period in frequency
The ratio between the integrated value in spectral domain and full k-space integrals value of frequency response in window particle reflection echo are responded, obtains geological radar
Objective body spectrum energy contribution rate intensity value in section per pass ray reflection echo window particle in different time periods,
In formula (3)：M is different detected target body spectrum energy contribution rate intensity values.
(3), judged whether to repeat " objective body and the destructing knowledge of medium physical property take " and " frequency spectrum energy according to the quality for the treatment of effect
The processing step and flow of amount contribution rate intensity value calculating "；
(4), the spectrum energy contribution rate intensity value for the different detected target bodies being calculated preserved into " * .txt " or
The data file of " * .dat " form with reference to the space coordinate of each particle, can be used for two dimensional inversion into figure, three-dimensional as needed
Modeling and imaging, anomaly analysis and interpretation.
Preferably, geological radar echo-signal physical property destructing according to the present invention and detection target number reconstruct intelligence
Change to know and take method, wherein, the step 4 specifically includes the following contents：
(1), geologic radar detection section or the coordinate of random measuring point and elevation terrain information are added, and is converted into required
Data format；
(2), after importing the spectrum energy contribution rate intensity value of corresponding coordinate and every one of different depth position, using in lattice
Grams interpolation method according to the continuous measurement method inverting of section and draws out the frequency spectrum energy of detection section underground different target body
Measure the two dimensional equivalent sectional view of contribution rate intensity distribution；
(3), the spatial dimension residing for the detected target body of underground is split into a series of netted cube cell, then
The physical parameter in web elements lattice to forming different target medium carries out linear homogenization processing, passes through dimensional Modeling Technology
Surface configuration and the spatial model of underground measured target body are rebuild, imports the coordinate of random eyeball, elevation information and every
The spectrum energy contribution rate intensity value of one measuring point different depth position, using Voxler or self-editing three-dimensional softwares using grams in lattice
It learns interpolation method and is finally inversed by the space geometry form of measured target body and the space point of underground medium spectrum energy contribution rate intensity
Cloth.
Preferably, geological radar echo-signal physical property destructing according to the present invention and detection target number reconstruct intelligence
Change to know and take method, wherein, the step 5 specifically includes the following contents：
(1), 2-d spectrum energy contribution rate intensity equivalence is cutd open according to existing geologic information, Use of Geophysical Data and borehole data
Face figure and three-dimensional spectrum energy contribution rate Strength Space distribution map are modified；
(2), according to the quality of modified result effect judge whether that two dimensional inversion and three-dimensional imaging result are repeated to correct, make
It matches and coincide with given data result to the greatest extent；
It (3), can output spectrum energy contribution rate intensity equivalence sectional view, different depth according to demand after final correct
The terrace cut slice figure of dropping cut slice figure, all directions takes the knowledge of its result intelligentization, and comprehensive in two dimensional surface and three dimensions realization
Close analysis and become more meticulous interpretation, can " * .jpg " or " * .GIFf " picture format output file, can also pass through three-dimensional result
" video is arrested " function obtains the video output File of " * .avi " form of different visual angles；
(4), by establishing the spectrum energy strength characteristic database of different target body and medium, with reference to abnormal edge search
Technology is taken with the knowledge of target body characteristics, appliance computer program intelligence knowledge takes discrimination off-note, according to borehole data or has
The comprehensive analysis such as other Use of Geophysical Data interpret, and then infer the detailed geological information of deep layer and faint useful geologic anomaly.
The beneficial effects of the invention are as follows：The present invention is using geological radar reflection echo signal in underground medium communication process
There is frequency response and detection objective body own frequency response characteristic unique resonant, responded by calculating measured target body
The spectrum energy contribution rate intensity value of frequency carries out inverting and analysis, efficiently extracts identification deep, faint or minor reflective letter
Number, the interference of noise signal is removed, solution translates the spatial distribution of the spectrum energy contribution rate intensity of buried target body or anomalous body.
Geologic radar detection depth and the accuracy explained inherently are improved, is reduced because of data interpretation personnel difference or experience not
It is different to solve geological radar Underground to the uncertainty of geological radar result interpretation to a certain extent for the restraining factors such as foot
The multi-solution often interpreted is, it can be achieved that objective body reconstruct and the intelligent automatic identification of off-note, so as to be geologic radar detection
Technology provides a kind of Digital Signal Processing and the new method of analysis interpretation.
Description of the drawings：
Fig. 1 takes method for geological radar echo-signal physical property of the present invention destructing with detecting the intelligent knowledge of target number reconstruct
Flow diagram.
Specific embodiment：
Geological radar echo-signal physical property of the present invention is deconstructed with detecting number of targets in the following with reference to the drawings and specific embodiments
The intelligent knowledge of word reconstruct takes method to be described further：
It is all respective with to one's name and uniquely in radar reflection echo signal based on underground difference electromagnetism physical property medium
Frequency resonance response this physical property structure attribute.Pass through the physical property to different target medium in geology radar reflection echo signal
Resonance frequency analysis, the frequency-response characteristic for taking out detected target body is known in destructing, and is being returned by detected target body medium
Contribution rate in wave signal spectrum quantifies its spectrum energy intensity in dielectric space particle, reconstructs the space of measured target body
The spatial distribution of geometric shape and underground medium spectrum energy contribution rate intensity is realized and detects interpretation to underground objective body and exception
Purpose.By the Ground Penetrating Radar Signal processing interpretation of the physical property intelligent recognition of this pair of of physical prospecting signal and quantum chemical method and objective body
Reconstruct intelligence, which is known, takes method, referred to as spectrum energy intensive analysis method.
The principle of the invention is sketched：
The echo-signal of most geophysical prospecting method records is rung by the different medium frequency for forming detected target body or exception
Information iteration is answered to form, by the frequency domain spectrum analysis of Fourier transform, the physical property structure of detected target body can be deconstructed, and set up
The frequency response physical property characteristic database of different target medium knows the objective body or different taken in detected medium by comparison automatically
Often.
For spectrum analysis window set on geological radar reflection echo recording ray, the window of measured target body is obtained
Particle spatial coordinate location；Collection of illustrative plates based on spectrum analysis calculates the spectrum energy contribution rate intensity of detected target body window particle
Value, with reference to the spatial distribution of detected target body window particle, inverting reconstructs geometry and the space of detected target body
Preservation situation, to reach the detection purpose of intelligent recognition.
According to above-mentioned principle, as shown in Figure 1, geological radar echo-signal physical property of the present invention destructing and detection target number weight
Structure intelligence knowledge takes method to include the following steps：
Step 1. Ground Penetrating Radar Signal acquires and data prediction.
Wherein, step 1 specifically includes the following contents：
(1), using conventional geological radar and low-frequency antenna, set rational acquisition parameter (pay special attention to time window and
The setting of stacking fold), Continuous Profiling is carried out at the scene and random point measures two ways and obtains reflection echo letter respectively
Number and initial data, surveyed in the case that region measurement section and measuring point are distributed as much as possible in design bore position or meeting
There are prospecting bore position or the clear and definite position of geologic anomaly to acquire corresponding data, while measure to survey in area using GPS and own more
The spatial coordinated information of geological radar measuring point (including X, Y coordinate and elevation H)；
(2), the primary reflection wave signal data for acquiring and preserving is converted using conventional Georadar Data processing software
Into the data file of forms such as " .DZT ", " * .DT " or " * .DAT ", and to its data file editor, arrangement, storage；
(3), the methods of being interfered by filtering, zero migration, background denoising, rejecting is to the original Georadar Data of acquisition
Conventional pretreatment is carried out, tentatively to reduce noise jamming, signal-to-noise ratio is improved and obtains effective geological radar reflection echo signal.
Step 2. is based on the spectrum analysis of Fourier transformation to being included not in the geological radar echo-signal of measured target body
Destructing knowledge is carried out with medium physical property to take.
Wherein, step 2 specifically includes the following contents：
(1), in view of the additivity of Fourier transformation and its reciprocity with inverse Fourier transform, pass through inverse Fourier transform
The composite signal of an aperiodic interior different sine waves, conversion are decomposed during can the original radar reflection echo of collection in worksite be recorded
For the spectrum information of the sine wave signal of multiple periodically continued variations in a time-domain, thus obtained intensional data is pre-
Handling result carries out physical property destructing and identification；
(2), the sampled signal after being decomposed to above-mentioned conversion carries out frequency spectrum discerning, and according to different target medium itself, institute is peculiar
Spectral response feature and wave phase Frequency Response determine detection objective body signal, and determine its frequency spectrum characteristic of field；
(3), after known Ground Penetrating Radar carries out gridding subdivision, and to radar reflection echo signal in its cell
Carry out frequency spectrum analysis and Analyzing the amplitude spectrum；
(4), underground different target body is gone out to electromagnetism according to geology radar instruments device configuration and center of antenna frequency demarcating
The frequency response parameter of wave polarization；
(5), the amplitude spectrum and frequency bandwidth for determining destination media standard reflection echo are analyzed, is finally separating out destination media
Corresponding basic frequency response characteristic peak value.
Step 3. is based on frequency of the spectrum energy analytic approach digitlization different target body in measured target medium institute's preservation space
Spectral energy intensity and ratio.
Wherein, step 3 specifically includes the following contents：
(1), according to the unique distinctive dominant frequency response characteristic (i.e. peak value) of the different target body isolated in step 2 itself and
The determining frequency response range (i.e. frequency bandwidth) of analysis carries out Fourier integral calculating in corresponding time window particle, adopts
With mathematics definite integral to the frequency spectral density function of the different target body reflection echo after separation in its frequency response bound model
The respective response spectrum energy contribution rate intensity value that inner product gets different target body is enclosed, is calculated by formula (1), it is close to frequency spectrum
Degree function gets detected target window in entire detected target medium window particle reflection echo range of frequency bandwidths inner product
The spectrum energy intensity total value of full range domain response, is calculated by formula (2) possessed by the space medium ingredient of mouth particle institute preservation,
In formula (1) and formula (2)：WnBe different detected target bodies respectively responds frequency in spectrum analysis window particle
Spectral energy intensity value, WTFor full range domain response frequency possessed by the space medium ingredient of detected target body window particle institute preservation
Spectral energy intensity total value, n are the composition medium composition of different detection objective body window particles, and ω is angular frequency, and (a, b) is detection
Medium frequency response bound；
(2), integral and calculating program is write derived from formula (3) using C++, MATLAB computer programming language, by when window
Roll mode batch processing simultaneously calculates the window particle different target medium on radar cross-section on each period in frequency
The ratio between the integrated value in spectral domain and full k-space integrals value of frequency response in window particle reflection echo are responded, obtains geological radar
Objective body spectrum energy contribution rate intensity value in section per pass ray reflection echo window particle in different time periods,
In formula (3)：M is different detected target body spectrum energy contribution rate intensity values.
(3), judged whether to repeat " objective body and the destructing knowledge of medium physical property take " and " frequency spectrum energy according to the quality for the treatment of effect
The processing step and flow of amount contribution rate intensity value calculating "；
(4), the spectrum energy contribution rate intensity value for the different detected target bodies being calculated preserved into " * .txt " or
The data file of " * .dat " form with reference to the space coordinate of each particle, can be used for two dimensional inversion into figure, three-dimensional as needed
Modeling and imaging, anomaly analysis and interpretation.
Step 4. reconstructs the space geometry feature of detected target body by two dimensional inversion and dimensional Modeling Technology and is assigned
Deposit the spatial distribution of underground medium.
Wherein, step 4 specifically includes the following contents：
(1), geologic radar detection section or the coordinate of random measuring point and elevation terrain information are added, and is converted into required lattice
Formula；
(2), the spectrum energy contribution rate intensity value of corresponding coordinate and every one of different depth position (time window) is imported
Afterwards, using grams interpolation method in lattice, according to the continuous measurement method inverting of section and detection section underground difference mesh is drawn out
The two dimensional equivalent sectional view of the spectrum energy contribution rate intensity distribution of standard type；
(3), the spatial dimension residing for the detected target body of underground is split into a series of netted cube cell (its
Size and resolution ratio are proportionate), then physical parameter in web elements lattice to forming different target medium is into line
Property homogenization processing, surface configuration and the spatial model of underground measured target body are rebuild by dimensional Modeling Technology, import
The spectrum energy contribution rate intensity of the coordinate of eyeball, elevation information and each measuring point different depth position (time window) at random
Value goes out the space geometry of measured target body using Voxler or self-editing three-dimensional softwares using mathematical interpolations inversion methods such as Ge Like
The spatial distribution of form and underground medium spectrum energy contribution rate intensity.
The intelligent knowledge of step 5. takes.
Wherein, step 5 specifically includes the following contents：
(1), 2-d spectrum energy contribution rate intensity equivalence is cutd open according to existing geologic information, Use of Geophysical Data and borehole data
Face figure and three-dimensional spectrum energy contribution rate Strength Space distribution map are modified；
(2), according to the quality of modified result effect judge whether that two dimensional inversion and three-dimensional imaging result are repeated to correct, make
It is matched with given data result to the greatest extent with coincideing；
It (3), can output spectrum energy contribution rate intensity equivalence sectional view, different depth according to demand after final correct
The terrace cut slice figure etc. of dropping cut slice figure, all directions takes the knowledge of its result intelligentization, and is realized in two dimensional surface and three dimensions
Comprehensive analysis and become more meticulous interpretation, can the plurality of picture form output file such as " * .jpg " or " * .GIFf ", can also pass through three-dimensional
As a result " video is arrested " function obtains the video output File of the multiple formats such as " the * .avi " of different visual angles；
(4), by establishing the spectrum energy strength characteristic database of different target body and medium, with reference to abnormal edge search
Technology is taken with the knowledge of target body characteristics, appliance computer program intelligence knowledge takes off-note, according to borehole data or has other
The comprehensive analysis such as Use of Geophysical Data interpret, and then infer the detailed geological information of deep layer and faint useful geologic anomaly.
Embodiment described above is only that the preferred embodiment of the present invention is described, not to the scope of the present invention
It is defined, under the premise of design spirit of the present invention is not departed from, those of ordinary skill in the art are to technical scheme of the present invention
The various modifications made and improvement should all be fallen into the protection domain that claims of the present invention determines.
Claims (6)
1. a kind of geological radar echo-signal physical property destructing is known with detection target number reconstruct intelligence takes method, feature exists
In：This method includes the following steps：
Step 1. Ground Penetrating Radar Signal acquires and data prediction；
Step 2. is situated between to the difference included in the geological radar echo-signal of measured target body based on the spectrum analysis of Fourier transformation
Matter physical property carries out destructing knowledge and takes；
Step 3. is based on frequency spectrum energy of the spectrum energy analytic approach digitlization different target body in measured target medium institute's preservation space
Measure contribution rate intensity；
Step 4. is with reconstructing space geometry feature and the institute's preservation of detected target body by two dimensional inversion and dimensional Modeling Technology
The spatial distribution of lower medium；
The intelligent knowledge of step 5. takes.
2. it is deconstructed according to geological radar echo-signal physical property described in claim 1 and knows the side of taking with detecting target number reconstruct intelligence
Method, it is characterised in that：The step 1 specifically includes the following contents：
(1), using conventional geological radar and low-frequency antenna, rational acquisition parameter is set, carries out Continuous Profiling at the scene
Two ways is measured with random point and obtains reflection echo signal and initial data respectively, and region measurement section and measuring point point are surveyed meeting
In the case of cloth as much as possible design bore position or it is existing prospecting bore position or the clear and definite position of geologic anomaly more adopt
Collect corresponding data, while the three-dimensional coordinate information of all geological radar measuring points in area is measured using GPS；
(2), the primary reflection wave signal data for acquiring and preserving is converted into using conventional Georadar Data processing software
The data file of " .DZT ", " * .DT " or " * .DAT " form, and to its data file editor, arrangement, storage；
(3), the original Georadar Data of acquisition is carried out by filtering, zero migration, background denoising, rejecting interference method normal
Rule pretreatment tentatively to reduce noise jamming, improves signal-to-noise ratio and obtains effective geological radar reflection echo signal.
3. it is deconstructed according to geological radar echo-signal physical property described in claim 1 and knows the side of taking with detecting target number reconstruct intelligence
Method, it is characterised in that：The step 2 specifically includes the following contents：
(1), in view of the additivity of Fourier transformation and its reciprocity with inverse Fourier transform, can be incited somebody to action by inverse Fourier transform
The composite signal of an aperiodic interior different sine waves, conversion are decomposed into one in the original radar reflection echo record of collection in worksite
The spectrum information of the sine wave signal of multiple periodically continued variations in a time-domain, thus obtained intensional data pretreatment
As a result physical property destructing and identification are carried out；
(2), the sampled signal after being decomposed to above-mentioned conversion carries out frequency spectrum discerning, the frequency according to specific to different target medium itself
Spectrum response characteristic and wave phase Frequency Response determine detection objective body signal, and determine its frequency spectrum characteristic of field；
(3), after known Ground Penetrating Radar carries out gridding subdivision, and radar reflection echo signal in its cell is carried out
Frequency spectrum analysis and Analyzing the amplitude spectrum；
(4), to go out underground different target body according to geology radar instruments device configuration and center of antenna frequency demarcating inclined to electromagnetic wave
The frequency response parameter shaken；
(5), the amplitude spectrum and frequency bandwidth for determining detected target body standard reflection echo are analyzed, is finally separating out objective body pair
The basic frequency response characteristic peak value answered.
4. it is deconstructed according to geological radar echo-signal physical property described in claim 1 and knows the side of taking with detecting target number reconstruct intelligence
Method, it is characterised in that：The step 3 specifically includes the following contents：
(1), according to the determining frequency of the unique distinctive dominant frequency response characteristic of the different target body isolated in step 2 itself and analysis
Rate response range carries out Fourier integral calculating in the window particle of corresponding time, using mathematics definite integral to separation after
The frequency spectral density function of different target body reflection echo integrates in the range of its frequency response bound and obtains different target body
Respective response spectrum energy intensity value, by formula (1) calculate, to frequency spectral density function in entire detected target medium window
Mouth particle reflection echo range of frequency bandwidths inner product gets the space medium ingredient of detected target body window particle institute preservation
The spectrum energy intensity total value of possessed full range domain response is calculated by formula (2),
In formula (1) and formula (2)：WnThe respective response spectrum energy for being different detected target bodies in spectrum analysis window particle
Measure intensity value, WTFor full frequency-domain response spectrum energy possessed by the space medium ingredient of detected target body window particle institute preservation
Intensity total value is measured, n is the composition medium composition of different detection objective body window particles, and ω is angular frequency, and (a, b) is detection medium
Frequency response bound；
(2), integral and calculating program is write derived from formula (3) using C++, MATLAB computer programming language, passes through time-window rolling
Mode batch processing simultaneously calculates the window particle different target medium on radar cross-section on each period in frequency response
The ratio between full k-space integrals value of frequency response in integrated value and window particle reflection echo in spectral domain, obtains Ground Penetrating Radar
Objective body spectrum energy contribution rate intensity value in per pass ray reflection echo window particle in different time periods,
In formula (3)：M is different detected target body spectrum energy contribution rate intensity values.
(3), judged whether to repeat " objective body and the destructing knowledge of medium physical property take " and " spectrum energy tribute according to the quality for the treatment of effect
Offer rate intensity value calculating " processing step and flow；
(4), the spectrum energy contribution rate intensity value for the different detected target bodies being calculated is preserved into " * .txt " or " *
.dat " the data file of form with reference to the space coordinate of each particle, can be used for two dimensional inversion to be built into figure, three-dimensional as needed
Mould and imaging, anomaly analysis and interpretation.
5. it is deconstructed according to geological radar echo-signal physical property described in claim 1 and knows the side of taking with detecting target number reconstruct intelligence
Method, it is characterised in that：The step 4 specifically includes the following contents：
(1), geologic radar detection section or the coordinate of random measuring point and elevation terrain information are added, and is converted into required data
Form；
(2), after importing the spectrum energy contribution rate intensity value of corresponding coordinate and every one of different depth position, using grams in lattice
Interpolation method is learned, according to the continuous measurement method inverting of section and draws out the spectrum energy tribute of detection section underground different target body
Offer the two dimensional equivalent sectional view of rate intensity distribution；
(3), the spatial dimension residing for the detected target body of underground is split into a series of netted cube cell, then to group
Linear homogenization processing is carried out into the physical parameter in the web elements lattice of different target medium, by dimensional Modeling Technology again
Surface configuration and the spatial model of underground measured target body are built, imports the coordinate, elevation information and each survey of random eyeball
The spectrum energy contribution rate intensity value of point different depth position is inserted using Voxler or self-editing three-dimensional softwares using grams in lattice
Value inversion method goes out the space geometry form of measured target body and the spatial distribution of underground medium spectrum energy contribution rate intensity.
6. it is deconstructed according to geological radar echo-signal physical property described in claim 1 and knows the side of taking with detecting target number reconstruct intelligence
Method, it is characterised in that：The step 5 specifically includes the following contents：
(1), according to existing geologic information, Use of Geophysical Data and borehole data to 2-d spectrum energy contribution rate intensity equivalence sectional view
It is modified with three-dimensional spectrum energy contribution rate Strength Space distribution map；
(2), according to the quality of modified result effect judge whether that two dimensional inversion and three-dimensional imaging result are repeated to correct, make its with
Given data result is matched and is coincide to the greatest extent；
It (3), can output spectrum energy contribution rate intensity equivalence sectional view, different depth be horizontal according to demand after final correct
The terrace cut slice figure of slice map, all directions takes the knowledge of its result intelligentization, and realizes comprehensive point in two dimensional surface and three dimensions
Analyse and become more meticulous interpretation, can " * .jpg " or " * .GIFf " picture format output file, can also pass through the " video of three-dimensional result
Arrest " function obtain different visual angles " * .avi " form video output File；
(4), by establishing the spectrum energy strength characteristic database of different target body and medium, with reference to abnormal edge search and mesh
The knowledge of standard type feature takes technology, and appliance computer program intelligence knowledge takes discrimination off-note, according to borehole data or has other
The comprehensive analysis such as Use of Geophysical Data interpret, and then infer the detailed geological information of deep layer and faint useful geologic anomaly.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810064500.XA CN108267722A (en) | 2018-01-23 | 2018-01-23 | Geological radar echo-signal physical property, which deconstructs to know with detection target number reconstruct intelligence, takes method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810064500.XA CN108267722A (en) | 2018-01-23 | 2018-01-23 | Geological radar echo-signal physical property, which deconstructs to know with detection target number reconstruct intelligence, takes method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108267722A true CN108267722A (en) | 2018-07-10 |
Family
ID=62776319
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810064500.XA Pending CN108267722A (en) | 2018-01-23 | 2018-01-23 | Geological radar echo-signal physical property, which deconstructs to know with detection target number reconstruct intelligence, takes method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108267722A (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109085584A (en) * | 2018-09-03 | 2018-12-25 | 电子科技大学 | Multiple-input and multiple-output borehole radar high efficiency imaging method based on high-freedom degree |
| CN109255100A (en) * | 2018-09-10 | 2019-01-22 | 河海大学 | A kind of Urban Rain inversion algorithm based on microwave attenuation characteristic response fingerprint recognition |
| CN110319916B (en) * | 2019-06-06 | 2020-06-30 | 浙江大学 | Limited pool low-frequency expanding method based on water area reflected wave interference measurement |
| CN112540412A (en) * | 2020-11-18 | 2021-03-23 | 清华大学 | Target detection method, device, equipment and storage medium |
| CN112731048A (en) * | 2020-12-08 | 2021-04-30 | 深圳供电局有限公司 | Voltage sag detection method and device, computer equipment and readable storage medium |
| CN112965119A (en) * | 2021-02-03 | 2021-06-15 | 安徽圭目机器人有限公司 | Prediction convolution processing method based on three-dimensional ground penetrating radar |
| CN113376612A (en) * | 2021-08-12 | 2021-09-10 | 成都众享天地网络科技有限公司 | Radar clutter generation method based on terrain matrixing and detection |
| CN113759422A (en) * | 2021-09-10 | 2021-12-07 | 刘俊伟 | Underground abnormal body detection method |
| CN114488151A (en) * | 2022-04-08 | 2022-05-13 | 中国科学院空天信息创新研究院 | Active and passive combined detection method, device, equipment and medium for observation ship |
| CN117092644A (en) * | 2023-10-17 | 2023-11-21 | 中国电建集团江西省电力设计院有限公司 | Underground pipeline detection method |
Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH01113685A (en) * | 1987-10-28 | 1989-05-02 | Nippon Telegr & Teleph Corp <Ntt> | Object detecting method and apparatus |
| US4965582A (en) * | 1984-11-24 | 1990-10-23 | Hellsten Hans O | Method for radar mapping an area and a radar equipment to carry out the method |
| US5796363A (en) * | 1996-03-01 | 1998-08-18 | The Regents Of The University Of California | Automatic position calculating imaging radar with low-cost synthetic aperture sensor for imaging layered media |
| JPH11258341A (en) * | 1998-03-13 | 1999-09-24 | Nec Corp | Underground exploring radar device and method |
| CA2456459A1 (en) * | 2003-06-16 | 2004-12-16 | Andrew M. Mcdermott | Method for enhancing depth and spatial resolution of one and two dimensional residual surfaces derived from scalar potential data |
| RU2010129261A (en) * | 2010-07-14 | 2012-01-20 | Евгений Васильевич Биряльцев (RU) | METHOD OF GEOLOGICAL EXPLORATION |
| CN102608223A (en) * | 2012-03-06 | 2012-07-25 | 张维平 | Digital signal analysis method based on spectral energy ratio |
| CN102866429A (en) * | 2012-04-28 | 2013-01-09 | 中国神华能源股份有限公司 | Method for determining groundwater occurrence |
| CN103424769A (en) * | 2013-04-26 | 2013-12-04 | 上海申丰地质新技术应用研究所有限公司 | Combined multi-wave seismic exploration method in gob |
| CN103675922A (en) * | 2013-12-13 | 2014-03-26 | 南京工业大学 | Operation period underground pipeline pipe diameter measuring method based on ground penetrating radar |
| CN103941254A (en) * | 2014-03-03 | 2014-07-23 | 中国神华能源股份有限公司 | Soil physical property classification recognition method and device based on geological radar |
| CN104020495A (en) * | 2014-06-24 | 2014-09-03 | 中国矿业大学(北京) | Automatic underground pipeline parameter recognizing method on basis of ground penetrating radar |
| US20150054670A1 (en) * | 2010-10-27 | 2015-02-26 | Jianqi Wang | Multichannel UWB-based radar life detector and positioning method thereof |
| CN105137498A (en) * | 2015-09-17 | 2015-12-09 | 鲁东大学 | Underground target detection and recognition system and method based on feature fusion |
| RU2014121964A (en) * | 2014-05-29 | 2015-12-20 | Федеральное государственное казенное военное образовательное учреждение высшего профессионального образования "Военная академия войсковой противовоздушной обороны Вооруженных Сил Российской Федерации имени Маршала Советского Союза А.М. Василевского" Министерства Обороны Российской Федерации | METHOD FOR EXPERIMENTAL CHECK OF INFORMATION AND IDENTIFICATION POSSIBILITIES OF DOPPLER PORTRAITS OF AIR OBJECTS |
-
2018
- 2018-01-23 CN CN201810064500.XA patent/CN108267722A/en active Pending
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4965582A (en) * | 1984-11-24 | 1990-10-23 | Hellsten Hans O | Method for radar mapping an area and a radar equipment to carry out the method |
| JPH01113685A (en) * | 1987-10-28 | 1989-05-02 | Nippon Telegr & Teleph Corp <Ntt> | Object detecting method and apparatus |
| US5796363A (en) * | 1996-03-01 | 1998-08-18 | The Regents Of The University Of California | Automatic position calculating imaging radar with low-cost synthetic aperture sensor for imaging layered media |
| JPH11258341A (en) * | 1998-03-13 | 1999-09-24 | Nec Corp | Underground exploring radar device and method |
| CA2456459A1 (en) * | 2003-06-16 | 2004-12-16 | Andrew M. Mcdermott | Method for enhancing depth and spatial resolution of one and two dimensional residual surfaces derived from scalar potential data |
| RU2010129261A (en) * | 2010-07-14 | 2012-01-20 | Евгений Васильевич Биряльцев (RU) | METHOD OF GEOLOGICAL EXPLORATION |
| US20150054670A1 (en) * | 2010-10-27 | 2015-02-26 | Jianqi Wang | Multichannel UWB-based radar life detector and positioning method thereof |
| CN102608223A (en) * | 2012-03-06 | 2012-07-25 | 张维平 | Digital signal analysis method based on spectral energy ratio |
| CN102866429A (en) * | 2012-04-28 | 2013-01-09 | 中国神华能源股份有限公司 | Method for determining groundwater occurrence |
| CN103424769A (en) * | 2013-04-26 | 2013-12-04 | 上海申丰地质新技术应用研究所有限公司 | Combined multi-wave seismic exploration method in gob |
| CN103675922A (en) * | 2013-12-13 | 2014-03-26 | 南京工业大学 | Operation period underground pipeline pipe diameter measuring method based on ground penetrating radar |
| CN103941254A (en) * | 2014-03-03 | 2014-07-23 | 中国神华能源股份有限公司 | Soil physical property classification recognition method and device based on geological radar |
| RU2014121964A (en) * | 2014-05-29 | 2015-12-20 | Федеральное государственное казенное военное образовательное учреждение высшего профессионального образования "Военная академия войсковой противовоздушной обороны Вооруженных Сил Российской Федерации имени Маршала Советского Союза А.М. Василевского" Министерства Обороны Российской Федерации | METHOD FOR EXPERIMENTAL CHECK OF INFORMATION AND IDENTIFICATION POSSIBILITIES OF DOPPLER PORTRAITS OF AIR OBJECTS |
| CN104020495A (en) * | 2014-06-24 | 2014-09-03 | 中国矿业大学(北京) | Automatic underground pipeline parameter recognizing method on basis of ground penetrating radar |
| CN105137498A (en) * | 2015-09-17 | 2015-12-09 | 鲁东大学 | Underground target detection and recognition system and method based on feature fusion |
Non-Patent Citations (2)
| Title |
|---|
| 李耀华等: "地质雷达信号谱能分析在隧道基底深层地质异常探查中的应用", 《第十五届全国工程物探与岩土工程测试学术大会论文集》 * |
| 李耀华等: "地质雷达信号谱能分析法在混凝土缺损探伤中的应用", 《物探化探计算技术》 * |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109085584A (en) * | 2018-09-03 | 2018-12-25 | 电子科技大学 | Multiple-input and multiple-output borehole radar high efficiency imaging method based on high-freedom degree |
| CN109255100A (en) * | 2018-09-10 | 2019-01-22 | 河海大学 | A kind of Urban Rain inversion algorithm based on microwave attenuation characteristic response fingerprint recognition |
| CN110319916B (en) * | 2019-06-06 | 2020-06-30 | 浙江大学 | Limited pool low-frequency expanding method based on water area reflected wave interference measurement |
| CN112540412A (en) * | 2020-11-18 | 2021-03-23 | 清华大学 | Target detection method, device, equipment and storage medium |
| CN112731048A (en) * | 2020-12-08 | 2021-04-30 | 深圳供电局有限公司 | Voltage sag detection method and device, computer equipment and readable storage medium |
| CN112965119A (en) * | 2021-02-03 | 2021-06-15 | 安徽圭目机器人有限公司 | Prediction convolution processing method based on three-dimensional ground penetrating radar |
| CN113376612A (en) * | 2021-08-12 | 2021-09-10 | 成都众享天地网络科技有限公司 | Radar clutter generation method based on terrain matrixing and detection |
| CN113376612B (en) * | 2021-08-12 | 2021-11-23 | 成都众享天地网络科技有限公司 | Radar clutter generation method based on terrain matrixing and detection |
| CN113759422A (en) * | 2021-09-10 | 2021-12-07 | 刘俊伟 | Underground abnormal body detection method |
| CN113759422B (en) * | 2021-09-10 | 2023-11-21 | 刘俊伟 | Underground abnormal body detection method |
| CN114488151A (en) * | 2022-04-08 | 2022-05-13 | 中国科学院空天信息创新研究院 | Active and passive combined detection method, device, equipment and medium for observation ship |
| CN114488151B (en) * | 2022-04-08 | 2022-06-24 | 中国科学院空天信息创新研究院 | Active and passive combined detection method, device, equipment and medium for observation ship |
| CN117092644A (en) * | 2023-10-17 | 2023-11-21 | 中国电建集团江西省电力设计院有限公司 | Underground pipeline detection method |
| CN117092644B (en) * | 2023-10-17 | 2024-01-26 | 中国电建集团江西省电力设计院有限公司 | Underground pipeline detection method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108267722A (en) | Geological radar echo-signal physical property, which deconstructs to know with detection target number reconstruct intelligence, takes method | |
| CN113759425B (en) | Method and system for evaluating filling characteristics of deep paleo-karst reservoir stratum by well-seismic combination | |
| Paasche et al. | Integration of diverse physical-property models: Subsurface zonation and petrophysical parameter estimation based on fuzzy c-means cluster analyses | |
| Zeng et al. | High-frequency sequence stratigraphy from seismic sedimentology: Applied to Miocene, Vermilion Block 50, Tiger Shoal area, offshore Louisiana | |
| de Matos et al. | Unsupervised seismic facies analysis using wavelet transform and self-organizing maps | |
| US5461562A (en) | Method and apparatus for detecting and quantifying hydrocarbon bearing laminated reservoirs on a workstation | |
| CN107356967B (en) | A kind of compacting seismic data shields by force the sparse optimization method of interference | |
| CN110095773A (en) | The two-parameter inversion method of the multiple dimensioned Full wave shape of Ground Penetrating Radar | |
| CN103792573A (en) | Seismic wave impedance inversion method based on frequency spectrum fusion | |
| Zhou et al. | Migration velocity analysis and prestack migration of common-transmitter GPR data | |
| Bruno et al. | Groundwater characterization in arid regions using seismic and gravity attributes: Al Jaww Plain, UAE | |
| Martorana et al. | Assessment of a geological model by surface wave analyses | |
| Schuelke et al. | Reservoir architecture and porosity distribution, Pegasus field, West Texas—an integrated sequence stratigraphic-seismic attribute study using neural networks | |
| Zhao et al. | A Comprehensive Horizon‐Picking Method on Subbottom Profiles by Combining Envelope, Phase Attributes, and Texture Analysis | |
| CN104570090B (en) | The extraction of full waveform inversion noise filter operator and the method filtered using its noise | |
| Shaw et al. | Robust description of statistically heterogeneous seafloor topography through its slope distribution | |
| Hyndman et al. | Hydrogeophysical case studies at the local scale: The saturated zone | |
| Carpentier et al. | Semiautomated suppression of above-surface diffractions in GPR data | |
| CN115857028A (en) | Time-space quantitative expression method of hydrogeological parameters based on non-contact observation method | |
| CN110261905A (en) | Complex value based on pitch angle control is concerned with microfault recognition methods | |
| Saleh et al. | Seismic-electromagnetic projection attribute: Application in integrating seismic quantitative interpretation and 3D controlled-source electromagnetic-magnetotelluric broadband data inversion for robust ranking and sweet spotting of hydrocarbon prospects in offshore northwest Borneo | |
| Zhao et al. | Automatic events extraction in pre-stack seismic data based on edge detection in slant-stacked peak amplitude profiles | |
| CN105589096A (en) | Sedimentary facies belt dividing method based on D-S evidence theory | |
| Abdulsalam et al. | Evaluation of structural framework and depth estimates using High resolution Airborne magnetic data over some parts of Middle Benue Trough, Nigeria | |
| Song et al. | Improved lithospheric seismic velocity and density model of the Korean Peninsula from ambient seismic noise data using machine learning |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180710 |
|
| WD01 | Invention patent application deemed withdrawn after publication |