Бесплатный автореферат и диссертация по наукам о земле на тему

Алгоритмы и технология обработки совмещенных наземно-скважинных сейсмических наблюдений

ВАК РФ 25.00.10, Геофизика, геофизические методы поисков полезных ископаемых

Содержание диссертации, кандидата технических наук, Баранов, Константин Владимирович

Введение

Глава 1. Состояние вопроса и постановка основных задач

1.1 Приемы коррекции статических поправок в сейсмограммах 3D

1.2 Описание методики «Локальный проект 3D+BC11»

1.3 Построение скоростной модели среды по годографам ВСП

Глава 2. Методика оценки статических поправок по годографу первых вступлений прямой волны

2.1 Выбор модели поля времен первых вступлений прямой волны

2.2 Описание методики

2.3 Модельный эксперимент

Глава 3. Восстановление скоростной модели среды по годографам первых вступлений ВСП

3.1 Описание методики

3.2 Модельный эксперимент

Глава 4. Методика коррекции формы импульса возбуждения

4.1 Описание методики

4.2 Модельный эксперимент

Глава 5. Технология обработки совмещенных наземно-скважинных наблюдений

5.1 Кинематическая обработка данных ВСП

5.2 Обработка записей глубинного контрольного зонда ВСП

5.3 Обработка реальных данных 3D на поверхности с применением результатов обработки записей глубинного прибора

Введение Диссертация по наукам о земле, на тему "Алгоритмы и технология обработки совмещенных наземно-скважинных сейсмических наблюдений"

Основными геофизическими методами, используемыми при поисках, разведке и мониторинге нефтяных и газовых месторождений, являются геофизические исследования скважин (ГИС) и сейсморазведка на поверхности (2D и 3D).

Методами ГИС изучают как весь разрез, вскрытый скважиной, так и в особенности продуктивные интервалы. Разрешающая способность по глубине таких исследований составляет 10-20 см, и, в большинстве случаев, с необходимой точностью и детальностью оцениваются все основные параметры геологической среды, необходимые при подсчете запасов и составлении проекта разработки месторождения.

При этом количество глубоких скважин, как правило, недостаточно для оценки поведения продуктивных пластов в межскважинном пространстве. Для построения непрерывной модели строения продуктивного пласта обычно используется увязка результатов ГИС с результатами сейсморазведки на поверхности, которая позволяет коррелировать отражения на непрерывных 3D или 2D изображениях объекта. Однако разрешенность сейсморазведки на поверхности составляет в лучшем случае около 20 метров, что на два порядка хуже, чем ГИС.

Вертикальное сейсмическое профилирование (ВСП), будучи одновременно сейсмическим и скважинным методом, находится на стыке наземной сейсморазведки и ГИС. ВСП позволяет, с одной стороны, достичь при изучении разреза вдоль ствола скважины высокой разрешенности и, с другой - распространить результаты на окрестности скважины. Метод ВСП предложен Е. И. Гальпериным и развит им с участием большой группы исследователей (Г. А. Шехтман, Л. Л. Худзинский, А. А. Табаков, В. А. Теплицкий, Г. Е. Руденко и др.).

Основными направлениями применения ВСП в настоящее время являются:

- изучение упруго-плотностных характеристик среды на сейсмических частотах;

- динамическая привязка отражений, регистрируемых на поверхности, к литологическому разрезу;

- детальное изучение околоскважинного пространства с использованием более широкого спектра частот, чем при сейсморазведке на поверхности;

- прогнозирование геологического разреза ниже забоя скважины [54].

Область околоскважинного пространства, которая может быть изучена ВСП, принципиально ограничена расстоянием порядка половины глубины залегания изучаемого горизонта из-за трудностей решения прямых и обратных задач при больших углах падения сейсмического луча на границу.

Таким образом, сейсморазведка на поверхности обладает возможностью получения связного изображения больших объектов, но недостаточной разрешенностью и точностью для решения задач высокой детальности. ВСП обеспечивает большую детальность изучения околоскважинного пространства, но в ограниченном диапазоне удалений от скважины.

