Физическое моделирование распространения упругих колебаний в анизотропных средах для задач сейсморазведки

Размещено на http://www.allbest.ru/

Физическое моделирование распространения упругих колебаний в анизотропных средах для задач сейсморазведки

Д.А. Попов, Т.И. Чичинина, В. В. Половков, И.О. Корсакова, Д. Пинбо, Б.М. Каштан, В.Н. Троян Санкт-Петербургский государственный университет, Российская Федерация, 19904, Санкт-Петербург, Университетская наб., 7-9

Мексиканский нефтяной институт, Мексика, 07730, Мехико, Сан Бартоло Атепэуацан,

Циудад де Меxико, Ейе Централ Лазаро Царденас Норте, 152

Китайский нефтяной институт,

Китай, 266580, Циндао, Чангйианг Вест Роад, Хуангдао Дистрикт, 66

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

Установка позволяет производить измерения скоростей и амплитуд продольных волн при угле 0° и углах в интервале от 15 до 90° относительно оси симметрии горизонтально-слоистой модели с минимальным шагом 8,1'. В настоящий момент идет модернизация оборудования с целью расширения возможностей -- изготовления многокомпонентных преобразователей, которые дадут возможность регистрировать P, SV и SH-волны в каждой точке измерения. Было выполнено физическое моделирование распространения волн в искусственном образце трансверсально-изотропной среды, которая представляет собой модель трещиноватой породы. В качестве тестовой среды используется модель из пластин оргстекла с вертикальной осью симметрии бесконечного порядка. По оценке значений скоростей продольных волн, были вычислены параметры анизотропии Томсена и определена зависимость скорости продольных волн от угла между осью симметрии модели и направлением распространения луча. Полученные экспериментальные данные показывают хорошее соответствие теоретической модели. скорость анизотропия симметрия волна

Ключевые слова: сейсморазведка, анизотропия, физическое моделирование, параметры анизотропии, модель линейного проскальзывания, модель Шонберга, трансверсально-изотропная среда.

A new approach in physical modeling for velocity anisotropy studies

D.A. Popov¹, T.I. Chichinina², V.V. Polovkov¹, 0. Korsakova¹, D. Pinbo³, B.M. Kashtan¹, V.N. Troyan¹

St. Petersburg State University, 7-9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation Instituto Mexicano del Petroleo, Mexico City, Mexico, Eje Central Lazaro Cardenas Norte 152,

San Bartolo Atepehuacan, Ciudad de Mйxico. C.P.07730

China University of Petroleum, Qingdao, China,

No. 66, Changjiang West Road, Huangdao District, 266580

Special automated equipment was constructed for ultrasonic velocity measurements at various incidence angles of waves that is necessary for estimation of anisotropy parameters of rocks. The equipment includes a hydraulic press for simulation of vertical pressure in the

Earth at a range up to 40 MPa, and a computer-controlled positioning system that provides precise independent positioning of the ultrasonic source and receiver. This setup enabled us to perform reliable velocity measurements at the angle interval from 15° to 90° and an additional single measurement at 0°. The minimum step of a transducer position in the 15°-90° range is 8.1'. As an example for the anisotropy study, we performed physical modeling of wave propagation in an artificial sample of transversely isotropic (TI) symmetry, which represents a model of fractured rock. Using an experimental setup, we acquired P-wave velocity data for a plate-stack model with a vertical axis of symmetry formed by acrylic 1-mm-thick sheet plates. Using measured velocities, we estimated Thomsen's anisotropy parameters Ј and 6, and got dependence between the velocity of P-waves and the symmetry axis of the model. In our oncoming research, the upgraded Ultrasonic Measuring System will be supplemented by multicomponent transducers similar to those described in Chichinina et al. (2009), which allow simultaneous recording of P, SV and SH waves at each data point.

Keywords: seismic, anisotropy, physical modeling, anisotropy parameters, linear-slip model, Schoenberg's model, VTI.

Введение

