RU2169382C1 - Method of multilevel vibration-seismic prospecting - Google Patents
Method of multilevel vibration-seismic prospecting Download PDFInfo
- Publication number
- RU2169382C1 RU2169382C1 RU2000113639/28A RU2000113639A RU2169382C1 RU 2169382 C1 RU2169382 C1 RU 2169382C1 RU 2000113639/28 A RU2000113639/28 A RU 2000113639/28A RU 2000113639 A RU2000113639 A RU 2000113639A RU 2169382 C1 RU2169382 C1 RU 2169382C1
- Authority
- RU
- Russia
- Prior art keywords
- level
- filters
- depth
- seismic
- reference signal
- Prior art date
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
Изобретение относится к вибросейсморазведке и предназначено для повышения качества и информативности сейсмического материала. The invention relates to vibroseismic exploration and is intended to improve the quality and information content of seismic material.
Известен способ сейсмической разведки в соответствии, с которым функцию вибратора пропускают через фильтр, непрерывно настраиваемый на мгновенную частоту вибратора после перемножения ее с регистрируемыми колебаниями фильтруют на полосовом фильтре, частотную характеристику которого, перемещают по оси частот на величину, пропорциональную отклонению скорости излучения частоты вибратора от своего среднего значения [1]. There is a known method of seismic exploration in accordance with which the vibrator function is passed through a filter that is continuously tuned to the instantaneous frequency of the vibrator after it is multiplied with registered vibrations and filtered on a band-pass filter, the frequency response of which is moved along the frequency axis by an amount proportional to the deviation of the frequency of the vibrator from the vibrator its average value [1].
Недостатком этого способа является то, что он учитывает лишь искажения вибрационных сигналов, возникающие в системе вибратор-грунт, а также искажающее влияние паразитных гармоник. Искажения же вибрационного сигнала в среде, связанные с частотно-зависимым затуханием и дисперсией скорости, не учитываются. Это снижает эффективность способа. Кроме того, он требует изготовления сложной аппаратуры на базе двух-трехдорожечного магнитофона специального типа, перемножающих устройств и линий задержки. The disadvantage of this method is that it takes into account only distortion of vibration signals that occur in the vibrator-soil system, as well as the distorting effect of spurious harmonics. Distortions of the vibrational signal in the medium associated with frequency-dependent attenuation and velocity dispersion are not taken into account. This reduces the effectiveness of the method. In addition, it requires the manufacture of sophisticated equipment based on a two-three-track tape recorder of a special type, multiplying devices and delay lines.
Известен способ формирования сейсмических данных с использованием сейсмического вибратора и его варианты. Согласно этому способу формируют первую последовательность каскадированной развертки, содержащую n идентичных сегментов, и вторую такую же последовательность с дополнительным сегментом, который расположен и фазирован так, чтобы подавлять гармонические помехи во время корреляции [2]. A known method of generating seismic data using a seismic vibrator and its variants. According to this method, the first cascaded scan sequence is formed containing n identical segments, and the second same sequence with an additional segment that is located and phased so as to suppress harmonic interference during correlation [2].
Недостатком этого способа является то, что он, подавляя искажения коррелограмм, связанные с гармоническими помехами при вычислении интеграла корреляции, не учитывает отрицательное влияние на результаты корреляции искажений свип-сигналов в среде. Это снижает эффективность выделения импульсных сигналов на фоне помех. Кроме того, способ характеризуется сложной технологией и большой трудоемкостью, требуя применения нескольких разверток при отработке одной физической точки. The disadvantage of this method is that, while suppressing correlogram distortions associated with harmonic noise in the calculation of the correlation integral, it does not take into account the negative effect on the correlation results of sweep distortions in the medium. This reduces the efficiency of the allocation of pulsed signals against interference. In addition, the method is characterized by complex technology and high complexity, requiring the use of multiple scans when working out one physical point.
Наиболее близким к изобретению по технической сущности является способ вибросейсмической разведки и устройство для его осуществления. В соответствии с этим способом используемый при корреляции опорный сигнал дополнительно обрабатывают путем формирования на его основе первой, второй и т.д. гармоник, которые суммируют с исходным опорным сигналом и сумму этих сигналов коррелируют с виброграммами. При этом информация гармоник определяется их весовыми коэффициентами, в результате перебора которых отбирается информация об исследуемой геологической среде [3]. Closest to the invention in technical essence is a method of vibroseismic exploration and a device for its implementation. In accordance with this method, the reference signal used in correlation is further processed by forming the first, second, etc., based on it. harmonics that sum with the original reference signal and the sum of these signals correlate with vibrograms. In this case, the information of harmonics is determined by their weight coefficients, as a result of which enumeration information on the studied geological environment is selected [3].
Недостатком этого способа является то, что он не учитывает проявление эффектов частотно-зависимого затухания и зависимость скорости распространения гармоник сигналов от частоты. В результате нарушаются предпосылки метода вибрационной сейсморазведки, основанной на представлении о том, что регистрируемый сигнал подобен по форме излучаемому. Вследствие же частотно-зависимого затухания и дисперсии скоростей регистрируемые и излученные сигналы существенно отличаются, что ведет к снижению геологической эффективности вибросейсморазведки. The disadvantage of this method is that it does not take into account the manifestation of the effects of frequency-dependent attenuation and the dependence of the propagation velocity of the harmonics of the signals on the frequency. As a result, the premises of the vibrational seismic survey method are violated, based on the notion that the recorded signal is similar in shape to the radiated one. Due to the frequency-dependent attenuation and velocity dispersion, the recorded and emitted signals are significantly different, which leads to a decrease in the geological efficiency of vibroseismic exploration.
Технической задачей, решаемой изобретением, является повышение качества и информативности сейсмического материала за счет обеспечения оптимальных условий выделения сейсмических сигналов для различных уровней глубин. The technical problem solved by the invention is to improve the quality and information content of seismic material by providing optimal conditions for the seismic signals for various levels of depth.
Поставленная задача решается следующим образом:
В способе многоуровневой вибросейсморазведки, включающем возбуждение сейсмических колебаний виброисточником, их регистрацию, обработку и корреляцию опорного сигнала с виброграммой, выделяют подлежащие изучению уровни глубин и соответствующие им временные окна на виброграммах, перебирают для каждого из уровней значения характеристик частотно-зависимого затухания и дисперсии скоростей, вычисляют по ним характеристики фильтров, описывающих искажения сейсмических сигналов в среде, фильтруют этими фильтрами излученный вибратором сигнал, после чего вычисляют пробные коррелограммы и определяют параметры фильтров, оптимальных для каждого из уровней, по максимальной интенсивности подлежащих изучению волн, затем осуществляют окончательную корреляцию выброграммы с опорным сигналом после его предварительной обработки оптимальными для каждого уровня глубин фильтрами.The problem is solved as follows:
In a multi-level vibroseismic exploration method, including excitation of seismic vibrations by a vibration source, their registration, processing and correlation of the reference signal with the vibrogram, the depth levels and the corresponding time windows on the vibrograms to be studied are selected, the frequency-dependent attenuation and velocity dispersion characteristics are selected for each level, they calculate the characteristics of filters describing the distortion of seismic signals in the medium, filter the sig emitted by the vibrator with these filters al, after which the calculated tentative correlogram and determining filter parameters optimal for each of the levels of the maximum intensity of the waves to be researched and then carried vybrogrammy final correlation with the reference signal after pre-processing optimum for each level of depth filters.
Существенными отличиями изобретения являются следующие. Significant differences of the invention are as follows.
На виброграммах выделяют подлежащие изучению уровни глубин и соответствующие временные окна. Уровни глубин выбирают исходя из поставленных задач в соответствии с целями исследований, что позволяет сократить затраты труда и машинного времени на вычисления коррелограмм. On vibration programs, the depth levels and the corresponding time windows are selected for study. Depth levels are selected on the basis of the tasks in accordance with the objectives of the study, which reduces the cost of labor and computer time for calculating correlograms.
Перебирают для каждого из уровней значения характеристик частотно-зависимого затухания и дисперсии скоростей, вычисляют по ним характеристики фильтров, описывающих искажения сейсмических сигналов в среде, фильтруют этими фильтрами излученный вибратором сигнал, после чего вычисляют пробные коррелограммы и определяют параметры фильтров, оптимальных для каждого из уровней, по максимальной интенсивности подлежащих изучению волн. Это позволяет определить характеристики среды для всех подлежащих изучению уровней глубин, необходимые для вычисления параметров фильтров, использующихся для предваряющей корреляцию обработки свип-сигналов. For each of the levels, the values of the characteristics of the frequency-dependent attenuation and velocity dispersion are selected, the characteristics of the filters describing the distortion of seismic signals in the medium are calculated from them, the signal emitted by the vibrator is filtered with these filters, then the test correlograms are calculated and the filter parameters are determined that are optimal for each of the levels , according to the maximum intensity of the waves to be studied. This allows you to determine the characteristics of the environment for all the depth levels to be studied, which are necessary for calculating the filter parameters used to precede the correlation of the processing of sweep signals.
Осуществляют окончательную корреляцию виброграммы с опорным сигналом после его предварительной обработки оптимальными для каждого уровня глубин фильтрами, что позволяет получить максимальную амплитуду функции взаимной корреляции для подлежащих изучению волн во всех интервалах глубин. Тем самым обеспечивается высокое качество и информативность сейсмического материала для всех изучаемых уровней глубин. The final correlation of the vibrogram with the reference signal is carried out after its preliminary processing by filters that are optimal for each depth level, which allows one to obtain the maximum amplitude of the cross-correlation function for the waves to be studied in all depth intervals. This ensures high quality and informativeness of seismic material for all studied depth levels.
Из изученной научно-технической и патентной литературы автору не известно техническое решение с перечисленной совокупностью отличительных признаков, это дает основание сделать вывод о соответствии заявляемого объекта критериям изобретения. From the studied scientific, technical and patent literature, the author does not know the technical solution with the listed set of distinctive features, this gives reason to conclude that the claimed object meets the criteria of the invention.
На фиг. 1 и 2 приведены примеры искажений формы и спектра отраженного опорного (свип) сигнала: а) исходный опорный сигнал, б) амплитудный спектр исходного опорного сигнала, в) отраженный от границы опорный сигнал, г) амплитудный спектр отраженного сигнала. In FIG. Figures 1 and 2 give examples of distortions in the shape and spectrum of the reflected reference (sweep) signal: a) the original reference signal, b) the amplitude spectrum of the original reference signal, c) the reference signal reflected from the boundary, d) the amplitude spectrum of the reflected signal.
На фиг. 3 и 4 примеры изменения отраженного сигнала на коррелограммах. In FIG. 3 and 4 are examples of changes in the reflected signal in correlograms.
На фиг. 5 относительная интенсивность сигнала при значении декремента затухания в среде δ = 0,02 (а) и 0,05 (б). In FIG. 5, the relative intensity of the signal with the attenuation decrement in the medium δ = 0.02 (a) and 0.05 (b).
Обоснование реализации описанного способа базируется на следующих теоретических положениях. В вибросейсморазведке возбуждение упругих колебаний осуществляется с помощью вибратора, посылающего в течение интервала времени T в среду опорные (в иностранной литературе свип) сигналы. Наибольшее распространение в практике нашли линейные частотно-модулированные сигналы. В общем случае они записываются в виде
S(t) = A(t)cos{2π[fн+(fк-fн)t/2T]t+Φ(t)}, (1)
Спектр этого сигнала определяется прямым преобразованием Фурье
Здесь A(t) - амплитуда сигнала; fн, fк - начальная и конечная частота; ΔF - полоса частот; T - длительность сигнала; Φ(t) - функция, определяющая фазу сигнала.The rationale for the implementation of the described method is based on the following theoretical provisions. In vibroseismic exploration, the excitation of elastic vibrations is carried out using a vibrator that sends reference signals (in the foreign literature, sweep) during the time interval T to the medium. The most common in practice are linear frequency-modulated signals. In general, they are written as
S (t) = A (t) cos {2π [f n + (f to -f n ) t / 2T] t + Φ (t)}, (1)
The spectrum of this signal is determined by the direct Fourier transform
Here A (t) is the signal amplitude; f n, f k - initial and final frequency; ΔF is the frequency band; T is the signal duration; Φ (t) is a function that determines the phase of the signal.
Обязательной процедурой при обработке записей, полученных с помощью вибрационных источников (виброграмм), является вычисление функции взаимной корреляции виброграммы и опорного сигнала. В результате получают записи - коррелограммы, сходные с записями от импульсных источников. Функция взаимной корреляции принимает максимальное значение, когда сдвиг виброграммы оказывается равным времени пробега одной из волн, т.е. когда колебания одной из волн, содержащихся в полевой записи, совпадают по времени с колебаниями свип-сигнала. Вычисления функции взаимной корреляции осуществляются по формуле:
где FBK(t) - функция взаимной корреляции (ФBK); SΣ(t) - результирующая запись - виброграмма.A mandatory procedure for processing records obtained using vibration sources (vibrograms) is to calculate the cross-correlation function of the vibrogram and the reference signal. The result is records - correlograms similar to those from pulsed sources. The cross-correlation function takes its maximum value when the shift of the vibrogram is equal to the travel time of one of the waves, i.e. when the oscillations of one of the waves contained in the field recording coincide in time with the oscillations of the sweep signal. The cross-correlation function is calculated according to the formula:
where F BK (t) is the cross-correlation function (ФBK); S Σ (t) - resulting record - vibrogram.
Практическое использование вышеприведенных формул базируется на предположении о постоянстве формы опорного сигнала для всех изучаемых границ раздела среды и всех интервалов разреза. Поэтому при вычислении коррелограмм используется неизменный для всех временных интервалов опорный сигнал. Однако сигналы, распространяясь в реальных средах, испытывают значительные искажения. В общем виде импульс после прохождения пути длиной l представляется следующим интегралом Фурье:
где ψ(l,f) - спектральная характеристика среды; ψи(f) - спектр импульса при 1 = 0.The practical use of the above formulas is based on the assumption that the shape of the reference signal is constant for all studied medium interfaces and all section intervals. Therefore, when calculating correlograms, the reference signal that is invariable for all time intervals is used. However, signals propagating in real environments experience significant distortion. In general terms, a pulse after passing a path of length l is represented by the following Fourier integral:
where ψ (l, f) is the spectral characteristic of the medium; ψ and (f) is the momentum spectrum at 1 = 0.
Спектральная характеристика среды может быть записана в виде
ψ(l,f) = Φ(l,f)eiX(l,f). (5)
А. Г. Авербухом [4] показано, что в средах с частотно-зависимым затуханием и сопутствующей ему дисперсией скоростей функция ψ(l,f) определяется следующим соотношением
Здесь δ декремент поглощения.The spectral characteristic of the medium can be written as
ψ (l, f) = Φ (l, f) e iX (l, f) . (5)
A. G. Averbukh [4] showed that in media with frequency-dependent damping and the accompanying velocity dispersion, the function ψ (l, f) is determined by the following relation
Here δ is the absorption decrement.
Подставляя (2) и (6) в (4), получим
Формула (7) позволяет количественно оценить искажения формы опорного сигнала после прохождения им некоторого пути 1 в разрезе. Эти искажения увеличиваются с ростом увеличения пути пробега волны. Поэтому использование в формуле (4) исходной формы опорного сигнала вместо истинной искаженной (см. формулу 7) ведет к снижению эффективности вибросейсморазведки, особенно для глубоких горизонтов.Substituting (2) and (6) in (4), we obtain
Formula (7) allows you to quantify the distortion of the shape of the reference signal after it passes a
Для обеспечения оптимальных условий выделения сигналов для всех уровней глубин и, соответственно, повышения геологической эффективности вибросейсморазведки предлагается при вычислении коррелограмм изменять опорный сигнал в соответствии с искажениями, определяемыми частотно-зависимыми свойствами среды и дисперсией скорости. Необходимые для оптимального изменения опорного сигнала характеристики среды предлагается определять путем подбора величины δ V(f0) в формуле (3) по критерию максимальной интенсивности отражений на пробных коррелограммах, вычисляемых для различных уровней.To ensure optimal conditions for the selection of signals for all depth levels and, consequently, to increase the geological efficiency of vibroseismic exploration, it is proposed to calculate the reference signal in calculating correlograms in accordance with distortions determined by the frequency-dependent properties of the medium and velocity dispersion. It is proposed to determine the environmental characteristics necessary for the optimal change in the reference signal by selecting δ V (f 0 ) in formula (3) according to the criterion for the maximum intensity of reflections on test correlograms calculated for different levels.
Эффективность заявляемого способа доказывается результатами сейсмогеологического моделирования разреза с частотно-зависимым затуханием. Анализировалась пятислойная сейсмогеологическая модель разреза с горизонтальными границами раздела и однородными слоями. Коэффициенты отражения для всех границ выбраны равными и составляли 0,5. Глубины границ раздела, отвечающих серединам четырех выбранных уровней глубин, составляли 1 км, 2 км, 3 км и 4 км. Первый уровень отвечал уровню глубин 0,5 км - 1,5 км, второй 1,5 км - 2,5 км и т.д. Значения декремента затухания в среде варьировались в пределах от 0,02 до 0,05, что типично для осадочных пород [4]. В качестве исходного был выбран опорный ЛЧМ сигнал длительностью 1 с, начальной частотой 10 Гц и конечной частотой 50 Гц. На фиг. 1 и фиг 2 приведены примеры искажения формы и амплитудного спектра исходного опорного сигнала после его отражения от границы, расположенной на глубине 1 км при значениях δ = 0,02 и δ = 0,05 соответственно. Сопоставляя форму сигнала до и после прохождения среды можно отметить значительные его искажения, особенно в хвостовой части, проявляющиеся в значительном затухании интенсивности высокочастотных составляющих сигнала, сдвиге одноименных экстремумов во времени, особенно заметные на фиг. 2. Следствием этих искажений является снижение эффективности вибросейсморазведки, что показывается на примере изменения максимальных амплитуд отражений на коррелограммах (фиг. 3,4). Как следует из анализа фиг. 3 и фиг. 4 неучет искажений формы опорного сигнала в среде при вычислении коррелограмм приводит к значительному уменьшению интенсивности выделяемых сейсмических сигналов, особенно, для глубокозалегающих уровней. Так, с увеличением глубины с 1-го км до 4-х км при δ = 0,02 амплитуда основного максимума уменьшилась в 1,55 раза (см. фиг. 3), а при δ = 0,05 это уменьшение составило 6,5 раза (фиг. 4). В случае учета эффекта частотно-зависимого затухания и дисперсии скоростей при корреляции наблюдается значительное увеличение интенсивности отраженных сейсмических сигналов на коррелограммах (фиг. 5). Относительный выигрыш в интенсивности варьирует от 1,1 для границы с глубиной 1 км до 11 для границы с глубиной 4 км, что значительно повышает качество и информативность сейсмического материала за счет увеличения соотношения сигнал/помеха. The effectiveness of the proposed method is proved by the results of seismic-geological modeling of the section with frequency-dependent attenuation. A five-layer seismic-geological model of the section with horizontal interfaces and homogeneous layers was analyzed. Reflection coefficients for all boundaries were chosen equal and amounted to 0.5. The depths of the interfaces corresponding to the midpoints of the four selected depth levels were 1 km, 2 km, 3 km, and 4 km. The first level corresponded to a depth level of 0.5 km - 1.5 km, the second 1.5 km - 2.5 km, etc. The attenuation decrement in the medium ranged from 0.02 to 0.05, which is typical for sedimentary rocks [4]. The reference LFM signal with a duration of 1 s, an initial frequency of 10 Hz, and a final frequency of 50 Hz was chosen as the initial one. In FIG. 1 and 2, examples of distortion of the shape and amplitude spectrum of the original reference signal after its reflection from the boundary located at a depth of 1 km at values of δ = 0.02 and δ = 0.05, respectively, are given. Comparing the waveform before and after the passage of the medium, significant distortions can be noted, especially in the tail, which are manifested in a significant attenuation of the intensity of the high-frequency components of the signal, a shift of the same extremes in time, especially noticeable in FIG. 2. The consequence of these distortions is a decrease in the efficiency of vibroseismic exploration, which is shown by the example of a change in the maximum amplitudes of the reflections in the correlograms (Fig. 3.4). As follows from the analysis of FIG. 3 and FIG. 4 neglecting distortions in the shape of the reference signal in the medium when calculating correlograms leads to a significant decrease in the intensity of the extracted seismic signals, especially for deep-seated levels. So, with increasing depth from 1 km to 4 km at δ = 0.02, the amplitude of the main maximum decreased 1.55 times (see Fig. 3), and at δ = 0.05 this decrease was 6, 5 times (Fig. 4). In the case of taking into account the effect of frequency-dependent attenuation and velocity dispersion during correlation, a significant increase in the intensity of the reflected seismic signals in the correlograms is observed (Fig. 5). The relative gain in intensity varies from 1.1 for a border with a depth of 1 km to 11 for a border with a depth of 4 km, which significantly increases the quality and information content of seismic material by increasing the signal-to-noise ratio.
Литература
1. А.С. СССР N 185503, публ. 09.10.67 г. БИ N 21 - аналог.Literature
1. A.S. USSR N 185503, publ. 10/09/67, BI N 21 - analogue.
2. Патент РФ N2126983, публ. 27.02.99 г. БИ N 6 - аналог. 2. RF patent N2126983, publ. 02/27/99, BI N 6 - analogue.
3. Патент РФ N2143713, публ. 27.12.99 г. БИ N 36 - прототип. 3. RF patent N2143713, publ. 12/27/99, BI N 36 - prototype.
4. Авербух А.Г. Изучение состава и свойств горных пород при сейсморазведке. - М.: Недра, 1982. - 232 с. 4. Averbukh A.G. Study of the composition and properties of rocks during seismic exploration. - M .: Nedra, 1982. - 232 p.
Claims (1)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2000113639/28A RU2169382C1 (en) | 2000-05-29 | 2000-05-29 | Method of multilevel vibration-seismic prospecting |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2000113639/28A RU2169382C1 (en) | 2000-05-29 | 2000-05-29 | Method of multilevel vibration-seismic prospecting |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2169382C1 true RU2169382C1 (en) | 2001-06-20 |
Family
ID=20235449
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2000113639/28A RU2169382C1 (en) | 2000-05-29 | 2000-05-29 | Method of multilevel vibration-seismic prospecting |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2169382C1 (en) |
-
2000
- 2000-05-29 RU RU2000113639/28A patent/RU2169382C1/en not_active IP Right Cessation
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lang et al. | Estimating slowness dispersion from arrays of sonic logging waveforms | |
Park et al. | Multichannel analysis of surface waves | |
RU2579164C1 (en) | Handling method for determining quality of geologic environment | |
US3895343A (en) | Apparatus for producing adaptive pilot signals | |
EA007911B1 (en) | Shaped high frequency vibratory source | |
CN105785440A (en) | Dispersion curve extracting method for mine channel wave double-component seismic signal | |
US20050152220A1 (en) | Seismic data processing | |
US3350683A (en) | Optimum detector type and depth in marine seismic exploration | |
Baradello et al. | Fast method to transform chirp envelope data into pseudo-seismic data | |
RU2593782C1 (en) | Method for vibration seismic survey | |
US4543632A (en) | Robust estimation method for determining when subsequent data processing can include sign-bit representations of full-waveform seismic traces | |
US4933911A (en) | Method for determining seismic velocities | |
CA1163353A (en) | Method of reducing monochromatic interference in continuous wave seismic prospecting | |
RU2169382C1 (en) | Method of multilevel vibration-seismic prospecting | |
SU1056100A1 (en) | Vibro-seismic prospecting method | |
RU2122220C1 (en) | Process of seismic prospecting | |
RU2686514C1 (en) | Method for seismic micro-zoning | |
Sivaji et al. | Relationship between fluctuations of arrival time and energy of seismic waves and scale length of heterogeneity: an inference from experimental study | |
GB2190746A (en) | Improvements in or relating to methods of collecting data using seismic vibrators | |
RU2650718C1 (en) | Method of vibration seismic survey | |
Hoover et al. | Vibrator signals | |
Carrion et al. | Direct estimation of shear-wave interval velocities from seismic data | |
Naskar | Effectiveness of stacking in filtering noise from surface wave data | |
RU2809938C1 (en) | Method for vibration seismic exploration | |
Gritto et al. | Seismic Modelling of Seam Waves Excited by Energy Transmission Into a SEAM1 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20050530 |