CC BY
66
Дю. При ДПФ для уменьшения шага по частоте необходимо увеличивать длительность сигнала. В такой ситуации анализируемый участок сейсмотрассы может быть дополнен нулями. Однако значения ФЧХ в этом случае будут коррелированными. Кор-релированность значений ФЧХ, как показано выше, приведет к изменению весовых коэффициентов в алгоритмах обработки, которые на практике найти сложно. Однако, учитывая, что сама процедура обработки не меняется, переход к алгоритмам с равновесной обработкой позволяет избежать возникших трудностей и получить оценки временного положения сигналов с достаточно высокой точностью.
СПИСОК ЛИТЕРАТУРЫ
1. Иванченков В.П., Кочегуров А.И. Определение временного положения сейсмических сигналов по оценкам их фазочастотных характеристик // Геология и геофизика. - 1988. - № 9. -С. 77-83.
2. Иванченков В.П., Вылегжанин О.Н., Орлов О.В., Кочегуров А.И. Методы фазочастотного анализа волновых полей и их применение в задачах обработки данных сейсморазведки // Известия Томского политехнического университета. - 2006. -Т 309. - № 7. - С. 65-70.
3. Иванченков В.П., Кочегуров А.И. Фазочастотные алгоритмы оценки местоположения пространственно- временных сигналов в условиях априорной неопределенности // Известия высших учебных заведений. Физика. - 1995. - Т. 38. - № 9. -С. 100-104.
4. Поиск, обнаружение и измерение параметров сигналов в радионавигационных системах / Под ред. Ю.М. Казаринова. -М.: Советское радио, 1975. - 296 с.
При проведении сейсморазведочных работ методом вертикального сейсмического профилирования (ВСП) скорости продольных волн оценивают по времени первого вступления на наблюденном волновом поле с ближнего пункта возбуждения (ПВ) [1, 2]. Так как анализ проводится по одно-
Проведенный анализ способов развертывания ФЧХ показал, что каждый из них наряду с несомненными достоинствами, имеет и недостатки. В целом, предпочтение следует отдать методу Шафера, т. к. он прост в реализации, а возникающие погрешности при развертывании ФЧХ в анализируемой полосе частот можно контролировать путем анализа исходной записи.
Таким образом, фазочастотные методы обеспечивают достаточно высокую точность оценок временного положения сигналов даже при наличии корреляции в выборке ФЧХ сейсмической записи.
5. Долгополов Д.В., Пасторов А.И. О разделении двух наложившихся импульсов // Применение ЭВМ в сейсмологической практике. Методические работы ЕССН. - М.: Наука, 1985. -С. 86-91.
6. Кочегуров А.И., Быстров В.Н. Определение временного положения сложных сигналов в среде с дисперсией и поглощением // Известия высших учебных заведений. Радиоэлектроника. -2002. - Т. 45. - № 3-4. - С. 50-54.
7. Гольдин С.В. Смещение нулей и экстремумов сейсмических сигналов под воздействием помех // Геология и геофизика. -1964. - № 10. - С. 130-144.
8. Tribolet J.M. A new phase unwrapping algorithm // IEEE Transaction on acoustics, speech and signal processing. - 1977. - V. 25 (2). - P. 170-177.
Поступила 13.11.2009 г.
кратному наблюдению, полученному с помощью перестановки приемников и многократного возбуждения, и наблюденное поле осложнено помехами, оценка интервальных скоростей обычно обладает значительными погрешностями. Последующие процедуры обработки полей ВСП (приведение
УДК 550.8.053:519.2
КОРРЕКЦИЯ СКОРОСТНОГО ЗАКОНА ПО ДАННЫМ НЕПРОДОЛЬНОГО ВЕРТИКАЛЬНОГО СЕЙСМИЧЕСКОГО ПРОФИЛИРОВАНИЯ
Д.Ю. Степанов, М.С. Речкин
Томский политехнический университет E-mail: w00x@sibmail.com
Предложен алгоритм коррекции скоростного закона по данным непродольного вертикального сейсмического профилирования. Рассмотрены модели ошибок в определении статических поправок, показано, что данный алгоритм позволяет минимизировать влияние погрешности определения статических поправок и тем самым повысить точность оценки интервальных скоростей.
Ключевые слова:
Вертикальное сейсмическое профилирование, интервальные скорости, статическая поправка.
Key words:
Vertical seismic profiling, interval velocities, static correction.
поля к вертикали, построение временных и глубинных разрезов и т. п.) предъявляют повышенные требования к используемому скоростному закону. Поэтому на практике интервальные скорости не применяют, а используют пластовые скорости, полученные путем усреднения интервальных скоростей или из усредненных времен по глубине. Тем самым повышается точность скоростного закона, но теряется разрешенность. Улучшение качества оценки без потери разрешенности можно добиться путем разработки алгоритма коррекции скоростей, учитывающего результаты наблюдений не только с ближнего ПВ (продольное ВСП), но и с удаленных ПВ (непродольное ВСП) [1].
Система наблюдений прямого ВСП состоит из каротируемой скважины, пунктов приема (ПП) и пунктов возбуждения (рис. 1). Пункты приема располагаются в скважине на глубинах К1, Н2, ..., Ня. Пункты взрыва выносят на расстояния от каротируемой скважины и закладывают в неглубокую скважину (8-15 м). При удалении ПВ, не сопоставимым с глубиной забоя скважины, пункт возбуждения называют ближним, остальные - удаленными.
При первичной обработке материалов ВСП осуществляют ввод статической поправки, позволяющей получить сейсмограмму ВСП для одинаковых условий возбуждения [2]. Определение величины поправки осуществляется по записи контрольного канала и в ряде случаев не учитывает все изменения в системе наблюдений. Исходя из этого, наблюденное время прихода прямой продольной волны можно записать как:
Г(Н)= ¡(К)+Д/(К), (1)
где ¡н(К) - наблюденное время, ¡(К) - теоретическое время прихода волны, Д(К) - величина, определяющая отклонение наблюденного времени от теоретического, в общем случае имеет случайный характер. Пусть величина Д/(К) представима в виде: Д1(Н)=х(Н)+8(Н), где т(К) - ошибка ввода статической поправки, 8(Н) - случайная величина с нулевым математиче-
ским ожиданием, обусловленная ошибками измерений. Далее будем считать, что т(К) не зависит от глубины К.
Рассмотрим проведение ВСП с М пунктами возбуждения, удаленными от скважины на расстояния: Ьь Ьг, ..., ЬМ, и N пунктами приема. Пусть в результате получено М наблюденных времен первых вступлений прямых продольных волн, которые согласно (1) можно записать в виде:
¡кв(К)=1к(К)+тк, к=1,2, ..., М, где ¡к(К) - теоретическое время прихода на глубину К с к-го ПВ, тк - ошибка в определении статической поправки для к-го ПВ.
Для анализа ошибок статических поправок будем использовать матрицу невязок наблюденных данных и принятой скоростной модели среды. Построение такой матрицы предлагается проводить по следующей схеме:
1.
2.
Решается обратная задача: рассчитываются скоростные законы Им для каждого пункта возбуждения по наблюденным временам первых вступлений [3].
Скоростной закон V принимается как базовый (/=1).
Решается прямая задача, используя базовый закон: находятся теоретические времена ¡к(К) для всех пунктов возбуждения.
В /-й столбец матрицы Т записываются величины, определяемые выражением
Т, = Ní^Ъ^&)], I = 1,2,
N п=1
- средние невязки наблюденного и теоретического времени для к-ого ПВ.
5. Пункты 2-4 повторяются, в предположении истинности других скоростных законов.
Матрицу невязок Т можно использовать для коррекции наблюденных времен и, соответственно, скоростных законов с разных ПВ. Для этого предлагается следующая последовательность действий:
ПВ4
/а.
\ ¿4
ПВ\
&
£¿1
пв5
/
/
/
/
/
/
ь
3
ПВ3
Ж'
ь
пв2
Рис. 1. Пример схемы расположения ПВ при проведении ВСП
3
1. Рассчитывается матрица невязок Т по наблюденным временам.
2. Выбираются 2 столбца и АШх матрицы Т по следующим критериям
Ап1ах = а*ё тах
Ае[ 1, М ]
1 М
—Ет
М -1 £
¡*к
1 М
м -1£
3. Находятся скорректированные наблюденные времена:
¡к"<ШШ)(К) = ¡к"(К) - ТкАш|п,
¡кН<1ШК)(К) = 4"(К - ТкАшах, и по скорректированным временам рассчитываются новые матрицы невязок: Т^ и Тшах.
4. Из выбранных столбцов АШп и Ашах выбирается тот, использование которого при коррекции времен приводит к меньшей норме матрицы Т. За норму ||7Ц можно принять среднее значение модулей элементов матрицы.
5. Коррекция времен (пункты 2-4) происходит до выполнения одного из двух условий:
• условие достижения требуемой нормы: ||ТЦ<£, где е - порог;
• условие неэффективности дальнейших итераций: ||Т)||-||Т)-1||>0, где Т - матрица невязок, рассчитанная на у-ой итерации.
В результате выполнения вышеописанного алгоритма, оценки интервальных скоростей по скорректированным временам будут содержать только ошибку, обусловленную ошибками измерения времени первого вступления. Для анизотропной и горизонтально-слоистой среды следует ожидать, что данные оценки являются несмещенными. Поэтому коррекция скоростного закона может быть выполнена за счет усреднения по ансамблю:
1 М
укор (И) = —Еук°р (И),
Мк=1
где Уккор(Н) - интервальная скорость, рассчитанная по скорректированному времени к-го ПВ.
Данный алгоритм использован и апробирован на модельных материалах. Параметры моделей выбирались из условий проведения работ ВСП на территории Западной Сибири, а в качестве скоростного закона выбраны интервальные скорости, полученные на реальных материалах одной из скважин Томской области (далее - модельная скорость). Система наблюдения задавалась вертикальной скважиной с глубиной забоя 3220 м, расстановкой сейсмоприемников равномерно по глубине с шагом 10 м и пятью пунктами возбуждения с удалениями от устья скважины на расстояния: Х1=150 м, Х2=500 м, Х3=1000 м, Х5=1500 м.
На первом этапе проводились исследования алгоритма оценки скоростного закона с целью определения влияния ошибок ввода статических поправок и удаления ПВ. Теоретические времена ¡к(К) рассчитывались по модельным скоростям, а наблюденные времена ¡н(К) согласно (1) [3]. На рис. 2, а, представлены графики погрешностей интервальных скоростей при различных ошибках в статических поправках и постоянном удалении ПВ Х=500 м. На рис. 2, б, представлены погрешности интервальных скоростей для различных удалений ПВ при одинаковой ошибке в статической поправке т=8 мс. Как видно из результатов эксперимента ошибки в определении статической поправки приводят к нелинейному искажению интервальных скоростей. Наибольшие отклонения наблюдаются в верхней части разреза и усиливаются с уменьшением удаления ПВ, а с глубиной влияние статических поправок в среднем уменьшается.
На втором этапе исследований ставилась задача проверки работоспособности алгоритма коррекции наблюденных времен на различных моделях ошибок ввода статических поправок. Ошибка вво-
Отклонение, м/с Отклонение, м/с
а б
Рис. 2. Отклонение оценок интервальной скорости от модельной скорости в зависимости от: а) величины т; б) удаления ПВ
да статической поправки определялась согласно табл. 1 для каждого ПВ и модели ошибки.
Таблица 1. Модели ошибок ввода статических поправок
№ модели ошибки Ошибка в статических поправках для ПВ, мс/ Среднее отклонение оценок скоростей, м/с
L1 L2 Lз L4 L5
1 4/130 0/0 0/0 0/0 0/0
2 0/0 4/72 4/37 4/31 4/25
3 4/130 -4/72 0/0 -4/31 4/25
4 -4/130 -4/72 0/0 -4/31 -4/25
5 0/0 -2/35 8/73 -2/15 -4/25
6 2/66 -2/35 4/37 -8/64 4/25
Отклонения скорректированных времен, мс
L1 L2 Lз и и5
1,8 -2 2,4 -2,5 2,8
Скорость, м/с
1400 1600 1800
§ 200 <в
§ 250
I
С 300
£
■ Модельная -Результирующая Усредненная Интервальная ПВ1
Рис. 3. Модельная и интервальные скорости
На третьем этапе исследовалась помехоустойчивость алгоритма коррекции скоростного закона. Случайная величина 8(Н) задавалась датчиком случайных чисел с равномерным законом распределе-
ния на интервале [-0,5; 0,5] мс. Принятая модель ошибок ввода статических поправок представлена в табл. 3.
Таблица 3. Модельные ошибки в оценке статических поправок
Ошибка в статических поправках для ПВ, мс / Среднеквадратическая ошибка оценок скоростей, м/с
и А и-3 А А5
16/29 24/14 12/2 18/25 0/21
В результате использования алгоритма коррекции наблюденных времен для моделей № 1-5 удалось скорректировать времена с точностью алгоритма решения прямой задачи [3]. Таким образом, разработанный алгоритм позволяет скорректировать наблюденные времена при отсутствии ошибок хотя бы в одном из наблюденных времен. В шестой модели заложено условие наличия ошибок ввода статической поправки в разрезы всех ПВ. Для данной модели достичь вычислительной точности не удалось, однако для удаленных ПВ удалось уменьшить отклонения наблюденных времен от теоретических (табл. 2).
Таблица 2. Погрешности коррекции времен для модели № 6
На рис. 3 показаны графики модельной скорости, усредненной оценки интервальных скоростей до коррекции и результирующего скоростного закона (усреднение после коррекции) в верхней части разреза. На рис. 4 приведены графики отклонений расчетных скоростей от модельной. Усреднение оценок скоростей без коррекции наблюденных времен приводит к получению смещенной оценки. Из рисунков видно, что после коррекции наблюденных времен с помощью представленного алгоритма результирующая скорость является несмещенной оценкой.
Отклонение, м/с -100 -60 -20 20 60 100
200
250
I 3
I
С 300
Усредненной Результирующей
Рис. 4. Отклонение оценок скоростей от модельной скорости
Таким образом, постоянное отклонение наблюденных времен от теоретических влечет смещение оценки скоростного закона. Случайное отклонение времен влечет за собой случайную неоднородную ошибку оценивания скорости. Неоднородность данной ошибки обусловлена рекуррентным алгоритмом оценки скоростей [3]. Разработанный алгоритм коррекции наблюденных времен позволяет минимизировать влияние ошибки ввода статической поправки и, тем самым, получить несмещенные оценки интервальных скоростей, а алгоритм коррекции скоростного закона - уменьшить разброс скорости, обусловленный случайной величиной 8(Н).
0
50
2000
0
50
350
100
400
150
450
500
350
400
СПИСОК ЛИТЕРАТУРЫ
1. Гальперин Е.И. Вертикальное сейсмическое профилирование. 2-е изд., доп. и перераб. - М.: Недра, 1982. - 344 с.
2. Сейсморазведка. Справочник геофизика / Под ред. И.И. Гур-вича, В.П. Номоконова. - М.: Недра, 1981. - 464 с.
3. Речкин М.С. Алгоритм оценки скоростного закона по данным ВСП // Молодежь и современные информационные технологии: Сбо. трудов VII Всеросс. научно-практ. конф. студентов, аспирантов и молодых ученых - Томск, 25-27 февраля 2009. -
Ч. 1. - Томск: Изд-во «СПб Графикс», 2009. - С. 179-180.
Поступила 29.10.2009
УДК 622.276.05-192:519.6
АНАЛИЗ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ДВИГАТЕЛЬНОЙ СЕКЦИИ ВИНТОВОГО ЗАБОЙНОГО ДВИГАТЕЛЯ В ПРОГРАММНОМ КОМПЛЕКСЕ ANSYS
К.В. Сызранцева, В.А. Аришин
Тюменский государственный нефтегазовый университет, г. Тюмень E-mail: kv.syzr@gmail.com
Проведен прочностной конечно-элементный анализ двигательной секции винтового забойного двигателя в программном комплексе ANSYS. Подробно описаны все этапы анализа. Приведены результаты оценки напряженно-деформированного состояния ротора и статора, а также картины распределений контактных давлений между ними, необходимые для оценки работоспособности двигателя. Показано изменение распределений контактных давлений для двигателей с изношенным статором, приводящее к нарушению условия герметичности рабочих камер.
Ключевые слова:
Винтовой забойный двигатель, компьютерный анализ, метод конечных элементов, ANSYS, напряженно-деформированное состояние, контактные давления.
Key words:
Screw downhole motor, computer analysis, finite element method, ANSYS, stress-strain state, contact pressures.
Бурение скважин необходимо для изучения геологической структуры и свойств горных пород, а также для добычи жидких и газообразных полезных ископаемых. Процесс бурения заключается в постепенном углублении ствола скважины путем разрушения горных пород долотом с одновременным выносом жидкостью или газом выбуренной породы по затрубному пространству на поверхность.
Винтовые забойные двигатели (ВЗД) играют заметную роль в буровой технике и технологии последних десятилетий. С усложнением в последнее время буровых работ и работ по ремонту скважин (горизонтальное бурение, разбуривание дополнительных стволов в пробуренной скважине) и использованием технологий с применением гибкой трубы, роль винтовых забойных двигателей ввиду своих уникальных характеристик постоянно возрастает [1].
Уникальность характеристик ВЗД заключается в том, что при своей простоте конструкции и относительно небольших размерах они обладают высоким крутящим моментом на низких оборотах, а это необходимо для более качественного и экономически выгодного бурения. К тому же ВЗД имеют высокий КПД по сравнению с другими гидравлическими двигателями.
При всех своих достоинствах винтовые забойные двигатели обладают одним существенным недостатком - быстрым износом рабочих органов и, соответственно, недостаточно большим ресурсом работы.
Конструкция двигателя предусматривает несколько основных частей: двигательная секция, шпиндельная секция, переливной клапан, гибкий вал, переводники.
Двигательная секция состоит из ротора и статора, поперечные сечения которых изображены на рис. 1. Стальной статор имеет внутри привулкани-зированную резиновую обкладку с винтовыми зубьями, на стальном роторе нарезаны наружные зубья, причем число зубьев ротора на единицу меньше числа зубьев статора. Специальный профиль зубьев ротора и статора обеспечивает непрерывный контакт между собой, образуя на длине шага статора единичные рабочие камеры [2].
Ротор, совершая планетарное движение, поворачивается по часовой стрелке, в то время как геометрическая ось ротора перемещается относительно оси статора против часовой стрелки. За счет разности в числах зубьев ротора и статора переносное движение редуцируется в абсолютное с передаточным числом, равным числу зубьев ротора, что обеспечивает сниженную частоту вращения и высокий крутящий момент на выходе.
