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

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

Повний опис

Збережено в:

Бібліографічні деталі
Дата:	2015
Автори:	Тульчинский, П.Г., Роганов, В.Ю., Лавренюк, А.Н., Лавренюк, С.И.
Формат:	Стаття
Мова:	Russian
Опубліковано:	Інститут кібернетики ім. В.М. Глушкова НАН України 2015
Назва видання:	Компьютерная математика
Теми:	Системный анализ
Онлайн доступ:	http://dspace.nbuv.gov.ua/handle/123456789/168379
Теги:	Додати тег Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:	Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:	Оптимизированный метод обработки данных при определении азимутов трещиноватостей в волновом акустическом каротаже скважин / П.Г. Тульчинский, В.Ю. Роганов, А.Н. Лавренюк, С.И. Лавренюк // Компьютерная математика. — 2015. — № 2. — С. 40-50. — Бібліогр.: 12 назв. — рос.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine

id	irk-123456789-168379
record_format	dspace
spelling	irk-123456789-1683792020-05-02T01:28:17Z Оптимизированный метод обработки данных при определении азимутов трещиноватостей в волновом акустическом каротаже скважин Тульчинский, П.Г. Роганов, В.Ю. Лавренюк, А.Н. Лавренюк, С.И. Системный анализ Предложен оптимизированный метод обработки данных при определении азимутов трещиноватостей в волновом акустическом каротаже скважин. Этот подход дает возможность получать точный результат и уменьшить время работы программы. Приведены теоретические расчеты и результаты обработки реальных данных волнового акустического каротажа, проведенного прибором АВАК-11. Запропоновано оптимізований метод обробки даних при визначенні азимутів тріщинуватості у хвильовому акустичному каротажі свердловин. Цей підхід дає можливість отримувати точний результат і зменшити час роботи програми. Наведено теоретичні розрахунки та результати обробки реальних даних хвильового акустичного каротажу, проведеного приладом АВАК-11. Optimized data processing method in the azimuth fracturing determination in the wave acoustic borehole logging is presented. This approach provides obtaining accurate results and decreases the program runtime. Theoretical calculations and results of the real wave acoustic borehole logging data processing conducted by the «АВАК-11» device are presented. 2015 Article Оптимизированный метод обработки данных при определении азимутов трещиноватостей в волновом акустическом каротаже скважин / П.Г. Тульчинский, В.Ю. Роганов, А.Н. Лавренюк, С.И. Лавренюк // Компьютерная математика. — 2015. — № 2. — С. 40-50. — Бібліогр.: 12 назв. — рос. 2616-938Х http://dspace.nbuv.gov.ua/handle/123456789/168379 004.4, 004.9, 004.42, 004.67, 004.915, 004.427 ru Компьютерная математика Інститут кібернетики ім. В.М. Глушкова НАН України
institution	Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection	DSpace DC
language	Russian
topic	Системный анализ Системный анализ
spellingShingle	Системный анализ Системный анализ Тульчинский, П.Г. Роганов, В.Ю. Лавренюк, А.Н. Лавренюк, С.И. Оптимизированный метод обработки данных при определении азимутов трещиноватостей в волновом акустическом каротаже скважин Компьютерная математика
description	Предложен оптимизированный метод обработки данных при определении азимутов трещиноватостей в волновом акустическом каротаже скважин. Этот подход дает возможность получать точный результат и уменьшить время работы программы. Приведены теоретические расчеты и результаты обработки реальных данных волнового акустического каротажа, проведенного прибором АВАК-11.
format	Article
author	Тульчинский, П.Г. Роганов, В.Ю. Лавренюк, А.Н. Лавренюк, С.И.
author_facet	Тульчинский, П.Г. Роганов, В.Ю. Лавренюк, А.Н. Лавренюк, С.И.
author_sort	Тульчинский, П.Г.
title	Оптимизированный метод обработки данных при определении азимутов трещиноватостей в волновом акустическом каротаже скважин
title_short	Оптимизированный метод обработки данных при определении азимутов трещиноватостей в волновом акустическом каротаже скважин
title_full	Оптимизированный метод обработки данных при определении азимутов трещиноватостей в волновом акустическом каротаже скважин
title_fullStr	Оптимизированный метод обработки данных при определении азимутов трещиноватостей в волновом акустическом каротаже скважин
title_full_unstemmed	Оптимизированный метод обработки данных при определении азимутов трещиноватостей в волновом акустическом каротаже скважин
title_sort	оптимизированный метод обработки данных при определении азимутов трещиноватостей в волновом акустическом каротаже скважин
publisher	Інститут кібернетики ім. В.М. Глушкова НАН України
publishDate	2015
topic_facet	Системный анализ
url	http://dspace.nbuv.gov.ua/handle/123456789/168379
citation_txt	Оптимизированный метод обработки данных при определении азимутов трещиноватостей в волновом акустическом каротаже скважин / П.Г. Тульчинский, В.Ю. Роганов, А.Н. Лавренюк, С.И. Лавренюк // Компьютерная математика. — 2015. — № 2. — С. 40-50. — Бібліогр.: 12 назв. — рос.
series	Компьютерная математика
work_keys_str_mv	AT tulʹčinskijpg optimizirovannyjmetodobrabotkidannyhpriopredeleniiazimutovtreŝinovatostejvvolnovomakustičeskomkarotažeskvažin AT roganovvû optimizirovannyjmetodobrabotkidannyhpriopredeleniiazimutovtreŝinovatostejvvolnovomakustičeskomkarotažeskvažin AT lavrenûkan optimizirovannyjmetodobrabotkidannyhpriopredeleniiazimutovtreŝinovatostejvvolnovomakustičeskomkarotažeskvažin AT lavrenûksi optimizirovannyjmetodobrabotkidannyhpriopredeleniiazimutovtreŝinovatostejvvolnovomakustičeskomkarotažeskvažin
first_indexed	2025-07-15T03:10:18Z
last_indexed	2025-07-15T03:10:18Z
_version_	1837680839937228800
fulltext	40 Компьютерная математика. 2015, № 2 Предложен оптимизированный метод обработки данных при определении азимутов трещино- ватостей в волновом акустиче- ском каротаже скважин. Этот подход дает возможность полу- чать точный результат и уменьшить время работы про- граммы. Приведены теоретиче- ские расчеты и результаты обработки реальных данных вол- нового акустического каротажа, проведенного прибором АВАК-11. © П.Г. Тульчинский, В.Ю. Роганов, А.Н. Лавренюк, С.И. Лавренюк, 2015 УДК 004.4, 004.9, 004.42, 004.67, 004.915, 004.427 П.Г. ТУЛЬЧИНСКИЙ, В.Ю. РОГАНОВ, А.Н. ЛАВРЕНЮК, С.И. ЛАВРЕНЮК ОПТИМИЗИРОВАННЫЙ МЕТОД ОБРАБОТКИ ДАННЫХ ПРИ ОПРЕДЕЛЕНИИ АЗИМУТОВ ТРЕЩИНОВАТОСТЕЙ В ВОЛНОВОМ АКУСТИЧЕСКОМ КАРОТАЖЕ СКВАЖИН Введение. Определение параметров трещи- новатых зон – это важный этап при выполнении разведочных работ и изучении свойств резервуаров. Наличие трещиноватых зон приводит к появлению азимутально- зависимой анизотропии и изменению свойств сейсмических волн (поляризации, амплитуды и времени распространения), проходящих через эти зоны [1]. При наличии одной вертикальной системы трещин, разме- щенной в изотропной среде, возникает горизонтальная трансверсально-изотропная (ГТИ) среда. Квазипоперечная волна, распро- страняющаяся вертикально, расщепляется на быструю S1 и медленную S2 волны. Скорость этих волн зависит от их поляриза- ции. Быстрая волна имеет поляризацию вдоль системы трещин, а медленная – вдоль нормалей к плоскостям трещин [2]. Исполь- зуя эти волны, можно определить азимут трещиноватости [3 – 8]. Поскольку, медлен- ная волна отстает от быстрой, вектор поляризации волнового поля будет отличным от нуля на кросс компонентах XY и YX во всех случаях кроме случая, когда источник генерирует волну с поляризацией в плоско- сти трещиноватости или ей перпендикуляр- ной. Прямоугольная система координат, для которой ось OX расположена в плоскости трещиноватости или ей перпендикулярна, называется натуральной системой координат. ОПТИМИЗИРОВАННЫЙ МЕТОД ОБРАБОТКИ ДАННЫХ ПРИ ОПРЕДЕЛЕНИИ ... Компьютерная математика. 2015, № 2 41 Для определения направления трещин рассчитываются сигналы при различных поворотах системы координат вокруг вертикальной оси OZ и среди них выбирается тот поворот, для которого на кросс компонентах находится минимальная энергия. Повороты системы координат выполняются преобразова- нием Алфорда [4]. Для возбуждения квазипоперечных волн при волновом акустическом каротаже (ВАК) используются источники типа силы, направленные вдоль горизонтальных осей OX и OY. Два приемника располагаются на некотором расстоянии по оси OZ от источников и регистрируют колебания вдоль направлений OX и OY. В результате, для каждого положения прибора в скважине, регистрируются 4 трассы, XX, XY, YX и YY, которые содержат два типа квазипоперечных волн, распространяющихся с разными скоростями в около- скважинном пространстве. В данной работе представлен метод определения азимута трещиноватостей на основе поворотов Алфорда, адаптированный для зондов ВАК. Метод позво- ляет синтезировать сигналы в натуральной системе координат и определять азимут трещиноватости. Обработка данных обменных волн для определения азимутов вертикальных трещиноватостей. Обработка данных основана на одномерной модели, в которой источниками являются силы, направленные вдоль осей OX и OY, а синтезированные поперечные волны распространяются вертикально вниз. Единичные векторы вдоль осей OX и OY обозначим соответственно ˆ ˆ,x y . Источники генерируют сигналы ˆ( )x s ts x и ˆ( )y s ts y , которые распространя- ются в вертикальном направлении и регистрируются приемниками, находящи- мися на глубине z и направленными вдоль осей OX и OY. Предполагается, что среда содержит одну вертикальную систему трещин с азимутом , который требуется определить. Сигналы Xs и Ys порождают быстрые и медленные волны с поляризациями  cos , sin ,0  и  sin , cos ,0 ,   которые распространяются в вертикальном направлении со скоростями 1Sv и 2.Sv Следовательно, на расстоянии z от источника будут зарегистрированы волны, описанные соотношениями (1): 2 2 1 2 ˆ( ) cos sin ,XX S S z zt s t s t v v                    s x 1 2 ˆ( ) sin cos ,XY S S z zt s t s t v v                     s y 1 2 ˆ( ) sin cos ,YX S S z zt s t s t v v                     s x 2 2 1 2 ˆ( ) sin cos .YY S S z zt s t s t v v                    s y (1) П.Г. ТУЛЬЧИНСКИЙ, В.Ю. РОГАНОВ, А.Н. ЛАВРЕНЮК, С.И. ЛАВРЕНЮК 42 Компьютерная математика. 2015, № 2 После прохождения расстояния ,z между быстрой и медленной волнами возникает временной сдвиг 1 2/ /S St z v z v   и на кросс компонентах XY и YX при 1 2S Sv v и / 2k   появляется ненулевая энергия. Сигналы на этих компонентах должны быть одинаковы. Их отличие между собой свиде- тельствует о наличии помех или неприменимости одномерной модели. Если время t мало, что может произойти при маленьких расстояниях источник- приемник ,z либо при слабой анизотропии, когда 1 2 ,S Sv v сигналы на компонентах XX и YY становятся одинаковыми. В этом случае азимут трещиноватостей восстанавливается со значительными погрешностями. Имея сигналы ( ), ( ), ( ), ( ),XX XY YX YYt t t ts s s s возбужденные и зарегистрирован- ные вдоль осей OX и OY, можно рассчитать аналогичные сигналы ( ), ( ), ( ), ( )XX XY YX YYt t t t   s s s s для источников и приемников, расположенных вдоль осей системы координат, повернутой на угол  вокруг вертикальной оси. Очевидно, это будут сигналы, описываемые формулой (2), с заменами угла  на угол  и ортов ˆ ˆ,x y на орты  ˆ cos ,sin ,0 ,   x  ˆ sin ,cos ,0 .    y Следовательно,     1 2 / 0cos sin cos sin 0 /sin cos sin cos SXX XY SYX YY s t z vs s s t z vs s                              , (2) XX XY YX YY s s s s           XX XY YX YY s s s s                               1 2 cos sin / 0 cos sin , sin cos 0 / sin cos S S s t z v s t z v                         cos sin cos sin sin cos sin cos XX XYXX XY YX YYYX YY s ss s s ss s                                . (3) Выполняя умножения матриц, находим искомые соотношения, которые называются вращением Алфорда [4]: 2 2cos ( )sin cos sin ,XX XX XY YX YYs s s s s         2 2sin ( )sin cos cos ,YY XX XY YX YYs s s s s         2 2cos sin ( )sin cos ,XY XY YX YY XXs s s s s         2 2cos sin ( )sin cos .YX YX XY YY XXs s s s s         (4) ОПТИМИЗИРОВАННЫЙ МЕТОД ОБРАБОТКИ ДАННЫХ ПРИ ОПРЕДЕЛЕНИИ ... Компьютерная математика. 2015, № 2 43 Значению азимута  трещиноватости соответствует угол поворота , при котором достигается минимум энергии на компонентах , .XY YXs s  Если минимум энергии на этих компонентах достигается при разных углах, то исходные данные не соответствуют принятой модели и требуется детально разбираться с причиной такого расхождения.     max max min min 2 2 min , min . t t XY XY YX YX t t E s dt E s dt       (5) Предложенный подход определения направлений трещин проверен на модельных данных, полученных 2.5D моделированием [9]. Затем откорректиро- ванные и оптимизированные алгоритмы с программной реализаций были применены для обработки данных, полученных при волновом акустическом каротаже скважин. Оптимизированный метод обработки данных, полученных при волно- вом акустическом каротаже скважин прибором АВАК-11. В отличие от обычного акустического каротажа, волновой акустический каротаж – это метод, оперирующий с параметрами не только продольной волны, но и других упругих волн, распространяющихся в скважине. Скважинный прибор АВАК-11 предназначен для исследования акусти- ческим методом нефтяных и газовых скважин, с открытым стволом и обсаженных, для определения акустической анизотропии и модулей упругости пород. Скважинный прибор АВАК-11 является цифровым, т. е. волновые картины от приемников зонда оцифровываются в скважинном приборе и передаются на поверхность кодами. Прибор содержит 4 дипольных зонда: 2 длинных и 2 коротких, которые генерируют и принимают сигналы вдоль взаимно- перпендикулярных осей OX и OY. В процессе каротажа на жесткий диск регистратора записываются оци- фрованные волновые картины, синхронизированные с моментом излучения. Для решения приведенных в статье задач используются волновые поля с компонентами XX, XY, YX, YY по коротким и длинным зондам. На рис. 1 показан пример каротажных данных по коротким зондам. Сигнал имеет незначительные шумы, которые убираются предварительной обработкой. Описанные далее приемы обработки и алгоритмы реализованы в програм- мном комплексе «ГЕОПОИСК» [10]. В настоящее время, в практике волнового каротажа информативными волнами, т. е. нашедшими практическое применение при решении геолого- технических задач, являются три типа волн: – Р – продольная волна; – S – поперечная волна; – L – волна Лемба – Стоунли. П.Г. ТУЛЬЧИНСКИЙ, В.Ю. РОГАНОВ, А.Н. ЛАВРЕНЮК, С.И. ЛАВРЕНЮК 44 Компьютерная математика. 2015, № 2 РИС. 1. Файлы с компонентами XX, XY, YX, YY по коротким дипольным зондам прибора АВАК-11 Эти волны, при проведении измерений АВАК-11, регистрируются на компонентах XX, XY, YX, YY. Сигналы P – продольная волна, S – поперечная волна и L – волны Лемба – Стоунли показаны на рис. 2 и 3. РИС. 2. Первые вступления волн, зарегистрированные приемником: P – продольной, S – поперечной и L – Лемба – Стоунли ОПТИМИЗИРОВАННЫЙ МЕТОД ОБРАБОТКИ ДАННЫХ ПРИ ОПРЕДЕЛЕНИИ ... Компьютерная математика. 2015, № 2 45 РИС. 3. Первые вступления волн на временном срезе: P – продольной, S – поперечной и L – Лемба – Стоунли Одна из первоочередных задач обработки данных волнового каротажа – это определение первых вступлений волн S и L, так как эта обработка ведется в пределах временного окна [S..L]. Если это окно определить непра-вильно, то возникнут значительные погрешности в результате обработки. Временные окна [S..L] определяются отдельно для длинных [S1..L1] и коротких зондов [S2..L2]. Временные окна в программе Геопоиск можно найти в полуавтоматическом режиме, когда программа по одному из выбранных алгоритмов (режим слежения за интервалом, пороговый, режим слежения за фазами, упрощенный пороговый), выполняет расчет первых вступлений для указанных типов волн. Имеется возможность проанализировать полученные кривые S, L и, при необходимости, внести исправления. Затем, по откорректиро-ванным данным, выполняются вычисления параметров трещиноватых зон. Рассмотрим численный пример обработки данных дипольного зондирова- ния прибором АВАК-11 на конкретной глубине. На рис. 4 показаны графики зарегистрированных волн на компонентах XX, YY и YX. Из рис. 4 видно, что на компоненте YX содержится достаточно большой сигнал, что свидетельствует о наличии анизотропии на глубине размещения прибора 1456z  метров. П.Г. ТУЛЬЧИНСКИЙ, В.Ю. РОГАНОВ, А.Н. ЛАВРЕНЮК, С.И. ЛАВРЕНЮК 46 Компьютерная математика. 2015, № 2 РИС. 4. Компоненты зарегистрированного волнового поля, содержащие квазипоперечные волны Временной интервал обработки  1400, 1750 мс выбран таким, чтобы в нем на всех компонентах присутствовали сигналы одинаковой формы и отно- сились только к квазипоперечным волнам. Энергия на YX-компоненте, вычи- сленная для разных азимутов системы координат, показана на рис. 5. РИС. 5. Энергия на кросс-компоненте YX при различных азимутах системы координат ОПТИМИЗИРОВАННЫЙ МЕТОД ОБРАБОТКИ ДАННЫХ ПРИ ОПРЕДЕЛЕНИИ ... Компьютерная математика. 2015, № 2 47 Минимум энергии достигается при 071 .  Применяя поворот Алфорда для найденного угла 071 ,  мы добиваемся минимума энергии на кросс- компоненте YX. В результате, получаем волновое поле, показанное на рис. 6. РИС. 6. Компоненты волнового поля после поворота на угол 0φ 71 На XX компоненте находится медленная квазипопречная волна, а на YY-компоненте – быстрая квазипоперечная волна. Следовательно, азимут трещиноватости, в данном случае, равен 0 090 161 .     На кросс-компоненте YX сигнал имеет небольшую амплитуду, что свидетельствует о том, что найден- ная система координат является натуральной, т. е. в ней ось OY расположена вдоль трещин, а ось OX ей перпендикулярна. Оценивая время регистрации по второй фазе (времени смены положительной фазы сигнала на отрицательную), находим время 1500 мс для быстрой квазипоперечной волны и 1530 мс – для медленной квазипоперечной волны. Коэффициент сдвиговой анизотропии оценивается соотношением:  2 1 2 1 2 ,S S S S t t a t t    (6) где 1,st 2st  времена распространения быстрой и медленной волн, соответ- ственно. В приведенном примере, коэффициент сдвиговой анизотропии, рассчитан- ной по формуле (6), равен 0.02. Аналогичные вычисления проводятся для всех глубин. Для вычисления сдвиговой анизотропии, при известном угле поворота , выполняется вращение Алфорда для получения компонент быстрой (S1) и медленной волн (S2). П.Г. ТУЛЬЧИНСКИЙ, В.Ю. РОГАНОВ, А.Н. ЛАВРЕНЮК, С.И. ЛАВРЕНЮК 48 Компьютерная математика. 2015, № 2 Значения компонент быстрой и медленной поперечных волн рассчи- тываются по формулам: 2 21 cos ( ) ( )sin( )cos( ) sin ( ),i i i i iS XX XY YX YY        2 22 sin ( ) ( )sin( )cos( ) cos ( ),i i i i iS XX XY YX YY        (7) где i = 1..n, n – количество отсчетов волновой картины. По полученным данным S1 и S2, рассчитываются времена первых вступлений всех типов волн и времена 1st и 2st для формулы (6). Временное окно волновой картины имеет w = 768 точек, длина каротажа h > 1500 точек. Так как современные рабочие станции имеют несколько процессоров (CPU) и процессорных ядер (CORE), обработку нескольких компо- нент данных оптимально проводить параллельно. Используя подходы по опре- делению оптимального метода распараллеливания, аналогично приведенным в [11], было определено, что распараллеливание целесообразно выполнять по w в рамках временного окна обработки для получения промежуточного результата – файла после преобразования Алфорда, в котором разделены быстрая и мед- ленная волны. Для оптимизации расчета и получения ускорения использованы директивы OpenMP [12] для языка C++. OpenMP (MP – сокращение «multiprocessing») – это кросс-языковый, кросс-платформенный интерфейс прикладного программирова- ния (API) для организации параллелизма на CPU. Он существует с 1997 года, поддерживает языки Fortran, C и C++ , реализован в Windows и на ряде других платформ. Директивы OpenMP позволяют выполнить распараллеливание части программного модуля для получения ускорения вычислений за счет распре- деления вычислений на имеющиеся процессоры и процессорные ядра. Данный механизм является универсальным для разных платформ. Также отметим, что не требуется дополнительное программирование при использовании данного подхода на компьютерах с одним процессором и одним процессорным ядром. В таблице приведены результаты ускорения в зависимости от количества процессорных ядер для численного примера, описанного выше. ТАБЛИЦА. Ускорение счета в зависимости от ядер процессора Количество ядер Ускорение счета, раз 2 1.08 4 1.20 Результатом вычислений является угол между направлением поляризации быстрой волны и осью OX. Этот угол необходимо скорректировать на угол между осью OX и вертикальной плоскостью, проходящей через скважину. В результате получим угол азимутальной трещиноватости, т. е. угол между направлением поляризации быстрой волны и вертикальной плоскостью в диа- пазоне (0 – 180)°, по часовой стрелке, при виде сверху. Также результатом вычислений является коэффициент сдвиговой анизотропии. Результат можно отобразить, как в виде кривых, так и в виде диаграмм Шмидта (рис. 7). ОПТИМИЗИРОВАННЫЙ МЕТОД ОБРАБОТКИ ДАННЫХ ПРИ ОПРЕДЕЛЕНИИ ... Компьютерная математика. 2015, № 2 49 РИС. 7. Отображение результата вычислений в программе «ГЕОПОИСК» Полученные данные – это один из важных параметров трещино- ватых зон. Выводы. 1. Разработан эффективный и устойчивый метод определения азимута вер- тикальной системы трещин по данным волнового акустического каротажа. Метод позволяет разделить зарегистрированное волновое поле на быстрые и медленные квазипоперечные волны. 2. Выполнена программная реализация предложенного метода для обработ- ки реальных данных волнового акустического каротажа. 3. Выполнена оптимизация вычислений на многоядерных или многопроцес- сорных компьютерах. П.Г. Тульчинський, В.Ю. Роганов, А.М. Лавренюк, С.І. Лавренюк ОПТИМІЗОВАНИЙ МЕТОД ОБРОБКИ ДАНИХ ПРИ ВИЗНАЧЕННІ АЗИМУТІВ ТРІЩИНУВАТОСТІ У ХВИЛЬОВОМУ АКУСТИЧНОМУ КАРОТАЖІ СВЕРДЛОВИН Запропоновано оптимізований метод обробки даних при визначенні азимутів тріщинуватості у хвильовому акустичному каротажі свердловин. Цей підхід дає можливість отримувати точний результат і зменшити час роботи програми. Наведено теоретичні розрахунки та результати обробки реальних даних хвильового акустичного каротажу, проведеного приладом АВАК-11. P. Tulchinsky, V. Roganov, A. Lavreniuk, S. Lavreniuk OPTIMIZED DATA PROCESSING METHOD IN THE AZIMUTH FRACTURING DETERMINATION IN THE WAVE ACOUSTIC BOREHOLE LOGGING Optimized data processing method in the azimuth fracturing determination in the wave acoustic borehole logging is presented. This approach provides obtaining accurate results and decreases the program runtime. Theoretical calculations and results of the real wave acoustic borehole logging data processing conducted by the «АВАК-11» device are presented. П.Г. ТУЛЬЧИНСКИЙ, В.Ю. РОГАНОВ, А.Н. ЛАВРЕНЮК, С.И. ЛАВРЕНЮК 50 Компьютерная математика. 2015, № 2 1. Bakulin A., Grechka V., Tsvankin I. Estimation of fracture parameters from reflection seismic data. – Part 1: HTI model due to a single fracture set // Geophysics. – 2000. – 65, N 6. – P. 1788 – 1802. 2. Schoenberg M., Helbig K. Orthorhomic media: Modeling elastic wave behavior in a vertically fractured eath // Geophysics. – 1997. – 62, N 6. – P. 1954 – 1974. 3. Роганов В.Ю., Роганов Ю.В. Моделирование и использование полей обменных волн для определения азимутов трещиноватостей // Геофиз. журн. – 2011. – 33, № 2. – С. 64 – 79. 4. Alford R. Shear data in the presence of azimuthal anisotropy // 56th Ann. Internal. Mtg., Soc. Expl. Geophys. Expanded Abstracts. – 1986. – P. 476 – 479. 5. Carcione J.M. Wave fields in real media. Wave propagation in anisotropic, anelastic, porous and electromagnetic media. Vol. 38. – Amsterdam: Elsevier, 2007. – 515 p. 6. Dellinger J., Nolte B., Etgen J. Alford rotation, ray theory, and crossed-dipole geometry // Geophysics. – 2001. – 66, N 2. – P. 637 – 647. 7. Gaiser J. Application for vector coordinate system of 3-D converter-wave data // The Leading Edge. – 1999. – 18, N 11. – P. 1290 – 1300. 8. Winterstein D. Velocity anisotropy terminology for geophysicist // Geophysics. – 1990. – 55, N 8. – P. 1070 – 1088. 9. Tulchinsky V.G., Roganov Y.V. Acceleration of 2.5D elastic anisotropic modelling // Copenhagen, Denmark: Society of Petroleum Engineers – 74th European Association of Geoscientists and Engineers Conference and Exhibition, 2012/ 10. Красножон М.Д., Косаченко В.Д., Тульчинський В.Г. Комплексна інтерпретація матеріалів ГДС з використанням комп’ютерної технології «Геопошук» // Матеріали міжнародної науково-технічної конференції «Прикладна геологічна наука сьогодні: здобутки та проблеми». – К.: УкрДГРІ, 2007. – С. 104 – 105. 11. Лавренюк А.Н., Лавренюк С.И., Назаренко Е.В. Подход к оптимизации программы обработки больших объемов сейсмических данных на примере 3D миграции дуплексных волн // Компьютерная математика. – 2014. – № 1. – C. 69 – 75. 12. OpenMP //[Электронный ресурс] – http://openmp.org/. Получено 06.05.2015 Об авторах: Тульчинский Петр Григорьевич, кандидат физико-математических наук, старший научный сотрудник Института кибернетики имени В. М.Глушкова НАН Украины, Роганов Вячеслав Юрьевич, младший научный сотрудник Института кибернетики имени В.М. Глушкова НАН Украины, Лавренюк Алла Николаевна, кандидат технических наук, доцент Физико-технического института НТУ Украины «КПИ», Лавренюк Сергей Иванович, кандидат физико-математических наук, научный сотрудник Института кибернетики имени В.М. Глушкова НАН Украины.

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

Репозитарії

Схожі ресурси