Процедура RADON
Преобразование Радона
Вычисляется прямое или обратное
пребразование Радона или подавление многократных отражений с помощью параболического
преобразования Радона.
Эффективное
подавление кратных волн осуществляется параболическим преобразованием Радона по
сейсмограммам с введенной кинематикой для однократных целевых отражений. В
прямом преобразовании Радона по таким сейсмограммам однократные отражения
сосредоточены вблизи точки p=0. Кратные отражения отображаются на больших значениях параметра p.
По прямому
преобразованию Радона можно определить произвольный
фильтр подавления (выделения) отдельных отражений. Этот фильтр применяется
процедурой RADON_F.
См. также: Методические
рекомендации по применению преобразованияРадона.
Параметры:
Start SP - Начальный пикет SP в метрах.
Значение должно быть умножено на делитель расстояний, указанный в Geom2D.
Stop SP - Конечный пикет SP в метрах.
Значение должно быть умножено на делитель расстояний, указанный в Geom2D.
Процедура будет обрабатывать только сейсмограммы в интервале
[Start SP - Stop SP].
Option
=1 прямое преобразование Радона.
=2 обратное преобразование Радона.
=3 прямое+обратное преобразование Радона.
=4 прямое+обратное преобразование с подавлением кратных волн .
=5 выделение кратных волн.
Mode - Способ преобразования.
1 - Параболическое преобразование g(x)=offset**2.
2 - Псевдо гиперболическое преобразование Фостера/Мошера
g(x) = sqrt(DepthRef**2 + offset**2)
3 - Абсолютное линейное tau-p: g(x) = abs(offset)
4 - Знаковое параболическое преобразование
g(x)=offset**2*sign(offset).
5 - Знаковое гиперболическое преобразование Фостера/Мошера
g(x) = sqrt(DepthRef**2 + offset**2)sign(offset)
6 - Линейное tau-p: g(x) = offset.
Знаковые Преобразования 4,5,6 модно применять при любой сортировке сейсмотрасс. Для отрицательных удалений будут отрицательные значения поправок.
Беззнаковые преобразования 1,2,3 можно применять только при сортировке сейсмограмм по общей глубинной точке. В этом случае сейсмограмма на выходе получится симметричной относительно удаления Intersept Offset.
Левосторонние и правосторонние трассы смешаются.
Offset Ref - Удаление, для которого заданы параметры преобразования.
Значение должно быть умножено на делитель расстояний, указанный в Geom2D.
Intersept Offset - Удаление, с которым ассоциируется параметр tau.
Обычно это нулевое удаление.
Значение должно быть умножено на делитель расстояний, указанный в Geom2D.
Pmin - Минимальная кинематическая поправка на удалении Offset Ref.
Значение должно быть умножено на делитель времен, указанный в Geom2D.
Pmax - Максимальная кинематическая поправка на удалении Offset Ref.
Значение должно быть умножено на делитель времен, указанный в Geom2D.
Все волны с наклонами за пределами [Pmin,Pmax] будут уничтожены.
Отрицательные наклоны присущи отрицательным удалениям.
Должно быть Pmax > Pmin.
Во многих случаях достаточно просто прямого и обратного преобразования для подавления волн-помех.
DP - Шаг перебора кинематических поправок на удалении Offset Ref.
Значение должно быть умножено на делитель времен, указанный в Geom2D.
Pmult1 - Минимальная кинематическая поправка для кратных волн
на максимальном времени трассы для удаления Offset Ref.
Значение должно быть умножено на делитель времен, указанный в Geom2D.
Pmult2 - Минимальная кинематическая поправка для кратных волн.
на минимальном времени трассы для удаления Offset Ref.
Значение должно быть умножено на делитель времен, указанный в Geom2D.
Параметры Pmult1 и Pmult2 применяются только при Option>3.
Будут подавлены все волны с параметрами, превышающими Pmult.
Значение Pmult на каждом времени определяется линейной интерполяцией от Pmult1 до Pmult2.
Должно быть Pmult1 < Pmult2
Depth Ref - Ссылочная глубина для гиперболического преобразования
Фостера/Мошера в метрах.
Nwindows - Number of
windows to use through the mute zone
Nwindows notes:
The parabolic transform runs with higher resolution if the
mute zone is honored. When "nwin" is specified larger than
one (say 6), then multiple windows are used
through the mute
zone. It is assumed in this case that the input data are
sorted by the offkey header item from small offset to large
offset. This causes the code to run 6 times longer. The
mute time is taken from the "muts" header word. Beware,
the SU mute module does not set this header word as one
would normally expect. You have to manually set it yourself.
Freq1 - Первая частота среза низкочастотного фильтра в гц
Freq2 - Вторая частота среза низкочастотного фильтра в гц
Prewhite - Фактор предварительного отбеливания в процентах (0.1).
Может быть задано целым и действительным числом.
Следующие параметры ассоциируются с проблемой пространственного элайсинга при линейном tau-p преобразовании. Не применяйте их при иных преобразованиях.
Ninterp - Параметр интерполяции трасс перед преобразованием.
=0 - без интерполяции
=1 - шаг по х сгущается в 2 раза.
=3 - шаг по х сгущается в 3 раза. и т.д
F1=3.0
low-end frequency in Hz for picking (good default: 3 Hz)
(Known bug: F1 cannot be zero)
F2=20.0
high-end frequency in Hz for picking (good default: 20 Hz)
LenAGC=400
length of AGC operator for picking (good default: 400 ms)
LenT=5
length of time smoother in samples for picker
(good default: 5 samples)
LenX=1
length of space smoother in samples for picker
(good default: 1 sample)
Xopt
1 = use differences for
spatial derivative (works with irregular spacing)
0 = use FFT derivative for spatial
derivatives
(more accurate but requires
regular spacing and
at least 16 input tracs--will
switch to differences
automatically if have less
than 16 input traces
В основу алгоритма положено
Cohen, J. K. and Stockwell, Jr. J. W., (1998), CWP/SU: Seismic Unix
Release 31:
A free package for seismic
research and processing,
Center for Wave Phenomena,
Colorado School of Mines.
Anderson, J. E., 1993, Parabolic and linear 2-D, tau-p transforms
using the generalized radon
tranform, in May 11-14, 1993
Project Review, Consortium Project
on Seismic Inverse methods
for Complex Structures, CWP-137,
Center for Wave Phenomena
internal report.