Модели эффективных сред, разработанные Хадсоном (Hudson, 1980) и Шон-бергом (Schoenberg, 1980; 1983), отмечают начало исследования анизотропии горных пород, приуроченной к естественной трещиноватости. Модель линейного проскальзывания Шонберга (linearslip; LS) является одной из простых моделей изотропной среды с одним набором параллельных трещин (Schoenberg and Sayers, 1995). Ее преимуществом является то, что количество независимых элементов матрицы модулей упругости уменьшено с пяти до четырех. Из-за удобства использования и меньшего числа параметров эта модель нашла широкое применение в сейсморазведке. В течение последних 20 лет было проведено множество лабораторных ультразвуковых экспериментов для проверки модели линейного проскальзывания, например: (Hsuand Schoenberg, 1993; Gikand Bobrov, 1996; Chichininaet al., 2007; Chichininaet al., 2009; Far, 2011; Faret al., 2014). Эксперименты проводились с использованием пластинчатых моделей, чаще всего представляющих собой пакет пластин из оргстекла, находящийся под воздействием внешнего давления.

Анализ опубликованных данных говорит о том, что эти исследования содержат оценки скоростей для небольшого количества углов прохождения волны и довольно ограниченных интервалов значений сдавливающего усилия. Фактически Хсу и Шонберг получили данные о скорости только для шести углов между 0 и 90° в диапазоне давлений от 6 до 24 МПа с шагом 6 МПа (Hsuand Schoenberg, 1993). Фар измерял скорости распространения волн для аналогичной пластинчатой модели только для трех углов 0, 45 и 90° в диапазоне давлений от 3 до 14 МПа (Far, 2011). В ряде экспериментов (Gikand Bobrov, 1996; Chichininaet al., 2009) использовалось невысокое значение сдавливающего усилия (максимум 4 МПа).

Для полноценного анализа имеющихся теоретических моделей анизотропных сред (и как следствие -- повышения достоверности изучения реальных трещиноватых пород) необходимо путем измерений на известной модели получить полный тензор модулей упругости. Наиболее сложно оценить модуль С₁₃, который необходим для надежных оценок параметра анизотропии Томсена ф (Thomsen, 1986), используемого во многих сейсмических и геомеханических приложениях. Одной из целей настоящей работы является повышение точности оценки параметра анизотропии Томсена ф из акустических измерений на образцах путем создания специальной ультразвуковой измерительной системы.

Существует ряд факторов, влияющих на точность определения модуля С₁₃ и, следовательно, параметра ф по измерениям скоростей ультразвуковых волн в образцах, например: (Sarout et al., 2015; Vernik, 2016; Yanet al., 2016). Точность оценки параметра С₁₃ зависит от точности измерений фазовой скорости продольной (Р) волны Vp,распространяющейся в интервале углов между 0 и 90° (чаще всего измеряется Vp45при ф = 45°). Фактически ограниченный диапазон доступных углов прохождения усложняет надежную оценку пятого компонента С₁₃ тензора модулей упругости трансверсально-изотропной среды (transverselyisotropic; TI). В работе (Chichinina, 2017) дается обзор теоретических и эмпирических ограничений параметра С₁₃ в трансверсально-изотропных породах, имеющих вертикальную ось симметрии (vertical transverselyisotropic; VTI). В статье 2015 г. (Chichininaet al., 2015) обнаружено, что нижняя граница С_13тт, установленная в работе (Yan et al., 2016), совпадает с формулой C₁₃для модели линейного проскальзывания. Также были определены верхние и нижние физические границы Сіз в органических сланцах, представленные в (Chichininaand Vernik, 2018):

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

1. Методология