При эксплуатационном бурении на нефть и газ часто возникает потребность детального изучения продуктивных пластов на расстояниях до 1.5 - 2 км от скважины, не достигаемых для исследований методом ВСП. Наличие глубокой скважины может существенно повысить точность и детальность исследований 2D- или ЗБ-сейморазведки в окрестностях этой скважины, если использовать одновременные наблюдения на поверхности и в скважине, то есть совместить преимущества сейсморазведки на поверхности и ВСП. Методика совмещенных наблюдений, при которой возбуждения на поверхности на сетке работ 3D регистрируются зондом ВСП, расположенным вблизи забоя глубокой скважины, названа "Локальный проект 3D+BC11". Также методика предполагает проведение ВСП из нескольких пунктов возбуждения.

Ожидаемые преимущества таких наблюдений можно сформулировать следующим образом:

- возможность компенсации влияния неоднородностей условий возбуждения за счет контроля формы импульса каждого воздействия и деконволюции по форме сигнала;

- использование достоверных статических поправок за неоднородности условий возбуждения;

- использование дополнительной информации о скоростях, полученной путем построения куба скоростей с точным скоростным законом на скважине по данным ВСП и с учетом горизонтального градиента скоростей, оцениваемого по вариации времен прихода прямой волны при возбуждении на сетке работ 3D.

Несмотря на то, что известны примеры совмещенных наблюдений (MDSeis, 1994), в целом отсутствуют опубликованные методики и технологии обработки и интерпретации таких наблюдений, а в используемых пакетах программ обработки данных наземной сейсморазведки не представлены специализированные блоки, обеспечивающие оптимальное использование преимуществ совместных наблюдений.

Поэтому, целью данной работы является разработка алгоритмов, программ, методики и технологии обработки совмещенных наземно-скважинных наблюдений для повышения детальности и точности изучения околоскважинного пространства.

Основные задачи исследования: 1. Разработка алгоритмов и программ определения статических поправок за различия условий возбуждения и оценки локальных аномалий скоростей верхней части разреза.

2. Разработка алгоритмов и программ оптимизационного решения обратных кинематических задач ВСП на прямых и головных волнах в горизонтально-слоистых трансверсально-изотропных средах одновременно по годографам от нескольких пунктов возбуждения.

3. Разработка алгоритмов и программ оценки эталонных импульсов и компенсации влияния неоднородностей условий возбуждения в записях наземных приемников.

4. Оценка эффективности разработанных алгоритмов и программ на модельных данных.

5. Опробование комплекса программ специализированной обработки наземно-скважинных наблюдений на реальных данных.

Защищаемые положения:

1. Технология кинематической обработки наземно-скважинных наблюдений «Локальный проект 3D+BCII» повышает достоверность результатов сейсморазведки 3D в окрестности глубокой скважины путем более точной оценки и учета статических поправок, связанных с локальными аномалиями скоростей ВЧР.

2. Решение обратной кинематической задачи на прямых и головных волнах в горизонтально-слоистых трансверсально-изотропных средах на основе одновременной оптимизации модели среды по годографам первых вступлений от нескольких пунктов возбуждения повышает точность и достоверность определения параметров модели.

3. Методика компенсации влияния неоднородностей условий возбуждения, опирающаяся на регистрацию колебаний зондом ВСП, позволяет корректировать изменения формы импульса возбуждения в записях наземных сейсмоприемников в окрестности глубокой скважины.

Работа состоит из пяти глав, первая из которых содержит описание проблем в обработке данных наземной сейсморазведки, связанных с неоднородностью ВЧР и неоднородностями условий возбуждения. Описана методика «Локальный проект ЗБ+ВСП», позволяющая оценить статические поправки за ВЧР без ограничения длины периода, компенсировать неоднородности условий возбуждения и получить скоростную модель среды с учетом вариаций скоростей ВЧР.

Во второй главе обосновывается выбор модели поля времен первых вступлений прямой волны, регистрируемой расположенным на забое глубокой скважины зондом при возбуждении на сетке работ 3D на поверхности. Описываются алгоритмы определения высокочастотных статических поправок за возбуждения и оценки локальных аномалий скоростей ВЧР.

В третьей главе поставлена и решена задача восстановления пластовой модели среды по годографам первых вступлений прямых и головных волн из нескольких пунктов взрыва. Решение задачи ищется в классе анизотропных параллельно-слоистых сред и сводится к нахождению минимума функции многих переменных.

Четвертая глава посвящена разработке алгоритмов оценки эталонных импульсов для компенсации влияния неоднородностей условий возбуждения в записях наземных приемников.

