Цитата:
|
Тогда непонятен dualrail Infiniband с FCA :)
|
У нас основная задача это prestack depth migration и inversion, для изотропных и анизотропных сред и разной acquisition geometry (я не знаю правильной русской терминологии здесь).
В самом упрощенном случае? сейсмический эксперимент порождает для каждого источника множество сейсмических трасс (traces), на которых записываются отраженные волны, порожденные скачками скорости в исследуемом объеме. Например, для случая плоской границы между двумя средами, отраженная волна будет иметь на сейсмограмме вид гиперболы в осях (time, offset (расстояние между источником и приемником)). Очевидно, что в случае сложной структуры внутри исследуемого объема, вид сейсмограмм будет тоже сложным.
Задача сейсмической миграции заключается в том, чтобы правильно переместить амплитуды с сейсмограмм в глубь среды (т.е. из координат time, offset в координаты x,z в простейшем двумерном случае. Если здесь все сделать правильно, то всплески на сейсмических трассах (с определенной степенью коррекции) лягут в правильные позиции, которые соответствуют разрывам скорости в исследуемом объеме, т.е. границами раздела, например, между солью или гранитом и осадочными породами.
В наипростейшем случае высокочастотного (лучевого) приближения, для выполнения миграции мы должны соединить каждую точку в исследуемом (двумерном или трехмерном) объеме c соответствующим источником и приемником лучами, отвечающими падающей отраженной волне (лучи будут нормальны к волновому фронту и строятся при помощи решения определенной системы дифференциальных уравнений). Затем мы должны вычислить суммарное время (traveltime) прохода волны через среду по составному пути (падающий + отраженный лучи) и, взяв с соответствующей трассы амплитуду, отвечающую этому traveltime, поместить ее в соответствующую позицию в глубине объема где строится миграция. На практике, таблицы traveltime вычисляются заранее, и независимо для разных подобъемов мигрируемого объема, и миграция сводится к накоплению амплитуд снимаемых с сейсмических трасс и их коррекции с учетом геометрического расхождения и т.п.
Чтобы алгоритмы сейсмической миграции работали, необходимо иметь не только набор трасс, но и хотя бы приблизительную модель среды, в которой желательно выполнить миграцию. Разумеется, этого у нас нет изначально. Поэтому миграция выполняется последовательно на множестве моделей, которые постепенно приближаются к истиной. Здесь происходит минимизация некоторого функционала, либо регулярными либо стохастическими методами. Это процесс называется seismic inversion, так как это фактически решение обратной задачи (замечу, что это изначально так называемая некорректная задача). Т.е. мы должны выполнять миграцию дикого количества трасс в объеме многих сотен км^3 и делать это сотни раз, чтобы построить адекватную модель.
В конечном результате мы должны получить правильно позиционированные границы раздела между различными геологическими слоями и значения скорости, плотности, коэффициентов анизотропии и т.д. в каждом слое. Это позволяет увидеть есть что-либо интересное в данном объеме или нет.
В действительности алгоритмы существенно более сложные, но общая идея прежняя, даже для reverse time migration. Очень много исходных трасс, сейчас используется многократное перекрытие (т.е. много источников) и длинные offset, а не миграция объема вертикальными полосами. Мигрируемый объем разбивается на сравнительно независимые части, которые легко можно раздать по нодам. Но эти ноды очень много вынуждены обращаться к трассам (т.е. сейсмическим данным, а это часто многие террабайты) и здесь чем быстрее обмен, тем лучше.
Что-нибудь получилось объяснить или нет? Нормально здесь очень долго описывать, я это сейчас писал, говоря по телефону, т.е. мог написать не очень понятно
Цитата:
|
Да и по североевропейской NorduGrid можно набрать вычислителей, не строя один большой кластер. ИМХО
|
У нас будут не только выполнятся тесты нового ПО и исследования новых алгоритмом, но и обрабатываться данные, полученные в результате геологоразведки. Это уже чисто "денежная" деятельность, которая, кроме того, связана с коммерческой тайной, т.е. здесь добровольно-общественные вычислительные ресурсы не будет возможно использовать.
|