Для проведения измерений на пластинчатых моделях с горизонтальной слоистостью под углами 0 и 90° к вертикальной оси симметрии достаточно расположить пары преобразователей на верхней/нижней и боковых гранях модели соответственно, однако получение данных для других углов требует использования более сложной методики. В ранее выполненных работах приводятся примеры перемещения приемника упругих колебаний на боковой грани в вертикальном направлении (Hsuand Schoenberg, 1993) (рис. 1А) либо движения приемника по окружности в плоскости боковой грани (Gikand Bobrov, 1996) (рис. 1Б). В обоих случаях источник является неподвижным. Второй способ является более корректным, так как расстояние между источником и приемником и их относительное расположение остаются неизменными, что избавляет от необходимости учета влияния диаграммы направленностей преобразователей и позволяет получать в процессе эксперимента данные, пригодные для динамического анализа. Нами было выбрано логическое продолжение второго подхода, при котором и источник, и приемник перемещаются в противофазе по одинаковым круговым траекториям, что позволяет уменьшить минимальное значение угла между линией измерения и осью симметрии, а также увеличить длину пути сигнала, что повышает точность определения значений скоростей (рис. 1В).

В основе установки (рис. 2) находится гидравлический пресс с максимальным усилием сжатия 50 тонн с ручным насосом и манометром (цена деления 2 МПа) для контроля давления в гидросистеме. Соотношение площадей основания штока пресса и верхней плоскости модели составляет примерно 10/6,5; точность поддержания постоянного давления при работе -- около 0,25 МПа. Таким образом, минимальный шаг изменения сдавливающего усилия составляет 1,3 ± 0,16 МПа. На плотно зафиксированную раму основания установлена плита из 40-миллиметровой стали. Аналогичная плита меньшего размера помещена непосредственно под штоком пресса. На нижнюю стальную плиту установлено основание измерительной системы, изготовленное из 16-миллиметрового алюминия. К основанию прикручены два вертикальных элемента из конструкционного алюминиевого профиля, на боковых сторонах которого закреплены две независимых системы позиционирования. Перемещение источника и приемника колебаний (пары № 2) происходит при помощи шаговых двигателей относительно нулевого положения, определяемого при помощи щелевого оптического датчика (рис. 1). Каждая система позиционирования подключается к управляющей аппаратуре при помощи двух разъемов.

Рис. 1. Три варианта траекторий перемещения пары преобразователей источник -- приемник (№ 2; обозначены красным и синим цветами соответственно) на боковых плоскостях пластинчатой модели: а -- вариант (Hsu, Schoenberg, 1993); б-- вариант (Gik, Bobrov, 1996); в-- предлагаемая методология. Положение пары преобразователей № 1 (обозначены серым цветом) постоянно, S = 0°

Рис. 2. Измерительная установка: а -- внешний вид; б -- конструкция: 1 -- гидравлический пресс; 2 -- стальная плита; 3 -- фиксирующие стержни; 4 -- элементы из конструкционного алюминиевого профиля; 5 -- рама; 6 -- шаговый двигатель; 7 -- разъемы; 8 -- оптический датчик нулевого положения; 9 -- основание измерительной системы; 10 -- пара преобразователей № 2; 11 -- блок прижима преобразователя; 12 -- плечо с регулируемым радиусом вращения; 13 -- пластинчатая модель; 14 -- оптический флаг; 15 -- главная ось вращения; 16 -- алюминиевые блоки; 17 -- стальная плита

Пластинчатая модель размером 200 х 50 х 200 мм собрана из листового оргстекла толщиной 1 мм. Размеры модели выбраны из соображений компромисса между максимальным значением сдавливающего усилия, расширения диапазона углов при измерениях и механической прочностью. Модель помещена между двух 20-миллиметровых алюминиевых блоков, служащих для распределения усилия и имеющих профрезерованные полости для расположения двух пар (для продольных и поперечных волн) ультразвуковых преобразователей, использующихся при измерениях вдоль оси симметрии модели (пара преобразователей № 1). Для фиксации положения модели используются вертикальные стержни.