В пятой главе описан граф обработки совмещенных наземно-скважинных наблюдений и приведены результаты опробования на реальном объекте.

Автор выражает благодарность научному руководителю, доктору физико-математических наук, профессору МГГРУ Б. Е. Лухминскому за внимание, помощь и поддержку, оказанную автору в процессе подготовки данной работы. Автор выражает признательность и благодарность начальнику отдела ВСП ОАО «ЦГЭ» А. А. Табакову как автору многих идей, положенных в основу данной разработки, которые впоследствии были развиты автором и составляют неотъемлемую часть данной работы.

Заключение Диссертация по теме "Геофизика, геофизические методы поисков полезных ископаемых", Баранов, Константин Владимирович

Основные выводы по пятой главе:

- разработана технология обработки данных совмещенных наземно-скважинных сейсмических наблюдений по методике «Локальный проект 3D+BCn»;

- применение алгоритмов и технологии оптимизационного решения обратных кинематических задач для годографов первых вступлений из нескольких пунктов возбуждения решает задачу оценки правдоподобной модели среды при отсутствии измерений времен для части наблюдений из ближнего пункта возбуждения;

- применение алгоритмов и технологии кинематической обработки времен первых вступлений глубинного прибора к данным реального полевого эксперимента позволило выделить правдоподобные среднечастотные аномалии времен, хорошо коррелирующие с руслами речной системы;

- применение алгоритмов и технологии компенсации влияния неоднородностей условий возбуждения к данным реального полевого эксперимента позволило определить набор индивидуальных для каждого возбуждения устойчивых операторов, обеспечивающих выравнивание формы сигнала в записях наземных приборов без заметного усиления шумов.

Заключение

Выполнен комплекс исследований на модельных и реальных материалах совместных наземно-скважинных наблюдений, направленных на повышение точности, детальности и достоверности изучения околоскважинного пространства, и получены следующие основные результаты:

1) Предложены и реализованы алгоритмы и технология, позволяющие использовать преимущества комплексирования наблюдений сейсморазведки 3D и ВСП.

2) Для. повышения точности оценки скоростной модели среды усовершенствованы алгоритмы оптимизационного решения обратной кинематической задачи. В результате совместной оптимизации параметров модели по годографам нескольких пунктов возбуждения получены решения, существенно уменьшающие невязки расчетных и фактических времен первых вступлений для совокупности наблюдений ВСП из нескольких пунктов возбуждения.

3) Разработаны алгоритмы и технология оценки высокочастотных и среднечастотных статических поправок по записям времен первых вступлений глубинного контрольного прибора. На реальных данных выделены аномалии времен первых вступлений, хорошо коррелирующихся с речной сетью. Амплитуда аномалий ± 8 мс и их знак соответствует задержке в зонах растепления, приуроченных к поймам рек.

4) Разработаны алгоритмы и технология расчета операторов выравнивания формы импульса возбуждения по записям глубинного контрольного прибора. В результате достигнута стандартизация условий возбуждения, которая выразилась в существенном повышении коррелируемости осей синфазности в первичных сейсмограммах.

5) Разработана технология применения результатов обработки данных глубинного контрольного прибора к сейсмограммам, полученным на поверхности. В результате применения разработанной технологии к реальным данным получено заметное повышение разрешенности, качества отражений и закономерное уточнение структурного плана в кубе 3D+BCTI съемки.

6) На основании теоретических и экспериментальных исследований можно рекомендовать применение методики «Локальный проект ЗЭ+ВСП» в задачах детального изучения околоскважинного пространства на расстояниях до 2 км от скважины.

7) Комплекс программ специализированной обработки совмещенных наземно-скважинных наблюдений «Локальный проект 3D+BCn» включен в пакет обработки и интерпретации геофизических данных «ЮНИВЕРС».

8) Дальнейшие исследования в области обработки данных совмещенных наземно-скважинных наблюдений могут быть связаны с изучением поляризационных характеристик прямых и отраженных волн, зарегистрированных глубинным зондом ВСП в скважине.

Библиография Диссертация по наукам о земле, кандидата технических наук, Баранов, Константин Владимирович, Москва

1. Аладжанц Г.С., Марков В.Ф. Комплексирование вертикального сейсмо-профилирования с наземными наблюдениями MOB// Нефтегазовая геология и геофизика. М.: ВНИИОЭНГ, 1979. - Вып. 5. - С. 25-26.

