Способ обработки морских сейсмических данных и способ сейсмической разведки

005475 Область техники Настоящее изобретение относится к способу обработки морских сейсмических данных, и в частности к способу обработки морских сейсмических данных, который включает в себя распространение двух мод сейсмической энергии, одна из которых возникает из-за частичного преобразования мод на границе слоя земли. Способ по настоящему изобретению предусматривает оценку статических сдвигов, вызываемых поверхностным или подповерхностным слоем земли, который имеет сейсмические свойства, отличающиеся от свойств нижележащих слоев. Настоящее изобретение также относится к способу сейсмической разведки, который включает в себя обработку собранных сейсмических данных указанным выше образом. Предшествующий уровень техники На фиг. 1 представлена схематически система для морской сейсмической разведки донного грунта. В этой системе разведки сейсмическая энергия испускается сейсмическим источником 1, который свешен под поверхностью моря с буксировочного судна 2. Сейсмическая энергия испускается вниз, проходит во внутреннее пространство земли и отражается геологической неоднородностью, которая действует в качестве отражателя 3. Отраженная сейсмическая энергия проходит вверх через внутреннее пространство земли в море и детектируется сейсмическим приемником 4, расположенным на морском дне. Информацию о внутреннем пространстве земли получают путем определения времени прохождения сейсмической энергии от сейсмического источника 1 до приемника 4. На фиг. 1 изображен только один отражатель 3, но на практике в качестве частичных отражателей сейсмической энергии будет действовать некоторое количество геологических неоднородностей во внутреннем пространстве земли. На фиг. 1 изображены только один источник 1 и один приемник 4, но на практике система морской сейсмической разведки может иметь ряд источников и ряд приемников. Геологическая структура земли не является однородной. Одна из проблем при обработке морских сейсмических данных заключается в том, что часто существует слой 5 на поверхности или вблизи нее,свойства которого могут значительно отличаться от свойств лежащей под ним геологической структуры 6 (далее упоминается как фундамент). Это может иметь место, например, в том случае, если существует некоторый слой на поверхности земли или вблизи нее, который является менее консолидированным,чем фундамент. В частности, скорость сейсмической энергии может быть значительно более низкой в поверхностном или подповерхностном слое 5, чем в фундаменте 6, и такой поверхностный или подповерхностный слой повсеместно известен в качестве низкоскоростного слоя (или LVL от англ. lowvelocity layer). Эта разность в скорости будет вызывать сдвиг во времени прохождения сейсмической энергии по сравнению со временем прохождения, которое регистрировалось бы, если бы поверхностный или подповерхностный слой и фундамент имели идентичные сейсмические свойства, и эти сдвиги во времени прохождения известны как статические сдвиги или просто статика. Низкоскоростной слой 5 изображен на фиг. 1 как поверхностный слой, но он не должен обязательно простираться до поверхности, и поэтому может существовать дополнительный слой, лежащий выше низкоскоростного слоя. Статический сдвиг, генерируемый поверхностным или подповерхностным низкоскоростным слоем 5, зависит от толщины этого слоя и от скорости распространения сейсмической энергии через этот слой. Как правило, существуют горизонтальные вариации как по толщине низкоскоростного слоя 5, так и по скорости распространения через этот слой, так что статический сдвиг, наблюдаемый на сейсмическом приемнике в одном положении, вероятно, будет отличным от статического сдвига, наблюдаемого на приемнике в другом положении. В качестве первого приближения, все множество данных, регистрируемое на одном приемнике, будет иметь опережающий или опаздывающий статический временной сдвиг по отношению к данным, регистрируемым на другом приемнике. В высшей степени желательно учесть статический сдвиг при обработке сейсмических данных. До тех пор пока эти статические сдвиги не удалены из сейсмических данных, будет существовать неоднозначность относительно того, вызываются ли вариации времени поступления сейсмических событий от более глубоких слоев вариациями по глубине или по горизонтальным положениям этих более глубоких слоев или просто возникают из-за эффектов распространения в низкоскоростном слое 5. Сущность изобретения В своем первом аспекте настоящее изобретение предлагает способ обработки сейсмических данных, содержащих первую и вторую моды сейсмической энергии, причем первая и вторая моды сейсмической энергии распространялись через слой внутреннего пространства земли вдоль, по существу, одного и того же пути, включающий в себя стадию обработки сейсмических данных для получения разности времен прохождения через слой внутреннего пространства земли между сейсмической энергией, распространяющейся в первой моде, и сейсмической энергией, распространяющейся во второй моде. В предпочтительном варианте воплощения вторая мода генерируется путем частичного преобразования первой моды на границе раздела слоя донного грунта. На фиг. 1 показано, что сейсмическая энергия проходит через низкоскоростной слой 5 2 раза: один раз - когда она проходит от сейсмического источника 1 к отражателю 3, и снова - когда она проходит от отражателя 3 к приемнику 4. Оба прохождения низкоскоростного слоя 5 будут вызывать статические-1 005475 сдвиги; при этом сдвиг, вызываемый при прохождении распространяющейся вниз сейсмической энергии через низкоскоростной слой 5, известен как статический сдвиг со стороны источника, а сдвиг, вызываемый при прохождении распространяющейся вверх энергии через низкоскоростной слой 5, известен как статический сдвиг со стороны приемника. Настоящее изобретение предусматривает способ оценки статического сдвига со стороны приемника. Способ согласно настоящему изобретению может быть использован для обработки уже имеющихся сейсмических данных. Альтернативно, он может быть включен в способ сейсмической разведки для обработки полученных данных по мере их сбора или впоследствии. Согласно одному варианту воплощения настоящего изобретения предполагается, что контраст в физических свойствах между фундаментом 6 и низкоскоростным слоем 5 является достаточно большим,чтобы имело место значительное преобразование между р-модой распространения энергии и s-модой распространения энергии, или наоборот, при распространении сейсмической энергии вверх через границу раздела между фундаментом 6 и низкоскоростным слоем 5. В результате такого преобразования мод для каждого события р-моды, регистрируемого с помощью сейсмического приемника выше границы раздела между низкоскоростным слоем и фундаментом, будет существовать соответствующее событие sмоды. Событие р-моды и соответствующее событие s-моды будут происходить в различные моменты времени, поскольку скорости распространения р-моды и s-моды в низкоскоростном слое не равны друг другу. Однако, поскольку преобразование мод происходит на нижней границе низкоскоростного слоя, в результате различия скоростей двух мод в низкоскоростном слое должно возникать время задержки между событием р-моды и соответствующим событием s-моды. Время задержки между событием р-моды и соответствующим событием s-моды не будет подвергаться значительному влиянию со стороны свойств фундамента. Разность времен прохождения между р- и s-модами через низкоскоростной слой может быть легко определена путем идентификации события р-моды и соответствующего события s-моды в сейсмических данных и определения времени задержки между событием р-моды и соответствующим событием s-моды. Преобразование мод также может происходить тогда, когда распространяющаяся вниз волна претерпевает отражение на границе раздела между низкоскоростным слоем 5 и фундаментом 6. Настоящее изобретение может быть применено и в этом случае, поскольку, когда при отражении происходит преобразование мод, отраженный сигнал будет содержать р-компоненту и s-компоненту, и эти две компоненты будут иметь различные времена прохождения через низкоскоростной слой 5. Преобразование мод также может происходить, когда распространяющаяся вниз волна претерпевает критическое преломление на границе раздела между низкоскоростным слоем 5 и фундаментом 6, с генерацией сейсмической волны, распространяющейся вдоль границы раздела между низкоскоростным слоем 5 и фундаментом 6. Распространяющаяся вдоль границы раздела волна будет возбуждать в низкоскоростном слое волны как р-моды, так и s-моды, и настоящее изобретение может быть применено к этим р- и s-модам. Для любого параметра, являющегося показателем одного из аспектов сейсмических данных, таких как давление или компонента движения частицы (термин движение частицы включает в себя смещение частицы, скорость частицы, ускорение частицы и высшие производные от смещения частицы), события одной из р- и s-мод, как правило, будут в целом проявляться сильнее, чем события любой другой моды. Следовательно, в соответствии с предпочтительными вариантами воплощения настоящего изобретения для локализации события р-моды и соответствующего события s-моды используются два параметра, являющихся показателями двух различных аспектов сейсмических данных. В принципе, если на нижней границе низкоскоростного слоя происходит достаточно сильное преобразование мод, то соответствующие пары событий р-моды и s-моды могут быть локализованы путем анализа полученных сейсмических данных. Однако во многих случаях амплитуда событий, связанных с преобразованием мод, является низкой, и тогда предпочтительна методика взаимной корреляции или деконволюции (т.е. восстановления сигнала методом обращения свертки). Статический сдвиг для р-волны обычно является малым по сравнению со статическим сдвигом дляs-волны, поскольку скорость р-волны больше, чем скорость s-волны, и, следовательно, вариация во времени прохождения меньше. Кроме того, менее консолидированные осадки (отложения) донного грунта,находящиеся в низкоскоростном слое 5, обычно не вызывают настолько же значительных горизонтальных вариаций в скорости р-волн по сравнению с тем, какие горизонтальные вариации они вызывают в скорости s-волн. Во многих случаях статический сдвиг р-волны настолько мал, что с большой точностью он может рассматриваться как пренебрежимо малый. Настоящее изобретение позволяет оценить разность между статическим сдвигом s-волны и статическим сдвигом р-волны в каждом приемнике. Если статика р-волны известна, т.е. может быть оценена,или можно предположить, что она равна нулю, то настоящее изобретение позволяет получить статику sволны. Согласно второму аспекту настоящего изобретения предлагается способ сейсмической разведки,включающий в себя стадии направления сейсмической энергии, распространяющейся в первой моде, к границе слоя донного грунта таким образом, что на указанной границе происходит частичное преобразо-2 005475 вание моды сейсмической энергии во вторую моду; сбора сейсмических данных, содержащих первую и вторую моды сейсмической энергии, на одном или нескольких приемниках; и обработки сейсмических данных в соответствии с определенным выше способом первого аспекта изобретения. Согласно третьему аспекту настоящего изобретения предложено устройство для обработки сейсмических данных, включающих в себя первую и вторую моды сейсмической энергии, причем вторая мода сгенерирована за счет частичного преобразования первой моды на границе слоя донного грунта, в результате чего первая и вторая моды сейсмической энергии распространялись через слой донного грунта вдоль, по существу, одного и того же пути, содержащее средство для обработки сейсмических данных с получением разности времен прохождения через указанный слой между сейсмической энергией, распространяющейся в первой моде, и сейсмической энергией, распространяющейся во второй моде. В предпочтительном варианте воплощения устройство содержит программируемый процессор для обработки данных. Краткое описание чертежей Предпочтительные варианты воплощения настоящего изобретения описаны ниже со ссылками на прилагаемые чертежи, где фиг. 1 представляет собой схему морской сейсмической разведки; фиг. 2 представляет собой схематическую иллюстрацию преобразования моды распространяющейся вверх р-волны на границе раздела между фундаментом и низкоскоростным слоем; фиг. 3 представляет собой схематические примеры трасс сейсмических данных для ситуации согласно фиг. 2; фиг. 4 представляет собой кросс-коррелограмму для трасс сейсмических данных согласно фиг. 3; фиг. 5(а) и 5(b) демонстрируют искусственные сейсмические данные, генерируемые по модели согласно прилагаемой ниже таблице; фиг. 6(а) и 6(b) демонстрируют деконволюцию данных согласно фиг. 5(а) и 5(b); фиг. 7 демонстрирует теоретические и экспериментальные результаты для разности времен прохождения pp-ps; фиг. 8 представляет собой схематическую иллюстрацию системы сейсмической разведки, имеющей геометрию источника и приемника с общим смещением; фиг. 9(а) и 9(b) демонстрируют вертикальную и радиальную компоненты скорости частицы, полученные с использованием системы сейсмической разведки согласно фиг. 8; фиг. 10(а) и 10(b) демонстрируют вертикальную и радиальную компоненты скорости частицы по фиг. 9(а) и 9(b) после предварительной обработки; фиг. 11 представляет собой схематическую иллюстрацию преобразования мод распространяющейся вверх s-волны на границе раздела между фундаментом и низкоскоростным слоем; фиг. 12(а) представляет собой схематическую иллюстрацию преобразования мод распространяющейся вниз р-волны при отражении на границе раздела между низкоскоростным слоем и фундаментом; фиг. 12(b) представляет собой схематическую иллюстрацию преобразования мод распространяющейся вниз s-волны при отражении на границе раздела между низкоскоростным слоем и фундаментом согласно изобретению; фиг. 13(а) представляет собой схематическую иллюстрацию преобразования мод р-волны, распространяющейся вдоль границы раздела между низкоскоростным слоем и фундаментом; фиг. 13(b) представляет собой схематическую иллюстрацию преобразования мод s-волны, распространяющейся вдоль границы раздела между низкоскоростным слоем и фундаментом; фиг. 14 представляет собой блок-схему устройства для обработки морских сейсмических данных согласно изобретению. Описание предпочтительных вариантов воплощения изобретения Принципы настоящего изобретения будут далее описаны со ссылкой на тот вариант воплощения, в котором за счет частичного преобразования мод генерируется пара соответствующих событий. Фиг. 2 иллюстрирует схему преобразования мод при прохождении распространяющейся вверх р-волны через границу раздела между фундаментом 6 и низкоскоростным слоем 5. Распространяющаяся вверх р-волна 7 частично проходит как р-волна 7', но также претерпевает частичное преобразование в s-волну 7". Предполагается, что распространяющаяся вверх р-волна была сгенерирована путем отражения от глубокой границы раздела между двумя слоями внутри земли или путем преломления через более глубокие слои. Преломление также будет происходить на границе раздела между фундаментом 6 и низкоскоростным слоем 5, и угол преломления для проходящей р-волны 7' будет отличаться от угла преломления для преобразованной s-волны 7". Проходящая р-волна 7' и преобразованная s-волна 7" вместе будут давать сигнал в ряду приемников. Проходящая р-волна 7' и преобразованная s-волна 7" распространяются, как показано на фиг. 2, в различных направлениях и, следовательно, будут приходить в приемники, расположенные в разных местах, а не в один и тот же приемник. Однако подобная же преобразованная s-волна от соседней точки преобразования будет приходить в тот же самый приемник, что и проходящая р-волна 7'. Предположим, что-3 005475 для проходящей р-волны и преобразованной s-волны, полученных в конкретном приемнике, различие между точкой преобразования (у основания низкоскоростного слоя) преобразованной s-волны и точкой в основании низкоскоростного слоя, через которую прошла р-волна, значительным образом не влияет на оценки статических сдвигов, получаемых с помощью способа согласно изобретению. Проходящая р-волна 7 будет регистрироваться преимущественно на вертикальной компоненте сейсмической энергии, регистрируемой в приемнике, а преобразованная s-волна 7" будет регистрироваться преимущественно на радиальной компоненте сейсмической энергии, регистрируемой в приемнике 4 (радиальное направление представляет собой направление между источником и приемником, проектируемое на донный грунт, и это направление будет обозначено как х-направление). Это происходит вследствие того, что в практической системе сейсмической разведки волны, которые распространяются вверх из глубины земли, имеют тенденцию составлять угол 30 или менее относительно вертикали для типичной структуры внутреннего пространства земли и типичной глубины целевого отражателя. Волна сдвига регистрируется преимущественно на радиальной компоненте, поскольку движение частицы для волны сдвига перпендикулярно к направлению распространения, в то время как движение частицы для р-волны направлено вдоль направления распространения. На фиг. 3 представлены схематические иллюстрации двух параметров, являющихся показателями аспектов сейсмических данных, регистрируемых в приемнике 4 в том случае, когда происходит частичное преобразование мод (фиг. 2). Эти два параметра представляют собой измеряемую приемником радиальную компоненту (х-компоненту) скорости частицы и измеряемую приемником вертикальную компоненту (z-компоненту) скорости частицы. Поскольку s-волна 7" образуется при преобразовании распространяющейся вверх р-волны 7, то в том случае, если граница раздела между фундаментом и низкоскоростным слоем локально является приблизительно планарной, регистрируемая s-волна будет содержать импульс (сейсмический импульсоид), подобный поступившему р-сигналу. Поступление s-волны будет наблюдаться не в то же самое время, что и поступление р-волны, а через некоторое время задержки dt(по сравнению с поступлением р-волны), вызванное различием между скоростью р-волн и скоростью sволн в низкоскоростном слое 5. Время задержки dt совпадает с разницей между статическим сдвигом рволны и статическим сдвигом s-волны. Таким образом, если статический сдвиг р-волны известен или может быть точно оценен, статический сдвиг s-волны может быть получен путем простого добавления времени задержки dt к статическому сдвигу р-волны. Время задержки dt может быть получено с помощью любого способа сравнения времени поступления р-волн со временем поступления s-волн. Однако в предпочтительных вариантах воплощения изобретения время задержки dt определяется либо путем деконволюции (восстановления методом обращения свертки), либо путем взаимной корреляции (вычисления функции взаимной корреляции) вертикальной компоненты и горизонтальной компоненты сейсмической энергии, регистрируемых сейсмическим приемником. На фиг. 4 представлена кросс-коррелограмма, полученная путем взаимной корреляции горизонтальной и вертикальной компонент скорости частицы, изображенных на фиг. 3. Видно, что существует пик с нулевым временным сдвигом, который возникает ввиду того, что поступившая р-волна 7' имеет ненулевую амплитуду х-компоненты измеренной скорости частицы (хотя она имеет место преимущественно в z-компоненте измеренной скорости частицы). Любая энергия s-волны, которая проявляется в вертикальной компоненте измеренной скорости частицы, будет также вызывать появление пика с нулевым временным сдвигом на кросс-коррелограмме. Пик в момент времени dt на кросс-коррелограмме возникает от пика поступившей р-волны в вертикальной компоненте измеренной скорости частицы и пика поступившей s-волны в х-компоненте измеренной скорости частицы. Для точного разрешения этого пика и пика с нулевым временным сдвигом желательно (но необязательно), чтобы данные пики не перекрывались. Пики не будут перекрываться,если импульс на регистрируемых данных имеет достаточно высокую частоту для того, чтобы толщина низкоскоростного слоя 5 была больше, чем две длины импульса. Если невозможно быть уверенным, что толщина низкоскоростного слоя 5 является достаточно большой, чтобы предотвратить перекрывание пиков, точность определения временного сдвига dt может быть увеличена путем деконволюции хкомпоненты и z-компоненты измеренной скорости частицы. Альтернативно, регистрируемое волновое поле может быть полностью разложено на его р-волновую компоненту и его s-волновую компоненту,например, с использованием методики, которая раскрыта в заявках на патент Великобритании 0001355.7 или 0003406.6, или публикации Робертсона и др. (Robertsson et al.), "Wavefield separation and"Wavefield separation and estimation of near surface properties in land seismic", 62nd EAGE Conference Glasgow, Extended Abstracts, 2000. Такие методики разложения имеют эффект удаления утекающих пиков,изображенных на фиг. 3, из х-компоненты и z-компоненты измеренной скорости частицы. Если полученная разложением р-компонента и полученная разложением s-компонента регистрируемого волнового поля затем взаимно коррелируются (т.е. вычисляется функция их взаимной корреляции), пик с нулевым-4 005475 временем задержки существовать не должен, в то время как пик с временем задержки dt для p-s должен увеличиться по сравнению с тем, который изображен на фиг. 4. Если сейсмические волны распространяются преимущественно в радиальном направлении (вместо того, например, чтобы отражаться назад в отрицательном радиальном направлении), то взаимная корреляция или деконволюция вертикальной и радиальной компонент, как ожидается, будет давать положительный пик с временем задержки dt для p-s, поскольку скорость сейсмической энергии, как правило,возрастает по направлению вниз по всей границе раздела между низкоскоростным слоем 5 и фундаментом 6. Взаимная корреляция или деконволюция между вертикальной и поперечной компонентами регистрируемой сейсмической энергии должна детектировать распространение волны, которая не находится в радиально-вертикальной плоскости, например, вызываемое подповерхностным рассеянием. Согласно предпочтительному варианту воплощения кросс-коррелограммы, получаемые приемником от множества различных источников сейсмических волн, усредняются для увеличения отношения сигнала к шуму на итоговой кросс-коррелограмме, известной как суммарная кросс-коррелограмма. Это должно увеличить отношение между амплитудой пика с временным сдвигом на вертикально-радиальной кросс-коррелограмме или на диаграмме деконволюции к амплитуде на вертикально-поперечной кросскоррелограмме или диаграмме деконволюции, если граница раздела между низкоскоростным слоем и фундаментом локально является приблизительно горизонтальной, а низкоскоростной слой является приблизительно изотропным азимутально. Настоящее изобретение не ограничивается использованием методик деконволюции или взаимной корреляции для получения разности времен прохождения через низкоскоростной слой, и поэтому может быть использована любая методика. В принципе, может быть использован любой алгоритм, который устанавливает корреляцию между двумя трассами. Еще один из примеров пригодной для использования методики представляет собой метод оценки времени задержки бикогерентности, описанный Икелле (L.Ikelle), "Geophysics", том 62, стр. 1947, 1997 г., с использованием радиальной и вертикальной компонент скорости частицы (или, как будет обсуждаться ниже, давления и радиальной компоненты скорости частицы). На фиг. 2 представлено преобразование мод, происходящее для распространяющейся вверх рволны. Подобный же эффект наблюдается, если распространяющаяся вверх волна представляет собой sволну (фиг. 11). Когда распространяющаяся вверх s-волна 8 проходит через границу между фундаментом 6 и низкоскоростным слоем 5, она претерпевает частичное прохождение как s-волна 8' и частичное преобразование в р-волну 8". Как преобразованная р-волна 8", так и проходящая s-волна 8' затем распространяются через низкоскоростной слой 5 к приемникам 4. Следовательно, единственное различие между этим случаем и случаем, представленным на фиг. 2, заключается в том, что амплитуда поступающейs-волны будет теперь больше, чем амплитуда поступающей р-волны. Время задержки между поступлением р-волны и соответствующим поступлением s-волны по-прежнему будет составлять dt, как в примере на фиг. 2 и 3. Следовательно, взаимная корреляция или деконволюция вертикальной и горизонтальной компонент сейсмических данных должны опять демонстрировать пик с временным сдвигом для p-s (то есть для разности между временем прохождения р-волн через низкоскоростной слой и временем прохождения s-волн через низкоскоростной слой 5). Распространяющаяся вверх сдвиговая волна, как правило, будет преломляться в направлении нормали к границе раздела. Однако, учитывая обычное уменьшение скоростей s-моды при входе в низкоскоростной слой 5 снизу, преобразованная s-p волна будет преломляться в сторону от указанной нормали. Это должно быть принято во внимание при выборе диапазона смещений, в котором применяются методики деконволюции или взаимной корреляции. Один из видов применения настоящего изобретения заключается в извлечении разности между статическим сдвигом р-волны и статическим сдвигом s-волны из морских сейсмических данных, регистрируемых с использованием расположенного на донном грунте приемника. В этом случае статический сдвиг р-волны настолько мал, что его можно игнорировать так, что статический сдвиг p-s, извлекаемый из этих данных с помощью заявленного способа, приблизительно равен статическому сдвигу s-волны. Для иллюстрации способа согласно настоящему изобретению искусственные сейсмические данные генерировались для простой одномерной модели с использованием модели отражения того типа, который описан Кеннетом (B.L.H. Kennett) в "Seismic Wave Propagation in Stratified Media", Cambridge University Press, Cambridge, England, 1983. Детали используемой модели представлены в таблице, и из нее можно легко понять, что модель состоит из низкоскоростного слоя, имеющего толщину 100 м и расположенного между слоем воды и слоем фундамента. Слой воды предполагается имеющим бесконечную глубину, и слой фундамента также предполагается имеющим бесконечную глубину. Взрывной точечный источник сейсмической энергии, расположенный глубже в земле (представляющий отражения от более глубоких слоев), используется как источник сейсмической энергии и располагается на 200 м ниже границы раздела между низкоскоростным слоем и слоем фундамента. Линейный ряд приемников расположен на донном грунте. Слой Вода Низкоскоростной слой Фундамент Глубина расположения источника: 200 м ниже границы раздела низкоскоростной слой/фундамент. Тип источника: взрывной точечный источник, генерирующий только р-волны, 100 Гц импульсоид Рикера. Положение приемника: на границе раздела вода/низкоскоростной слой, смещение 0-100 м с промежутками 5 м. Искусственные данные, полученные из этой модели, изображены на фиг. 5(а) и 5(b). На фиг. 5(а) представлена вертикальная компонента скорости частицы, а на фиг. 5(b) - радиальная компонента скорости частицы. Предполагается, что источник генерирует р-волны. Сигнал, регистрируемый приемниками, содержит проходящую р-компоненту и s-компоненту, возникающую из-за преобразования мод, происходящего на границе раздела между фундаментом и низкоскоростным слоем. Поступление проходящей ркомпоненты видно преимущественно на вертикальной компоненте регистрируемых сейсмических данных, показанной на фиг. 5(а). На радиальной компоненте данных, изображенной на фиг. 5(b), также видны энергия проходящей р-волны и поступление преобразованной p-s волны, в особенности, при увеличении смещения. При типичной сейсмической отражательной разведке донного грунта большая часть проходящей рэнергии будет регистрироваться на вертикальной компоненте приемника. Как правило, скорость сейсмической энергии резко понижается при вхождении снизу в низкоскоростной слой, так что пути лучей проходящих волн будут сильно преломляться в направлении вертикали. В противоположность этому, sволна, генерируемая путем преобразования на границе раздела между фундаментом и низкоскоростным слоем, будет проявляться на радиальной компоненте сейсмической энергии, регистрируемой на приемнике. Предположим, что граница раздела между фундаментом 6 и низкоскоростным слоем 5 является горизонтальной, так что во время преобразования р-моды в s-моду на границе раздела волна sh-моды не генерируется. Тангенциальная компонента регистрируемых данных будет равна нулю, и преобразованная s-волна будет представлять собой чистую sv-волну. Тогда теоретический отклик смещения для низкоскоростного слоя может быть записан согласно Лангстону (С.А. Langston, "J. Geophys. Res." том 84,стр. 4749-4762, 1979 г.) в видеDR(t)=S(t)ER(t) В уравнениях (1) S(t) представляет собой временную функцию источника падающей р-волны, Dv(t) и DR(t) представляют собой вертикальную и радиальную компоненты данных, a Ev(t) и ЕR(t) представляют собой соответственно вертикальный и радиальный импульсные отклики. Путем деконволюции вертикальной компоненты регистрируемых данных из радиальной компоненты из этих данных может быть точно извлечена разность времен прохождения между проходящей р-волной и преобразованной sволной. Фундаментальная идея, на которой основывается этот способ, заключается в том, что вертикальная компонента регистрируемых данных может рассматриваться в качестве хорошей оценочного параметра для импульсоида источника, поскольку импульсный отклик низкоскоростного слоя 5 для вертикальной компоненты данных, для входящей плоской р-волны, близок к -функции (см. Langston, выше). Другими словами,Ev(t)(t) или S(t)Dv(t) Реверберации р-волны также будут присутствовать на вертикальной компоненте, но они сильно уменьшены по амплитуде. Результат деконволюции заключается в переходной функции (называемой в сейсмологии землетрясений функцией приемника) низкоскоростного слоя для энергии сдвиговой волны Результат во временной области может быть интерпретирован очень похожим образом в виде ряда отражений. Поскольку деконволюция в своей основе представляет собой разделение в частотной области, оно может стать нестабильным. По этой причине, предпочтительно избегать слишком малых значений в знаменателе. Это достигается путем заполнения спектральных дырок знаменателя до доли с от его максимума (техника предложена S.K. Dey-Sarkar and R.A. Wiggins, "J. Geophys. Res.", т. 8, стр. 3633-3641,1976). Другая проблема, которая может возникнуть, заключается в том, что шум может приводить к появлению в результатах деконволюции высоких частот, которые являются нефизическими и нечетко присутствуют в полученных данных. Это может быть предотвращено путем умножения результата в частотной области на функцию Гаусса с максимумом в нуле и использования параметра ширины кривой Гаус-6 005475 са, т.е. а, для управления его шириной и, таким образом, частотным диапазоном. Стабильная деконволюция в частотной области определяется выражением (см. Langston, выше) где ФSS=mахDD,с.mах[DD], и G=е/4a Здесь DR и Dv представляют собой Фурье-преобразования соответственно радиальной и вертикальной компонент данных, а звездочка обозначает комплексное сопряжение. ЕR представляет собой полученный деконволюцией радиальный отклик земли в частотной области, и он может быть непосредственно преобразован обратно во временную область с получением желаемой функции приемника. ФSS и G соответственно представляют собой описанные стабилизированный знаменатель и Гауссов фильтр. Из-за невертикального падения некоторая часть энергии р-волны будет присутствовать на радиальной компоненте данных (фиг. 5(b. Она будет также присутствовать в полученном деконволюцией результате, где она должна присутствовать точно в нулевой момент времени. Из-за стабилизации деконволюции ширина импульсов функции приемника становится конечной, и по этой причине часть их будет видна вблизи конца трассы функции приемника (фиг. 6(а. В предпочтительном варианте воплощения это предотвращается путем дополнения нулями - на отрезке n трассы - перед радиальной компонентой и путем дополнения таким же количеством нулей после вертикальной компоненты. После этого полученная деконволюцией трасса (длина 2n) будет иметь свой нулевой момент времени на отсчете n+1, и вокруг него не будет появляться никакого паразитного поступления сигнала (фиг. 6(b. В идеале, функция приемника может быть интерпретирована как отклик низкоскоростного слоя на падающую снизу плоскую Р-волну в виде сдвиговой волны. По этой причине амплитуды функций приемника содержат информацию о свойствах среды в низкоскоростном слое (то есть о скорости и плотности) и о различии скоростей в низкоскоростном слое и фундаменте. Однако стабилизация деконволюции путем заполнения спектральных дыр изменяет амплитуду функции приемника. Следуя предложению,сделанному Аммоном (C.J. Ammon, "Bull. Seismol. Soc. Am." том 81, стр. 2504-2510, 1991), могут быть извлечены истинные амплитуды функции приемника. Эффект стабилизации деконволюции может быть оценен путем деконволюции самой вертикальной компоненты (с использованием тех же значений для параметров а и с), зная, что максимальная амплитуда этой полученной деконволюцией трассы должна состоять из единственного пика в нулевой момент времени с единичной амплитудой. По этой причине,путем деления функции приемника на максимальное значение этой автодеконволюции, могут быть извлечены правильные значения амплитуд. На фиг. 7 проиллюстрирована теоретическая разность времен прохождения между проходящей рволной и преобразованной р-волной, вычисленная из параметров модели таблицы, в сравнении с измеренной разностью времен, взятой из искусственных трасс приемника на фиг. 5(а) и 5(b). На фиг. 7 теоретическая разность времен прохождения изображена в виде прерывистой линии, а измеренная разность времен, извлеченная из искусственных данных фиг. 5(а) и 5(b), изображена в виде сплошной линии. Следует отметить, что сплошная линия имеет ступеньки по 0,2 мс, и это происходит из-за того, что при генерировании искусственных сейсмических данных использовался временной интервал получения выборки в 0,2 мс. Можно увидеть, что время задержки, полученное путем деконволюции искусственных данных согласно фиг. 5(а) и 5(b), хорошо согласуется с теоретической разностью времен прохождения,вычисленной из параметров модели. В модели согласно таблице данные генерируются источником, расположенным во внутреннем пространстве земли под нижней границей низкоскоростного слоя 5. Однако методики, описанные выше в связи с моделью согласно таблице, могут быть использованы аналогичным способом для данных, полученных с использованием сейсмической энергии, которая испускалась из поверхностного источника и которая преломилась или отразилась от более глубоких слоев во внутреннем пространстве земли, так что она распространяется вверх сквозь часть модели внутреннего пространства земли, используемой в этом способе. Использование вертикальной и радиальной компонент скорости частицы в качестве двух параметров, являющихся показателями сейсмических данных, как описано выше, имеет то преимущество, что события р-моды и события s-моды разделяются естественным образом, при этом события р-моды проявляются преимущественно на вертикальной компоненте, а события сдвиговой волны проявляются преимущественно на радиальной компоненте. Однако настоящее изобретение не ограничивается использованием вертикальной и радиальной компонент скорости частицы в качестве двух параметров, являющихся показателями сейсмических данных. При деконволюции искусственных данных согласно фиг. 5(а) и 5(b) такая деконволюция осуществлялась для всех смещений с использованием всего временного промежутка искусственных трасс. Это было возможно благодаря исключительной простоте модели земли, для которой вычислялись искусственные данные, и простой геометрии сбора данных (система прохождения источник-приемник, где источник располагается в области фундамента указанной модели). Однако простой прямолинейный подход не может быть использован на практике. При практической разведке сейсмические данные собираются с-7 005475 использованием более сложных геометрий сбора данных, и поэтому деконволюция данных для всех смещений в течение всей длительности трасс потребовалa бы значительной вычислительной мощности. По этой причине, более предпочтительным является осуществление деконволюции или взаимной корреляции для выбранного диапазона времени и/или смещения с целью сокращения требуемой обработки. Также предпочтительно выбирать диапазон времени и/или смещения таким образом, чтобы уменьшать шум в результатах процедуры деконволюции или взаимной корреляции. Как указывалось в изложенном выше объяснении теоретических основ метода деконволюции, полученный путем деконволюции результат представляет собой, по существу, импульсный отклик в виде сдвиговой волны от низкоскоростного слоя при падении на него волны сжатия снизу. По этой причине деконволюция и взаимная корреляция могут рассматриваться как пути для определения некоторого фильтра, представляющего собой низкоскоростной слой, при этом в рассматриваемом случае особенный интерес представляет собой часть фильтра, осуществляющая задержку по фазе. Рассматривая все это,становится ясным, что в выбранном окне данных не должна присутствовать никакая распространяющаяся вниз энергия, такая как реверберации в водном слое и прямая волна от источника до приемника. Более того, в выбранных данных не должно присутствовать никаких волн, преобразованных в фундаменте в сдвиговую энергию (если только преобразование s-в-р на границе низкоскоростного слоя не используется для нахождения разности между р-статическим сдвигом и s-статическим сдвигом, как будет обсуждаться ниже со ссылками на фиг. 11). Эти типы сейсмической энергии не должны фильтроваться низкоскоростным слоем, и по этой причине они могут рассматриваться как шумы для этого способа. Данные,используемые в процессе деконволюции/взаимной корреляции, предпочтительно выбираются для максимально возможного удаления этой энергии. Длительность окна, используемого для выбора данных для процесса деконволюции/взаимной корреляции, предпочтительно является короткой, чтобы минимизировать требуемую обработку. Хотя теоретически низкоскоростным слоем могут быть отфильтрованы несколько последовательных событий распространяющихся вверх отраженных р-волн (то есть они частично преобразуются в сдвиговые волны) точно таким же способом, однако, в области спектральных оценок уже давно было установлено, что выбор более длительного окна данных не уменьшает разброс результата, хотя соотношение фильтра между обоими компонентами является одним и тем же. Кроме того, более длительные окна потенциально содержат больше нежелательных поступлений, что приводит к дополнительным шумам в результатах процесса взаимной корреляции или деконволюции. Временное окно для выбора радиальной компоненты данных должно быть выбрано таким образом, чтобы перекрывать поступление сдвиговой волны, которая имеет время задержки, соответствующее максимальному ожидаемому времени задержки. При выборе части данных для деконволюции/взаимной корреляции также предпочтительно выбрать данные по смещению. Данные с нулевым смещением для плоской слоистой среды не будут содержать никакой энергии сдвиговой волны, возникающей в результате преобразования s-волны в р-волну на границе низкоскоростного слоя, поскольку преобразование не будет иметь места для падающих вертикально волн сжатия. С другой стороны, данные с длинным смещением (то есть данные, у которых смещение существенно больше, чем глубина отражателя) не будут удовлетворять предположению о поверхностно-согласованной статике, и статические поправки станут динамическими поправками. По этой причине предпочтительно исключить данные с нулевым смещением и данные с длинным смещением из тех данных, которые используются для деконволюции/взаимной корреляции. Для такого места разведки, где существует известная информация о структуре низкоскоростного слоя и фундамента, возможно вычисление диапазона смещений, в котором ожидается появление максимального p-s преобразования. Если это так, диапазон смещений, выбранный для деконволюции/взаимной корреляции, должен включать в себя диапазон смещений, в котором ожидается появление максимального p-s преобразования. В другом варианте воплощения вместо вертикальной компоненты скорости частицы используются(масштабированные) данные по давлению. Как замечено Ли и др. (X. Li et al.) в работе "Lomond Data(Applied Seismic Anisotropy), Vol. 7, Converted waves II: Case examples, стр. 185-212, 1999/2000, утечка малых количеств (5%) энергии сдвиговой волны от радиальной компоненты скорости частицы к вертикальной компоненте скорости частицы может происходить из-за нежелательной взаимной корреляции между двумя этими компонентами скорости при использовании кабелей для сбора данных с донного грунта. Такая нежелательная взаимная корреляция может сопровождаться небольшим фазовым сдвигом порядка 6 мс, и по этой причине она будет проявляться в результатах деконволюции и взаимной корреляции в виде статического события с разницей времен прохождения порядка 6 мс. Между давлением и радиальной компонентой скорости частицы таких эффектов взаимной корреляции не существует, и по этой причине функции приемника, вычисляемые путем деконволюции или взаимной корреляции давления и радиальной компоненты скорости частицы, не страдают от этого фазового сдвига. Однако из-за скалярной природы регистрации давления соответствующая нулевому времени или проекционная часть результатов деконволюции или взаимной корреляции будет более значительной.-8 005475 В описанных выше вариантах воплощения временной сдвиг dt определяется путем деконволюции или взаимной корреляции параметров, являющихся показателями измеренных сейсмических данных. Эти методики являются преимущественными тогда, когда события преобразования мод имеют низкую амплитуду. Если разность свойств между низкоскоростным слоем 5 и фундаментом 6 является большей,чем та, которая является достаточно сильной для того, чтобы работал способ, события p-s преобразования могут иметь столь высокие амплитуды, что не будет необходимым осуществлять деконволюцию или взаимную корреляцию (то есть в явном виде вычислять функции приемника), по которым могут быть оценены разности времен прохождения. Вместо этого, изменение фазы при p-s преобразовании может наблюдаться непосредственно в сейсмических данных. На фиг. 8 представлена система сейсмической разведки, в которой данные собираются с помощью пар источник-приемник, все из которых имеют одинаковое смещение, т.е. смещение является равным,как по величине, так и по знаку, для каждой пары источник-приемник. Эта конфигурация известна как геометрия сбора с общим профилем смещения (по знаку). Источники испускают сейсмическую энергию в р-моде, которая отражается отражателем 3 назад к приемникам 4. Преобразование мод происходит при прохождении отраженной р-волны вверх через границу раздела между фундаментом 6 и низкоскоростным слоем 5, так что принимаемая на приемниках сейсмическая энергия содержит как р-моды, так и sмоды. Смещение пар источников и приемников выбирается таким образом, чтобы соответствовать диапазону, в котором ожидается появление максимального p-s преобразования (смещение не является нулевым смещением и не является слишком большим). На фиг. 9(а) и 9(b) представлены соответственно вертикальная компонента и радиальная компонента скорости частицы, измеренные с помощью приемников 4 в системе разведки согласно фиг. 8. Видно,что отраженное р-событие проявляется преимущественно в вертикальной компоненте скорости частицы,в то время как преобразованное s-событие проявляется преимущественно в радиальной компоненте скорости частицы. Происходит небольшая утечка, так что в радиальной компоненте скорости частицы существуют слабые события, соответствующие отраженному р-событию, а в вертикальной компоненте существуют слабые события, соответствующие преобразованному s-событию. Такие события p-s преобразования в радиальной компоненте скорости частицы имеют амплитуду,достаточную для того, чтобы их можно было видеть непосредственно. Путем просмотра радиальной компоненты скорости частицы соответствующее событие р-s преобразования может наблюдаться непосредственно при временах выборки, находящихся вскоре после того (через 0-100 мс), как сильное рсобытие регистрируется в вертикальной компоненте. Для извлечения разности времен прохождения между р-волной и соответствующей преобразованной s-волной это р-событие может быть отслежено по трассам, получаемым различными приемниками. Время его поступления может быть использовано как эталонное нулевое время для определения начала временного окна, в пределах которого ожидается появление соответствующего s-события в радиальной компоненте. Продолжительность временного окна соответствует максимуму ожидаемой сдвиговой статики. После того как соответствующее s-событие будет локализовано в пределах этого временного окна, непосредственно из данных может быть считана разность времен прохождения между р-событием и s-событием. Альтернативно, s-событие опять может быть отслежено по трассам. Такое р-событие не является плоским в трассах, изображенных на фиг. 9(а) и 9(b), поскольку оно не происходит в один и тот же момент времени в каждой трассе. Фиг. 10(а) и 10(b) демонстрируют данные согласно фиг. 9(а) и 9(b) после того, как они подверглись предварительной обработке для того, чтобы сделать р-событие плоским. Описанные выше варианты воплощения настоящего изобретения относятся к преобразованию мод,происходящему тогда, когда распространяющееся вверх волновое поле претерпевает частичное преобразование мод при прохождении через границу раздела между фундаментом 6 и низкоскоростным слоем 5. Однако преобразование мод может происходить и тогда, когда распространяющаяся вниз волна претерпевает отражение на границе раздела между низкоскоростным слоем 5 и фундаментом 6. Это иллюстрируется схематически на фиг. 12(а) и 12(b). На фиг. 12(а) р-волна 7 распространяется вниз через низкоскоростной слой 5 и отражается на границе раздела между низкоскоростным слоем 5 и фундаментом 6. Частичное преобразование мод при отражении происходит таким образом, что, в дополнение к распространяющейся вверх отраженной р-волне 7', существует также распространяющаяся вверх преобразованная s-волна 7". Описанные выше методики определения разности между статическим сдвигом р-волны и статическим сдвигом s-волны могут быть равным образом применены и для системы разведки, изображенной на фиг. 12(а). Способ может быть применен точно таким же образом, как он применен к распространяющемуся вверх проходящему (и частично преобразованному) волновому полю, за исключением того, что при рассмотрении отражений распространяющегося вверх волнового поля источника должны быть использованы данные, регистрируемые на более коротких смещениях. Это происходит потому, что волновое поле, отраженное глубоким отражателем в модели, используемой в системах с прохождением, при данном смещении будет распространяться ближе к вертикали по сравнению с волнами, отраженными на границе раздела между низкоскоростным слоем 5 и фундаментом и поступающими при том же смещении.-9 005475 Фиг. 12(b) в целом соответствует фиг. 12(а), но на ней показано частичное преобразование мод,происходящее при отражении распространяющегося вниз поля s-волны 8, с получением распространяющегося вверх s-волнового поля 8' и распространяющегося вверх преобразованного р-волнового поля 8". На фиг. 13(а) и 13(b) показаны дополнительные варианты воплощения настоящего изобретения. В этих вариантах воплощения преобразование мод вызывается волновым полем, которое распространяется вниз через низкоскоростной слой 5 и претерпевает критическое преломление на границе раздела между низкоскоростным слоем 5 и фундаментом 6. На фиг. 13(а) представлена р-волна 9, распространяющаяся вниз через низкоскоростной слой 5. Угол i обозначает угол падения, то есть угол между направлением распространения волны и нормалью к границе раздела (или, в случае неплоской границы раздела, мгновенной нормалью к границе раздела). Когда р-волна 9 падает на границу раздела между низкоскоростным слоем 5 и фундаментом, обычно ожидается, что р-волна 9 должна частично проходить в фундамент 6, как показано с помощью прерывистой линии 9', и частично отражаться на границе раздела. Угол r между проходящей волной 9' и нормалью к границе раздела связан с углом падения i законом Снеллиуса. Однако скорость сейсмической энергии в фундаменте 6 будет, как правило, большей, чем скорость сейсмической энергии в низкоскоростном слое 5, так что, если угол падения превышает определенный критический угол, то закон Снеллиуса будет предсказывать значение для sin r, которое является большим, чем 1. Когда это происходит,преломленная волна 9" будет распространяться вдоль границы раздела между низкоскоростным слоем 5 и фундаментом 6. Это явление известно как критическое преломление, и оно является подобным явлению полного внутреннего отражения в оптике. Критически преломленная волна может возбудить в низкоскоростном слое распространяющиеся вверх волны, известные как головные волны, и они регистрируются приемниками 4. Такие головные волны обнаруживаются в сейсмических данных, когда смещение источник-приемник является достаточно большим, так что угол падения i на соответствующую границу раздела превосходит критический угол для наступления критического преломления. В случае границы раздела между низкоскоростным слоем 5 и фундаментом 6 критическое смещение является низким, и головные волны регистрируются для большинства положений приемников и для большинства точек взрываний. В качестве примера рассмотрим низкоскоростной слой 5, состоящий из материала со скоростью распространения р-волны 1 и скоростью распространения s-волны 1, лежащий поверх фундамента 6 со скоростью распространения р-волны 2 и скоростью распространения s-волн 2, когда удовлетворяются следующие отношения: 12 12; и 21 На фиг. 13(а) показан случай, когда угол падения i удовлетворяет отношению sin i=1/2. Этот угол падения представляет собой первый критический угол, и он будет упоминаться как с. Когда угол падения распространяющейся вниз р-волны 9 удовлетворяет уравнению i=с, sin r принимает значение 1, а это означает, что r равен 90. Проходящей волны 9', распространяющейся в фундамент 6, по этой причине не существует, и вместо нее будет существовать критически преломленная р-волна 9", которая распространяется вдоль границы раздела между низкоскоростным слоем 5 и фундаментом 6. Поскольку эта преломленная р-волна распространяется вдоль границы раздела, она будет возбуждать р-волны 10 в низкоскоростном слое 5, и эти р-волны 10 будут давать начало s-волнам 11 в результате преобразования мод на границе раздела между низкоскоростным слоем 5 и фундаментом 6. На фиг. 13(а) изображены две пары из р-волны 10 и преобразованной s-волны 11, хотя возбуждение р- иs-волн 10, 11 будет происходить во всех точках вдоль пути преломленной волны 9". Возбужденная рволна 10 будет распространяться под углом с к нормали к границе раздела, в то время как преобразованная s-волна 11 будет распространяться под углом с к нормали, где sin(с)=1/2. Критически преломленная р-волна существует для всех углов падения, равных или больших, чем с. Фиг. 13(b) в целом соответствует фиг. 13(а), но на ней показано критическое преломление s-волны 12, распространяющейся вниз через низкоскоростной слой 5. Если угол падения i является равным или большим, чем второй критический угол с, то распространяющаяся вниз s-волна 12 будет частично преобразовываться в критически преломленную р-волну 12", которая распространяется вдоль границы раздела между низкоскоростным слоем 5 и фундаментом 6. Критически преломленная р-волна снова возбуждает р-волны 13 и преобразованные s-волны 14, которые распространяются вверх через низкоскоростной слой 5. В случаях, показанных на фиг. 13(а) и 13(b), критически преломленная р-волна 9", 12" излучает как энергию р-волны, так и энергию s-волны вверх в низкоскоростной слой 5. Детектирование энергии как рволны, так и s-волны приемниками 4 опять дает разницу между статическим сдвигом р-волн и статическим сдвигом s-волн. Такая s-волна 11 формируется путем преобразования мод на границе раздела между низкоскоростным слоем 5 и фундаментом 6, так что разность времен прохождения через низкоскоро- 10005475 стной слой между р-волной 10, возбужденной в одной из точек на границе раздела, и соответствующей преобразованной s-волной 11 будет равна разности между статическим сдвигом р-моды и статическим сдвигом s-моды. Таким образом, разность во времени между временем поступления р-волны 10 (рпоступление) и временем поступления s-волны 11 (s-поступление) в приемники 4 равна разности между статическим сдвигом р-волны и статическим сдвигом s-волны. Таким образом, путем детектирования времени задержки между р-поступлением и соответствующим s-поступлением возможно определить разность между р- и s-статическими сдвигами. Время задержки может быть определено с использованием любого из описанных выше способов, например путем деконволюции или взаимной корреляции вертикальной и радиальной компонент измеренной скорости частицы или путем деконволюции или взаимной корреляции измеренного давления и радиальной компоненты скорости частицы. Необходимо заметить, что существует и третий критический угол, т.е. с, где sin с=1/2. Когда распространяющаяся вниз s-волна падает на границу раздела между низкоскоростным слоем 5 и фундаментом 6 под углом, равным или большим, чем этот третий критический угол с, будет генерироваться критически преломленная s-волна, которая распространяется вдоль границы раздела между низкоскоростным слоем 5 и фундаментом 6. Однако эта критически преломленная s-волна не будет вызывать распространяющуюся вверх р-волну в низкоскоростном слое 5. В примерах на фиг. 13(а) и 13(b) предполагается, что как скорость р-волны, так и скорость s-волны увеличиваются при движении вниз через границу раздела между низкоскоростным слоем и фундаментом. Однако иногда обнаруживается, что не все скорости увеличиваются при движении вниз через границу раздела, и это приводит к ситуациям, которые являются более сложными, чем те, которые изображены на фиг. 13(а) и 13(b). Однако способ по настоящему изобретению можно применять всякий раз,когда на границе раздела между низкоскоростным слоем 5 и фундаментом 6 происходит критическое преломление, приводящее к возбуждению в низкоскоростном слое 5 как р-волн, так и преобразованныхs-волн. Сейсмическая энергия, собираемая в приемниках 4, может анализироваться с помощью любого описанного выше способа для извлечения разности времен прохождения испускаемой энергии р-волны и энергии s-волны и, таким образом, для получения разности между статическим сдвигом р-волны и статическим сдвигом s-волны. В приведенных выше примерах настоящего изобретения разность времен прохождения р-мод и sмод через низкоскоростной слой получают из двух компонент скорости частицы или из одной компоненты скорости частицы и давления. В альтернативных вариантах воплощения настоящего изобретения вместо скорости частицы используют измеренное смещение частицы. Таким образом, в альтернативном варианте воплощения разность времен прохождения р-мод и s-мод через низкоскоростной слой получают из двух компонент смещения частицы, например путем взаимной корреляции или деконволюции вертикальной и радиальной компонент смещения частицы. В дополнительном варианте воплощения разность времен прохождения р-мод и s-мод через низкоскоростной слой получают из одной компоненты смещения частицы и давления, например путем взаимной корреляции или деконволюции компоненты смещения частицы и давления. Эти варианты воплощения, в целом, соответствуют тем вариантам воплощениям, которые описаны выше в отношении использования скорости частицы для получения разности времен прохождения, и поэтому далее они описываться не будут. В дополнительных вариантах воплощения настоящего изобретения вместо скорости частицы используют ускорение частицы. Таким образом, в этих вариантах воплощения разность времен прохождения р-мод и s-мод через низкоскоростной слой получают из двух компонент ускорения частицы или из одной компоненты ускорения частицы и давления. Это может быть сделано, например, путем взаимной корреляции или деконволюции вертикальной и радиальной компонент ускорения частицы или путем взаимной корреляции или деконволюции одной из компонент ускорения частицы и давления. Эти варианты воплощения, в целом, соответствуют тем вариантам воплощения, которые описаны выше в отношении использования скорости частицы для получения разности времен прохождения, и поэтому далее они описываться не будут. В принципе, могли бы использоваться и высшие производные смещения частицы, так что разность времен прохождения р-мод и s-мод через низкоскоростной слой могла бы быть получена из двух компонент высшей производной смещения частицы или из одной из компонент высшей производной смещения частицы и давления. В описанных выше вариантах воплощения приемник(и) 4 располагается(ются) на донном грунте(т.е. морском дне). В принципе, однако, настоящее изобретение не ограничивается этим, и приемник(и) может/могут располагаться где угодно на пути сейсмической энергии после того, как произошло частичное преобразование мод. Таким образом, в описанных выше вариантах воплощения приемник(и) может/могут, в принципе, располагаться в любом месте над границей раздела между низкоскоростным слоем и фундаментом. Например, приемник(и) может/могут быть зарыт(ы) в донном грунте при 4-мерной цейтраферной сейсмической разведке. Необходимо заметить, однако, что для определения разности между статическими сдвигами двух мод потребуется дополнительная обработка в том случае, если приемник не расположен на донном грунте.- 11005475 В приведенном выше описании вариантов воплощения настоящего изобретения х-направление было определено как радиальное направление, а именно проекция направления от источника к приемнику на донный грунт. Сейсмический приемник будет регистрировать компоненты скорости частицы или смещения частицы в двух ортогональных горизонтальных направлениях, и эти направления могут быть обозначены как х- и у-оси приемника. Необходимо заметить, что приемник может и не быть установлен так, что его х-ось совмещена с радиальным направлением, так что для вычисления из выходного сигнала приемника радиальных и поперечных компонент может потребоваться дополнительное проецирование или вращение измеренных приемником исходных горизонтальных компонент. Измеренные приемником исходные х- и z-компоненты могут, в принципе, использоваться для оценки разности времен прохождения р-мод и s-мод через низкоскоростной слой с помощью любого из описанных выше способов, и это будет давать приемлемые результаты в том случае, если угол между х-осью приемника и радиальным направлением не является слишком большим. Способ обработки согласно настоящему изобретению может быть применен для обработки существовавших ранее сейсмических данных. Настоящее изобретение может также быть включено в способ сейсмической разведки, в котором собранные сейсмические данные включают в себя события получения преобразованных мод, возникающих из-за частичного преобразования мод на границе раздела между низкоскоростным слоем и фундаментом. Такие сейсмические данные могут быть получены путем направления сейсмической энергии к границе раздела таким образом, что при прохождении сейсмической энергии через границу раздела или при отражении от нее происходит частичное преобразование мод. Как указано выше, настоящее изобретение может быть применено для распространяющихся вверхs-волн, которые претерпевают частичное преобразование мод в р-волны на границе раздела между фундаментом и низкоскоростным слоем, а также к распространяющимся вверх р-волнам, которые претерпевают частичное преобразование мод в s-волны на указанной границе раздела. Оба этих случая могут быть использованы для нахождения разности времен прохождения s-волн и р-волн через низкоскоростной слой. В принципе, разность времен прохождения через низкоскоростной слой между р-волнами и sволнами будет зависеть от замедленности падающей волны даже для идеальной ситуации, в которой донный грунт и граница раздела между низкоскоростным слоем и фундаментом являются плоскими. На практике, однако, для всего диапазона смещений между источником и приемником, которые используются при типичных сейсмических изысканиях, и для типичных материалов низкоскоростного слоя зависимость разности времен прохождения от замедленности падающей волны является, как обнаружено,очень слабой и разность времен прохождения является фактически постоянной, безотносительно к замедленности падающей волны. Обнаружено, что для полного диапазона приемников типичный разброс составляет всего лишь несколько миллисекунд. По этой причине разумно предположить, что разность времен прохождения через низкоскоростной слой не зависит от замедленности падающей волны. Это предположение является эквивалентным поверхностно-согласованному статическому приближению, которое, как правило, рассматривается как правильное в большинстве статических задач. Это предположение означает, что волновое поле, регистрируемое на приемниках 4, не должно фильтроваться или выделяться, даже если оно представляет собой смесь волн, поступающих с различными замедленностями от различных отражателей внутри земли, и также представляет собой смесь как р-волн, так и s-волн. В дополнение к этому, для применения способа по настоящему изобретению не является необходимым извлекать отдельные события из данных. Это означает, что в процессе деконволюции или взаимной корреляции могут быть использованы относительно длительные временные окна (порядка секунд). Более того, трассы могут суммироваться в общей области приемника для улучшения отношения сигнала к шуму у этих данных. Выше отмечено, что приемник может и не использоваться таким образом, что его х-ось совмещена с радиальным направлением, и в этом случае измеряемые приемником исходные горизонтальные компоненты могут потребовать проецирования или вращения для вычисления из выходного сигнала приемника радиальных и поперечных компонент. Необходимо заметить, что эта процедура не является необходимой там, где сейсмические силы расположены на линии взрываний и приводятся в действие, когда линия взрыва расположена параллельно линии приемников и непосредственно над ней. В этом случае хкомпонента исходных данных приемника может рассматриваться как радиальная компонента, а укомпонента исходных данных приемника может рассматриваться как поперечная компонента. Необходимо заметить, однако, что при изменении смещения с отрицательного на положительное (приближение одномерной земли) изменяется знак данных для событий на радиальной компоненте сейсмических данных. Это необходимо принять во внимание при суммировании результатов деконволюции или взаимной корреляции от общего ряда приемников, который включает в себя как положительные, так и отрицательные смещения, поскольку пренебрежение этим может приводить к ослаблению и, возможно, даже к полной потере сигнала. Перед суммированием данных необходимо изменить на противоположный знак данных, собранных для одного смещения. Настоящее изобретение было описано выше со ссылками на взаимную корреляцию, деконволюцию и бикогерентность в качестве примеров способов установления корреляции двух трасс сейсмических- 12005475 данных. Однако настоящее изобретение не ограничено этими способами, и поэтому может быть использована любая пригодная методика или алгоритм для установления корреляции между трассами. Настоящее изобретение может также быть применено к событиям реверберации р- и s-волн. При таких событиях частичное преобразование мод происходит, например, при прохождении через или при отражении от границы раздела между ближней поверхностью 5 и фундаментом 6. При событии реверберации, однако, р- и s-волны не проходят к приемнику непосредственно от точки, в которой происходит преобразование мод. При событии реверберации р- и s-волны претерпевают одно или несколько отражений на границе подповерхностного слоя и, таким образом, осуществляют многократные прохождения через подповерхностный слой 5 перед попаданием на приемник. Сейсмические данные также могут содержать события, которые возникают из-за частичного преобразования мод, которое происходит не на границе подповерхностного слоя, а внутри этого подповерхностного слоя. Это может происходить, например, благодаря эффектам, вызываемым многослойностью в подповерхностном слое. В описанных выше вариантах воплощения одно из событий, обусловленных р- и s-волнами, генерируется за счет частичного преобразования мод. Однако, как отмечено выше, настоящее изобретение не ограничивается этим и может быть применено к любой паре соответствующих р- и s-событий, для которых различия в амплитуде и/или в форме сигнала двух событий возникают преимущественно из-за подповерхностного слоя 5. Например, при сейсмической разведке, которая использует источник, который испускает как р-, так и s-волны, собранные данные будут содержать событие, возникающее из-за отражения р-волн в конкретной точке на нижней границе подповерхностного слоя 5, а также будут содержать соответствующее событие, возникающее из-за отражения s-волн в этой точке на нижней границе подповерхностного слоя 5. Настоящее изобретение может быть применено к этим событиям. Необходимо отметить, что, поскольку столб воды не будет поддерживать распространение s-волны, испускающий какs-, так и р-моды источник должен быть размещен на донном грунте, т.е. на поверхности земли или в скважине. Фиг. 14 представляет собой блок-схему программируемого устройства 15 в соответствии с настоящим изобретением. Устройство содержит программируемый процессор 16 обработки данных с программируемой памятью 17, например, в виде постоянного запоминающего устройства ПЗУ, хранящего программу для управления процессором 16 обработки данных для осуществления любого из описанных выше способов обработки. Устройство дополнительно содержит нестираемую оперативную память 18 для хранения, например, любых данных, которые могут быть сохранены в отсутствие подачи питания. Рабочая или сверхоперативная память процессора обработки данных обеспечивается с помощью памяти 19 с произвольной выборкой (ОЗУ). Предусмотрен входной интерфейс 20, например, для получения команд и данных. Предусмотрен выходной интерфейс 21, например, для отображения информации, относящейся к ходу исполнения и к результатам способа. Сейсмические данные для обработки могут вводиться через входной интерфейс 20 или могут, альтернативно, извлекаться из машинно-считываемого носителя 22 данных. Программа для работы системы и для осуществления описанного выше способа хранится в памяти 17 программ, которая может быть воплощена в виде полупроводниковой памяти, например, хорошо известного типа ПЗУ. Однако программа может храниться и на любом другом пригодном для использования носителе информации, таком как магнитный носитель 17 а данных (дискета) или компакт-диск 17b. ФОРМУЛА ИЗОБРЕТЕНИЯ 1. Способ обработки сейсмических данных, содержащих первую и вторую моды сейсмической энергии, причем первая и вторая моды сейсмической энергии распространялись через слой внутреннего пространства земли вдоль, по существу, одного и того же пути, включающий в себя стадию обработки сейсмических данных для получения разности времен прохождения через слой внутреннего пространства земли между сейсмической энергией, распространяющейся в первой моде, и сейсмической энергией,распространяющейся во второй моде. 2. Способ по п.1, в котором вторая мода была сгенерирована за счет частичного преобразования первой моды на границе указанного слоя. 3. Способ по п.2, в котором указанная граница представляет собой нижнюю границу слоя. 4. Способ по п.2 или 3, в котором указанный слой представляет собой поверхностный слой или подповерхностный слой. 5. Способ по любому из пп.1-4, в котором стадия обработки сейсмических данных включает в себя определение первого и второго параметров, являющихся показателями первого и второго аспектов сейсмических данных; и обработку первого и второго параметров для получения разности времен прохождения через указанный слой между сейсмической энергией, распространяющейся в первой моде, и сейсмической энергией, распространяющейся во второй моде. 6. Способ по п.5, в котором стадия обработки включает в себя установление взаимной корреляции первого и второго параметров друг с другом для получения тем самым разности времен прохождения- 13005475 через указанный слой между сейсмической энергией, распространяющейся в первой моде, и сейсмической энергией, распространяющейся во второй моде. 7. Способ по п.5, в котором стадия обработки включает в себя деконволюцию первого и второго параметров для получения тем самым разности времен прохождения через указанный слой между сейсмической энергией, распространяющейся в первой моде, и сейсмической энергией, распространяющейся во второй моде. 8. Способ по п.5, в котором стадия обработки включает в себя обработку первого и второго параметров с использованием метода оценки времени задержки бикогерентности для получения тем самым разности времен прохождения через указанный слой между сейсмической энергией, распространяющейся в первой моде, и сейсмической энергией, распространяющейся во второй моде. 9. Способ по п.5, в котором стадия обработки включает в себя определение времени прохождения через указанный слой первой моды из первого параметра и определение времени прохождения через указанный слой второй моды из второго параметра. 10. Способ по любому из пп.5-9, в котором первый параметр представляет собой первую компоненту скорости частицы, а второй параметр представляет собой вторую компоненту скорости частицы, при этом первая и вторая компоненты скорости частицы не параллельны друг другу. 11. Способ по п.10, в котором первая компонента скорости частицы представляет собой, по существу, вертикальную компоненту, а вторая компонента скорости частицы представляет собой, по существу,горизонтальную компоненту. 12. Способ по любому из пп.5-8, в котором первый параметр представляет собой первую компоненту скорости частицы, а второй параметр представляет собой давление. 13. Способ по любому из пп.5-9, в котором первый параметр представляет собой первую компоненту смещения частицы, а второй параметр представляет собой вторую компоненту смещения частицы,при этом первая и вторая компоненты смещения частицы не параллельны друг другу. 14. Способ по п.13, в котором первая компонента смещения частицы представляет собой, по существу, вертикальную компоненту, а вторая компонента смещения частицы представляет собой, по существу, горизонтальную компоненту. 15. Способ по любому из пп.2-14, в котором первая мода представляет собой моду р-волны, а вторая мода, сгенерированная за счет частичного преобразования мод, представляет собой моду s-волны. 16. Способ по любому из пп.2-14, в котором первая мода представляет собой моду s-волны, а вторая мода, сгенерированная за счет частичного преобразования мод, представляет собой моду р-волны. 17. Способ по п.5, в котором стадия определения первого и второго параметров включает в себя разложение полученных сейсмических данных на р-компоненту и s-компоненту. 18. Способ по п.17, в котором стадия обработки включает в себя установление взаимной корреляции р-компоненты и s-компоненты сейсмической энергии друг с другом для получения тем самым разности времен прохождения через указанный слой между сейсмической энергией, распространяющейся в р-моде, и сейсмической энергией, распространяющейся в s-моде. 19. Способ сейсмической разведки, включающий в себя стадии направления сейсмической энергии, распространяющейся в первой моде, к границе слоя донного грунта таким образом, что на указанной границе происходит частичное преобразование моды сейсмической энергии во вторую моду; сбора сейсмических данных, содержащих первую и вторую моды сейсмической энергии, на одном или нескольких приемниках; и обработки сейсмических данных в соответствии со способом по любому из пп.1-18. 20. Способ сейсмической разведки по п.19, в котором стадия сбора сейсмических данных включает в себя сбор сейсмических данных на одном или нескольких приемниках, размещенных над указанным слоем. 21. Способ по п.19 или 20, в котором стадия сбора сейсмических данных включает в себя сбор сейсмических данных на одном или нескольких приемниках, расположенных на донном грунте. 22. Способ сейсмической разведки по любому из пп.19-21, включающий в себя направление сейсмической энергии вниз через указанный слой, в результате чего частичное преобразование мод происходит при отражении на границе сейсмической энергии первой моды, распространяющейся вниз через указанный слой. 23. Способ сейсмической разведки по любому из пп.19-21, включающий в себя направление сейсмической энергии вверх через указанный слой, в результате чего частичное преобразование мод происходит при прохождении через границу распространяющейся вверх сейсмической энергии первой моды. 24. Способ сейсмической разведки по любому из пп.19-21, включающий в себя направление сейсмической энергии к границе указанного слоя под таким углом, чтобы вызвать распространение сейсмической энергии в первой моде вдоль указанной границы. 25. Способ сейсмической разведки по п.24, в котором первая мода представляет собой моду рволны и угол i между направлением распространения падающей сейсмической энергии и нормалью к- 14005475 границе удовлетворяет условию sin i1/2, где 1 представляет собой скорость распространения р-волн в слое донного грунта, а 2 представляет собой скорость распространения р-волн в нижележащем слое. 26. Способ сейсмической разведки по п.24, в котором первая мода представляет собой моду рволны и угол i между направлением распространения падающей сейсмической энергии и нормалью к границе удовлетворяет условию sin i1/2, где 1 представляет собой скорость распространения s-волн в слое донного грунта, a 2 представляет собой скорость распространения р-волн в нижележащем слое. 27. Устройство для обработки сейсмических данных, включающих в себя первую и вторую моды сейсмической энергии, причем вторая мода сгенерирована за счет частичного преобразования первой моды на границе слоя донного грунта, в результате чего первая и вторая моды сейсмической энергии распространялись через слой донного грунта вдоль, по существу, одного и того же пути, содержащее средство для обработки сейсмических данных с получением разности времен прохождения через указанный слой между сейсмической энергией, распространяющейся в первой моде, и сейсмической энергией,распространяющейся во второй моде. 28. Устройство по п.27, содержащее средство для определения первого и второго параметров, являющихся показателями первого и второго аспектов сейсмических данных; и средство для обработки первого и второго параметров с получением разности времен прохождения через указанный слой между сейсмической энергией, распространяющейся в первой моде, и сейсмической энергией, распространяющейся во второй моде. 29. Устройство по п.28, в котором средство обработки приспособлено для установления взаимной корреляции первого и второго параметров друг с другом для получения тем самым разности времен прохождения через указанный слой между сейсмической энергией, распространяющейся в первой моде, и сейсмической энергией, распространяющейся во второй моде. 30. Устройство по п.28, в котором средство обработки приспособлено для деконволюции первого и второго параметров для получения тем самым разности времен прохождения через указанный слой между сейсмической энергией, распространяющейся в первой моде, и сейсмической энергией, распространяющейся во второй моде. 31. Устройство по п.28, в котором средство обработки приспособлено для обработки первого и второго параметров с использованием способа оценки времени задержки бикогерентности для получения разности времен прохождения через указанный слой между сейсмической энергией, распространяющейся в первой моде, и сейсмической энергией, распространяющейся во второй моде. 32. Устройство по п.28, в котором средство обработки приспособлено для определения времени прохождения через указанный слой первой моды из первого параметра и для определения времени прохождения через указанный слой второй моды из второго параметра. 33. Устройство по любому из пп.27-32, содержащее программируемый процессор обработки данных.