На данный момент система укомплектована двумя парами пьезокерамических преобразователей P-волн (№ 1 и 2). В будущем она будет дополнена многокомпонентными преобразователями, аналогичными описанным в (Chichinina et al., 2009), что позволит регистрировать продольную и две поперечные волны в каждой точке измерения. Основа пары преобразователей № 1 -- пьезокерамические диски диаметром 10 мм и толщиной 3 мм, где в качестве демпфера используется эпоксидная композиция с окисью железа. Для использования в составе пары № 2 первоначально были изготовлены сильнозадемпфированные преобразователи на основе пьезокерамических дисков диаметром 1,5 мм и толщиной 1 мм с латунным демпфером. Однако значения их чувствительности оказалось недостаточно для достижения необходимого соотношения сигнал/шум, и они были заменены на преобразователи на основе пьезокерамических дисков диаметром 10 мм и толщиной 3 мм, снабженных коническими латунными концентратами с пятном контакта около 3 мм. Данные преобразователи обеспечивают достаточное соотношение сигнал/шум при удовлетворительной форме сигнала.

В состав установки включена аппаратура для управления перемещением преобразователей, генерации и регистрации сейсмического сигнала (рис. 3). Большая часть используемых устройств аналогична описанным в (Popovet al., 2015; Polovkovet al., 2015). Основу составляет персональный компьютер с установленным программным обеспечением, написанным в среде C++, и подключенной многофункциональной платой ввода-вывода ADLinkDAQ-2010. Цифровые входы и выходы платы используются для подачи сигнала на драйверы шаговых двигателей при перемещении пары преобразователей № 2, а также для контроля состояния оптических датчиков. К входам АЦП платы через усилитель заряда с коэффициентом усиления 25-1200 и частотным диапазоном от 20 до 3 МГц подключены ультразвуковые приемники. Для выбора одного из двух приемников используется блок переключателя; в будущем данное устройство будет оснащено цифровым управлением для автоматического выбора регистрирующего канала.

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

Рис. 3. Блок-схема измерительной установки

Управляющая программа реализует полностью автоматический цикл измерений с записью результатов в стандартном сейсмическом формате SEG-Y после масштабирования сигнала (коэффициент масштабирования 1: 20 000) и с заполнением всех необходимых заголовков трасс. Реализована возможность гибкого задания параметров съемки -- возможно варьирование диапазона, шага перемещений, разницы азимутов приемника -- источника и накопления сигнала при регистрации.

Минимальный дискрет перемещения преобразователей составляет 8,1'. Таким образом, диапазон азимутов подвижного источника 2 (рис. 1В) от 0 до 90° может быть разбит более чем на 650 отдельных измерений, что обеспечивает накопление необходимого статистического объема данных. Также возможно измерение в четырех квадрантах (диапазонах 0-90, 90-180, 180-270 и 270-360°), что позволяет оценить возможную неоднородность используемой модели.

При проведении измерений для каждой трассы пары преобразователей № 2 вычисляется реальный угол между направлением источник -- приемник и осью симметрии модели (вертикальной осью). При движении преобразователей в противофазе (с разницей азимутов 180°) используется следующая формула:

где ф -- значение истинного угла, r-- радиус окружности, описываемой преобразователями, A-- толщина модели, а и Я -- азимуты преобразователей.

Исходя из размеров модели и радиуса окружности, равного 87,5 мм (выбранного исходя из минимизации влияния боковых отражений), диапазон значений истинных углов составляет 15-90°.

Аппаратура, формирующая путь сигнала, показала полное отсутствие изменения формы регистрируемого сигнала в процессе эксперимента.

Для оценки точности изготовления модели была произведена серия измерений при отсутствии сдвига преобразователей пары № 2 (при их расположении друг напротив друга). Среднее значение отклонения времени первого вступления составило 0,16 % (0,79 мс при времени пикировки 486 мс после масштабирования записей). Также была определена скорость распространения продольных волн в пластинах вдоль напластования, равная 2375 ± 18 м/с.

Рис. 4. Результаты двух независимых измерений с идентичными параметрами и разница между ними