2. Баранов К.В., Душутин А.К. Обратная кинематическая задача при обработке данных ВСП, Тезисы докладов IV международной конференции "Новые идеи в науках о Земле", Москва, 1999, С. 252.

3. Баранов К.В., Душутин А.К. Оценка разрешающей способности ВСП с использованием математического моделирования, Тезисы докладов научно-практической конференции «Состояние и перспективы развития метода ВСП», Москва, 2001, С. 67.

4. Баранов К.В., Душутин А.К., Ференци В.Н. Селекция волновых полей ВСП с использованием кинематического моделирования, Тезисы докладов международной конференции молодых ученых, специалистов и студентов «ГЕОФИЗИКА 2001», Новосибирск, 2001, С. 96.

5. Баранов К.В., Табаков А.А., Душутин А.К. Восстановление детальных скоростных характеристик среды методом ВСП, "Геоинформатика" № 4, Москва, 2001, С. 57-62.

6. Белоусов О. Н., Братчик Р. Ф. Осреднение вертикального годографа непрерывной ломаной по способу наименьших квадратов//Вопросы разведочной и промысловой геофизики. Саратов, 1971. - С. 88-92.

7. Берзон И. С., Епинатьева А. М., Карус Е. В. и др. Методика комплексных скважинных и наземных сейсмических исследований// Состояние и задачи разведочной геофизик. М.: 1970. - С. 459-464.

8. Берзон И. С., Епинатьева А. М., Парийская Г. Н., Стародубровская С. П. Динамические характеристики сейсмических волн в реальных средах. М., Изд. академии наук СССР, 1962.

9. Ю.Бляс Э. А. Обратная кинематическая задача для слоистых трехмерных сред. Геофиз. журн. 1992. Т. 14, № 5. С. 73-77

10. П.Бляс Э. А. Определение скоростной характеристики среды по данным МОВ-ВСШГеология и геофизика. 1987. - № з. С. 88-98.

11. Бляс Э. А., Луговкин С. Б. Прямые и обратные кинематические задачи ВСП в трехмерных сложно построенных средах // Методы расчета и интерпретации сейсмических волновых полей. Новосибирск: Наука, 1991. С. 188-211.

12. Бляс Э.А., Рожков А.С. Построение скоростных моделей по данным вертикального сейсмического профилирования, Вестник МГТУ, 2000, том 3, № 1, стр. 81-86.

13. М.Боголюбский А. Д. Методика выбора наилучшего разбиения вертикального годографа при обработке данных сейсмокаротажа // Разведочная геофизика. -М.: 1980. Вып. 88. - С. 86 - 97.

14. Вентцель Е.С. Теория вероятностей 7-е изд. стер., М., Высш. шк., 2001.

15. Гальперин Е. И. Вертикальное сейсмическое профилирование, М., Недра, 1971.

16. Гальперин Е. И., Вертикальное сейсмическое профилирование. Опыт и результаты. М., Наука, 1994.

17. Гальперин Е. И. К вопросу о влиянии дневной поверхности и верхней части разреза на характер и структуру сейсмограмм. Сб. «Вопросы динамической теории распространения сейсмических волн», сб. VII. Л., изд-во ЛГУ, 1964.

18. Гольдин С. В., 1979, Интерпретация данных сейсмического метода отраженных волн: М., Недра.

19. Гриншпун А. В. и др., 2000, Изучение геологического строения неокомских и юрских отложений на Сугмутском месторождении по результатам сейсмической съемки MOTT-3D сезона 1998-1999 и данных ГИС (отчет): М., ЦГЭ.

20. Гурвич И. И. Сейсмическая разведка. Изд. 2-е, перераб. и доп. - М., Недра, 1970.

21. Гурвич И. И., Боганик Г. Н., Сейсмическая разведка. 3 изд., перераб. - М., Недра, 1980.

22. Дойбель К., 1982, Исследование сейсмических сигналов в реальных средах с поглощением: Сборник докладов второго научного семинара стран-членов СЭВ по нефтяной геофизике, Том 1, Сейсморазведка, С.171-180.

23. Долгов Е. Ф., Шехтман Г. А. Определение сейсмических скоростей по вертикальному годографу без приведения к вертикали. Прикладная геофизика. М.: Недра, 1986. Вып. 114. С. 56-62.

24. Душутин А.К., Баранов К.В. Определение поглощения по данным вертикального сейсмического профилирования, Тезисы докладов V международной конференции "Новые идеи в науках о Земле", Москва, 2001, С. 380.

25. Душутин А.К., Баранов К.В. Совершенствование обработки данных ВСП, Тезисы докладов V международной конференции "Новые идеи в науках о Земле", Москва, 2001, С. 379.

26. Душутин А.К., Баранов К.В., Бикеев B.C., Симарчук А.В. Оперативная обработка и контроль качества полевых работ ВСП, Сборник докладов научно-практической конференции, посвященной 5-летнему юбилею ООО «КогалымНИПИнефть», Когалым, 2001.

27. Жданович В. В., Ознобихин Ю. В., Монастырев Б. В., 1997, Изучение и компенсация искажающих свойств верхней части разреза в сейсморазведке: Геофизика, 6, 22-36.

28. Канарейкин Б. А., Седлецкий Н. Н. Результаты решения прямой кинематической задачи скважинной сейсморазведки для некоторых сложнопостроенных сред // Геология и геофизика. 1983. - № 5. - С. 144-148.

29. Катников А. П. Стандартизация условий возбуждения на сейсмограммах вертикального сейсмического профилирования // Геофизические исследования на нефть и газ в Узбекистане. Ташкент, 1987. - С. 66-74.

30. Кашик А. С., Кивелиди В. X., Шевченко А. А., Шакиров Р. А. Патент РФ № 2065182 «Способ пространственной сейсморазведки». Приоритет 22 апреля 1994 года.

31. Козырев В. С., Королев Е. К., 1993, Интерактивная методика коррекции статических поправок для условий сложного строения ВЧР: Геофизика, 3, С.13-19.

32. Кондратьев О. К., 2000, Математика в геофизике: Геофизика, 2, С.6-15.

33. Лаврик А. С., Геништа А. Н., 2001, Интерпретационный подход к учету неоднородностей ВЧР при обработке 2D- и 3D- сейсморазведки ОГТ на территории Западной Сибири: Геофизика, 1, С. 61-63.

34. Лапин С. И. Годографы ВСП при наблюдениях в анизотропных средах // Геология и геофизика. -1978. № 8. - С. 97-104.

35. Левянт В.Б., Шустер В.Л., 2001, Некоторые результаты подтверждаемости и геолого-экономической эффективности применения трехмернойсейсморазведки ЗД. Сборник тезисов докладов научно-практической конференции «ГЕОМОДЕЛЬ -2001».

36. Милашин В. А., 2000, Почему необходимо переобрабатывать и переинтепретировать данные сейсморазведки, полученные в Западной Сибири в предыдущие годы: Геофизика, 4, С. 26-28.

37. Митяев С. В. Решение обратной кинематической задачи при интерпретации данных скважинной сейсморазведки// Геология и геофизика. 1983. - № 6. - С. 136-139.

38. Петрашень Г. И. Распространение объемных волн и методы расчета волновых полей в анизотропных упругих средах. Л., Наука, 1984.

39. Притчетт У. Получение надежных данных сейсморазведки. М., Мир, 1999.

40. Пузырев Н. Н. Измерение сейсмических скоростей в скважинах. М., Гостоптехздат, 1957.

41. Рапопорт М. Б. Вычислительная техника в полевой геофизике: Учеб. Для вузов. 2 изд., перераб. и доп. - М., Недра, 1993.

42. Реклейтис Г., Рейвиндран А., Рэгсдел К. Оптимизация в технике: Пер. с англ. -М.: Мир, 1986.

43. Рожков А.С. Способ автоматического определения точек излома вертикального годографа ВСП. ВИНИТИ РАН, № 2268-В99, 1999.

44. Савин И. В., Шехтман Г. А. Обратная кинематическая задача ВСП для сред с неплоскими границами раздела. Прикладаная геофизика. М.: Недра, 1995. -Вып. 129.-С. 34-46.

45. Справочник геофизика, т. IV, "Сейсморазведка", под ред. И. И. Гурвича и В. П. Номоконова, М., Недра, 1966.