На рис. 4 приведен пример зарегистрированных данных оценки повторяемости измерений. Сейсмические трассы регистрировались при одинаковых настройках аппаратуры и значении сдавливающего усилия (7 МПа) с промежутком в 3 часа и снятием давления между измерениями. Вычитание результатов съемок показывает очень высокую степень соответствия сигналов.

2. Пример использования оборудования. Оценка параметров анизотропии пластинчатой модели

На рис. 5 приведены результаты эксперимента по оценке зависимости фазовой скорости Vpот угла ф. Угол ф измеряется относительно вертикальной оси z(она является осью симметрии среды VTI), то есть ф = 0° показывает направление, перпендикулярное плоскости слоя пластин в образце. Согласно работе (Thomsen, 2002), фазовую скорость Vpдля случая малой анизотропии в зависимости от угла ф можно оценить следующим образом:

где Vp₀-- вертикальная скорость под углом ф = 0; угол ф измеряется относительно вертикальной оси симметрии среды VTI; Ј и б -- параметры анизотропии (Thomsen, 1986):

Формулу (1) для скорости Vp^) можно записать в виде:

гдеx =sin²ф, a = V_p0(s-ф), b = V_p0ф иc = V_p0.

По формуле (2) и на основании полученных экспериментальных данных, используя метод наименьших квадратов, вычислены коэффициенты a, b и с: Vp(ф) = 183,63sin⁴ф-165,51 sin²ф + 2353,5. По следующим формулам высчитываются параметры анизотропии Ј и б:

s = (a + b) / с, ф = b / с,

что дает Ј = 0,01, б = -0,07 и фазовую скорость V_p0= 2353,5 м/с. Теоретическая зависимость Vp^) хорошо описывает экспериментальные данные (рис. 5). Достаточная плотность данных в диапазоне углов от 0 до 90° позволяет с высокой точностью оценить б и превосходит обычный способ оценки по единичному значению скорости при 45° Vp45 (Sarout et al., 2015).

Рис. 5.Скорость по отношению к углу между траекторией волны и вертикальной осью симметрии образца VTI. Приложенное вертикально давление составляет 7 МПа

3. Оценка погрешностей измерений

Для общей оценки погрешности измерений необходимо оценить погрешности определения всех параметров, используемых в расчетах. Для начала оценим погрешность позиционирования пары преобразователей № 2 на исследуемой модели. Абсолютная погрешность положения вала шагового двигателя составляет 5 % от шага 1,8°, то есть 5,4 угловой минуты, погрешность определения нулевого положения преобразователя при комнатной температуре -- 0,002 мм, то есть 0,09 угловой минуты. За счет люфтов прижимного устройства преобразователь может незначительно менять свое положение при наличии неровностей на модели, ориентировочная оценка этих изменений составляет 0,25 мм, то есть 9 угловых минут. Линейные размеры модели измерены штангенциркулем с точностью измерений 0,1 мм, точность определения радиусов движения преобразователей и соосности осей вращения составляет примерно 0,2 мм. На данный момент нами не учтено изменение размеров модели под нагрузкой, поскольку полагаем, что в этом пока нет необходимости. Итого абсолютная точность определения углового положения преобразователей составляет 15 угловых минут, или 0,25°, относительная точность определения расстояний -- около 0,2 %.

Оценить точность определения времени прихода волн несколько сложнее. Дискретизация данных происходит с шагом 0,5 мкс, точность определения амплитуд +/-0,1 % от входного диапазона. Регистрируемые задержки прихода сигнала составляют примерно 80 мкс, однако точность пикировки превышает интервал дискретизации, так как определение точки первого срыва сигнала происходит по огибающей, и может быть оценена в 1/10 дискрета. Таким образом, точность определения значений скоростей может быть обозначена как величина, меньшая 0,3 %.

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

Заключение

На слоистой модели из оргстекла изготовлена и апробирована установка, предназначенная для проведения измерений на синтетических моделях вертикально-анизотропных сред (VTI). Она позволяет производить измерения высокой точности и повторяемости в широком диапазоне углов между направлением распространения волны и вертикальной осью симметрии, а также значений вертикального сдавливающего усилия. Для используемой на данный момент модели размером 200x50x200 мм эти значения составляют 0, 15-90° и до 40 МПа соответственно.

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

Литература/References

1. Chichinina, T. I., 2017. Physical constraints on C13 for transversely isotropic shales and their applications. Geophysics, 82 (4), 105-118.

2. Chichinina, T. I., Obolentseva, I. R., Dugarov,G., 2015. Effective-medium anisotropic models of fractured rocks of TI symmetry. Analysis of constraints and limitations in linear slip model. In:85th Annual International Meeting, SEG, New Orleans, 421-426.

3. Chichinina, T. I., Obolentseva, I. R., Ronquillo-Jarillo, G., Sabinin, V. I., Gik, L. D., Bobrov, B. A., 2007. Attenuation anisotropy of P- and S-waves: Theory and laboratory experiment. Journal of Seismic Exploration 16, 235-264.

4. Chichinina, T. I., Obolentseva, I. R., Gik, L. D., Bobrov, B. A., Ronquillo-Jarillo, G., 2009. Attenuation anisotropy in the linear-slip model: Interpretation of physical modeling data. Geophysics 74 (5), 165-176.

5. Far, M., 2011. Seismic characterization of naturally fractured reservoirs, PhD. thesis. Houston.

6. Far, M., Figueiredo, J. J. S., Stewart, R. R., Castagna, J. P., Han, D.-H., Dyaur, N., 2014. Measurements of seismic anisotropy and fracture compliances in synthetic fractured media. Geophysical Journal International 197, 1845-1857.

7. Gik, L. D., Bobrov, B. A., 1996. Experimental laboratory study of anisotropy for fine-layered media. Geology and Geophysics 37, 97-110. (In Russian)

8. Hsu, C.-J., Schoenberg, M., 1993. Elastic waves through a simulated fractured medium. Geophysics 58, 964977. https://doi.org/10.1190/L1443487

9. Hudson, J. A., 1980. Overall properties of a cracked solid. Mathematical Proceedings of the Cambridge Philosophical Society 88, 371-384.

10. Polovkov, V. V., Popov, D. A., Musin, M. V., Dzemchuzhnikov, E. G., 2015. Physical modelling in exploration seismology: Reasonability, limitations and perspectives of the method. In: Geomodel-2015. Gelendzhik (In Russian)

11. Popov, D. A., Musin, M. V., Molodtsov, D. M., 2015. A hardware-software complex for physical seismic ultrasound modeling. Tekhnologii seismorazvedki 2, 113-118. (In Russian)

12. Sarout, J., Delle-Piane, C., Nadri, D., Esteban, L., Dewhurst, D. N.,2015. A robust experimental determination of Thomsen's Ф parameter. Geophysics 80 (1), 19-24. .1Schoenberg, M., 1980.

13. Elastic wave behavior across linear slip interfaces. Journal of Acoustical Society of America 68, 1516-1521. Schoenberg, M., 1983.

14. Reflection of elastic waves from periodically stratified media with interfacial slip.

15. Geophysical Prospecting 31, 265-292.

16. Schoenberg, M., Sayers, C. M., 1995. Seismic anisotropy of fractured rock. Geophysics 60, 204-211.

17. Thomsen, L., 1986. Weak elastic anisotropy. Geophysics 51, 1954-1966.

18. Thomsen, L., 2002. Understanding seismic anisotropy in exploration and exploitation. 2002 Distinguished Instructor Short Course, Distinguished Instructor Series 5. SEG, Tulsa.

19. Vernik, L., 2016. Seismic petrophysics in quantitative interpretation. Investigations in Geophysics 18. SEG, Tulsa.

20. Yan, F, Han, D.-H., Yao, Q.,2016. Physical constraints on c13 and Ф for transversely isotropic hydrocarbon source rocks. Geophysical Prospecting 64, 1524-1536.

Размещено на Allbest.ru