46. Старобинец М., Васильков А., 1987, Расчет статических поправок по годографам волн в области первых вступлений: Исследования и разработки в области нефтяной геофизики в странах членах СЭВ, С.317-323.

47. Табаков А.А. Оценка геолого-геофизического разреза ниже забоя разведочной скважины по данным ВСП. Тр. СредазНИИ геол. и минер. Сырья. 1975. -Вып. 18.

48. Табаков А.А., Баранов К.В. и др. Патент РФ №2001611778 «Интегрированная система обработки и интерпретации геолого-геофизических данных» («ЮНИВЕРС») от 24 декабря 2001г.

49. Табаков А.А., Солтан И.Е., Боков П.М., Душутин А.К., Баранов К.В. Оценка трехмерной модели среды по трехкомпонентным наблюдениям ВСП, Тезисы докладов научно-практической конференции «Состояние и перспективы развития метода ВСП», Москва, 2001, С. 73.

50. Теплицкий В. А. Применение скважинной сейсморазведки для изучения структур в нефтегазовых районах. М., Недра, 1973.

51. Тонев К. С. Повышение разрешающей способности сейсмического метода путем комбинированного использования поверхностных и скважинных сейсмических исследований// Тр. XXX Междунар. геофиз. симпоз. Сер. А.- М., 1985.-Ч. 3-С. 18-38.

52. Трапезникова Н.А., Птецов Р.С., Ворожцов JI.H., Романов Ю.А., Птецов С.Н., 2000, Трехмерная сейсморазведка оценки точности и эффективности. Сборник тезисов докладов научно-практической конференции «ГЕОМОДЕЛЬ -2000».

53. Ференци В.Н., Яковлев И.В., Барков А.Ю., 2001, Редактирование гармонических и всплесковых шумов в записях ВСП: Тезисы докладов научно-практической конференции «Состояние и перспективы развития метода ВСП», Москва, 2001, С. 69.

54. Форсайт Дж., Малькольм М., Моулер К. Машинные методы математических вычислений. М., Мир, 1980.

55. Шехтман Г. А., Зернов А. Е. Развитие метода ВСП для поиской и разведки месторождений углеводородов в условиях сложнопостроенных сред // Абстракты и доклады 34-го Международного геофизического симпозиума. -Будапешт, 1989.

56. Chiu S. К. L., Kanasewich Е. R., Phadke S. Three dimensional determination of structure and velocity by seismic tomography // J. Geophysics. Vol.51, 1986, P. 15591571.

57. Chiu S. К. Г. and Stewart R. R. Tomographic determination of three-dimensional seismic velocity structure using well logs, Vertical seismic profiles, and surface seismic data // J. Geophysics. Vol.52, № 8, 1987, P. 1085-1098.

58. Hauge P. S. Measurement of attenuation from vertical seismic profile// Geophysics. 1981. Vol. 46. N 11. P. 1548-1558.

59. Lines L. R., Bourgeois A., Covey J. D. Traveltime inversion of offset vertical seismic profiles a feasibility study // J. Geophysical prospecting, 1984. - V. 49. - P. 250264.

60. More J.J., "The Levenberg-Marquardt Algorithm: Implementation and Theory," Numerical Analysis, ed. G. A.Watson, Lecture Notes in Mathematics 630, Springer Verlag, P. 105-116, 1977.

61. Murat M.E. and Rudman A.J. Automated first arrival picking: a neural network approach. Geophysical prospecting, v.40, P.587-605, 1992.

62. Rossi G., Vesnaver A., Petersen S. A., Anisotropy detection by tomography and polarisation analysis in a 3D three-component VSP. First Break, April, 2001. Vol. 19. Number 4. P. 191-200.

63. Schuster G. Т., 1988, An analytic generalized inverse for common-depth-point and vertical seismic profile traveltime equations: Geophisics, 53, P.314-325.

64. Schuster G. Т., Johnson D. P., and Trentman D. J., 1988, Numerical verification and extension for an analytic generalized inverse for common-depth-point and vertical-seismic-profile travel equations: Geophisics, 53, P.326-333.

65. Stainsby S. D., Worthington M. H. Q estimation seismic profile data and anomalous variation in the central North Sea. Geophisics, Vol. 50, № 4, 1885, P. 615-626.

66. Stewart R. R. VSP interval velocities from traveltime inversion // J. Geophysical prospecting, 1984.-V. 32. P. 608-628.