Русский
Русский
English
Статистика
Реклама

Геология

Анонс. Машинное обучение в геологии

17.06.2021 20:21:19 | Автор: admin
Завтра, 18 июня в 15:00 в наших соцсетях выступит Лейла Исмаилова, специалист машинного обучения в геологии и со-ведущая подкаста о геологах Про вулканы и людей

Лейла окончила геологический факультет МГУ им. М.В. Ломоносова. Поступила в аспирантуру Баварского Геологического Института в Германии. Во время обучения в аспирантуре опубликовала статьи в престижных научных журналах (Nature и Science Publishing group) и работала в разных лабораториях в Германии, Франции и США. С подробным списком публикаций можно ознакомиться по ссылке.

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




План выступления


Задавайте вопросы в комментариях и Лейла ответит на них во время прямого эфира.

  • Наука в геологии
  • Как стать ученым в России
  • Наука в России и наука зарубежом
  • Машинное обучение в нефтянке

До встречи в эфире!




Подробнее..

Как из днк-секвенатора сделать георадар и увидеть все под землей

13.01.2021 00:05:12 | Автор: admin
Какое-то время назад я опубликовал статью про то, как я занимался созданием ДНК-секвенатора (прибора для чтения ДНК) путём реверс-инжиниринга микросхемы. Но со временем, пока я это делал (а делал я это очень долго) этот секвенатор устарел, вышел новый, и мой проект отчасти потерял смысл, поскольку себестоимость секвенирования на новом секвенаторе Oxford Nanopore оказалась ниже, чем та стоимость, к которой стремился я. Проект заглох. С другой стороны, произошли события, которые сместили мой основной интерес в область такой вещи, как георадар. Точнее вернули его туда.

Лет 10 назад я наткнулся на ряд исследований, где применялось устройство, способное видеть, что находится под землёй георадар. Неделю не мог оторваться, просмотрел огромное количество видеоматериалов, запомнилось название того прибора, который использовали исследователи георадар ЛОЗА. Ниже пример исследования Саксайуамана (Перу), на котором найдены подземные склепы, котлованы, лежащая на глубине плита, а также можно предположить наличие засыпанного когда-то дна чашеобразной формы.


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

Из информации в интернете я составил представление о том, как устройство работает. Есть Передатчик и Приёмник. Передатчик посылает очень мощный, под несколько десятков киловольт, наносекундный электромагнитный импульс в глубь земли.

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

image

imageРазличные источники подчёркивали, что в георадарах используются так называемые резистивные (resist сопротивляться) антенны. Это специальные антенны, у которых нет звона они не резонируют. Обычная антенна (рис. 1) резонирует (начинает вынужденно колебаться) на собственной рабочей частоте. Её собственные колебания не дают ей качественно воспринять полезные сигналы, приходящие в это время. Импульсы, посылаемые при помощи резистивной антенны, отличаются от классического радиоимпульса отсутствием несущей частоты получаются ассиметричные непродолжительные импульсы (рис. 2). Если импульс отражается от среды с большей диэлектрической проницаемостью, то он инвертируется (рис. 3). Я начал соображать, как сделать такие антенны.

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

Представьте себе, что мощные 10-киловольтные импульсы длительностью в 1 наносекунду следуют один за другим, многократно отражаясь. Основная задача в том, чтобы всю приходящую на Приёмник информацию быстро оцифровать (перевести из аналогового сигнала в понятные электронике, возможные для дальнейшей обработки символы ноли и единицы), проанализировать и записать. А это гигабайты нолей и единиц в секунду непрерывным потоком.
Реализую-ка я этот девайс на компараторах для начала. Сущие копейки и многократное удешевление решил я. Компаратор это простейшая электронная схема, единственная задача которой, сравнить два входящих аналоговых сигнала. Она выдает на выходе 0 либо 1, в зависимости от того, какое из двух входящих напряжений больше. То есть компаратор это 1-битный аналогово-цифровой преобразователь (АЦП), позволяющий записать двоичную форму сигнала. Однако, чтобы восстановить всю амплитуду одного возвращённого из земли импульса, нужны сотни таких сравнений с разными порогами компаратора. Моё сверхбюджетное решение не могло похвастаться скоростью и эффективностью.

В общем, я медленно обдумывал подходы к реализации прибора, пока не познакомился с Владимиром Зубовым и не занялся всерьёз реверс-инжинирингом ДНК-секвенатора, о чём впоследствии опубликовал на habr.com статью. Всё это время я продолжал с интересом следить за георадарными исследованиями, но времени на всё не хватало. Каково же было моё удивление, когда цепь случайных событий привела меня в ту же точку с совершенно с другой стороны.

Три-четыре года назад я переехал жить под Троицк. Как раз тогда опубликовал статью про секвенатор, сидел разбирал входящие. Сам Jonathan Rothberg предложил мне start-up, но что-то не заладилось, было много других отнимающих время входящих в общем был полностью погружён в последствия публикации. Неподалеку от нашего дома на поле были огромные антенны на мачтах (ну очень здоровые), частенько ходили туда с квадрокоптером. Я заинтересовался ими, стал гуглить большие антенны, Троицк, наткнулся на сайт rk3b.ru c интригующим названием Школьный Центр космической связи, позвонил туда и напросился в гости сам радиолюбитель.

imageТак я познакомился с Александром Николаевичем Зайцевым, заслуженнейшим товарищем, который этот Центр возглавлял. Оказалось, что он много лет занимался исследованиями магнитосферы Земли в ИЗМИРАН и заодно знал весь научный бомонд г. Троицка. Разговорились, я упомянул, в том числе, про свой интерес к георадарам. Александр Николаевич в свою очередь рассказал мне про антенны на поле. Они оказались чисто любительскими: известный радиолюбитель В. Н. Комаров, будучи при этом успешным предпринимателем, собрал команду энтузиастов КВ связи, на свои деньги создал суперцентр и такие антенны, что они стали самыми лучшими в эфире. Подтверждение тому первые места в соревнованиях на первенство мира по КВ-связи.

А через некоторое время А. Зайцев познакомил меня с людьми из ООО Компания ВНИИСМИ, которые (сюрприз!:) и занимались георадарами ЛОЗА. Я встретился с П. Морозовым и А. Беркутом, её возглавлявшими. Они с радостью меня приняли, и мы во многом сошлись. Представьте себе моё изумление, когда они рассказали, что их устройство, георадар ЛОЗА, продаваемый по цене 25 000$ в базовой комплектации, работает на компараторах. И они как раз очень давно мечтали сделать устройство на аналого-цифровом преобразователе (АЦП), который позволил бы прибору работать быстро. И что попытки разные были, но никак успехом не увенчались. В общем, высказали полнейшую свою заинтересованность в разработке.

Впоследствии я получил полную информацию о том, как работает георадар ЛОЗА. Там и правда не было ничего кроме компараторов и нескольких микросхем, лучшая из которых, впрочем, была самой передовой в семействе FPGA лет 20 назад. Становилась понятна главная жалоба ВНИИСМИ на свой прибор: надоело пальцем жать на кнопку. Ведь что такое работа на компараторах: для одного измерения (измерения в одной точке) нужно последовательно послать в землю 128 импульсов, что занимает от секунд до 2 минут, в зависимости от передатчика. То есть прибор работает по принципу: нажал на кнопку, постоял-подождал, передвинул прибор дальше по рулетке на 10 см. (действительно расстилается по земле рулетка, чтобы двигаться с шагом в 10 см.) и так далее. Представляете себе скорость такой работы? Кстати, блок управления либо компьютер подключаются к Передатчику проводом, поэтому для съёмки необходимы минимум 2 человека: тот, кто перемещает георадар, и оператор, идущий за ним на проводе.

imageТакже мне представилась возможность изучить несколько разных георадаров серийных устройств и прототипов разных производителей и разработчиков.

Например, вот эта большая тяжёлая коробка весит килограммов 5 или около того; внутри там стоит АЦП на 1800 Мегасемплов в секунду, каждая плата в отдельном корпусе и с отдельным экраном в итоге много лишнего веса и большие габариты, что сделало прототип невозможным к использованию в принципе; он не работал.

Другой прототип был поудачнее, но проблема была в том, что его надо было очень долго настраивать; иногда он не включался, сбоил. Вообще, одна из проблем георадара это необходимость работать в очень большом динамическом диапазоне. Импульс, который отправляется под землю, по мере прохождения сквозь неё очень быстро затухает. Чтобы увидеть и первый отклик этого импульса, и отклик, пришедший глубоко из-под земли, нужны усилители и АЦП, способные принимать и различать как очень сильные сигналы, так и очень-очень слабые. Вот эта амплитуда и называется динамическим диапазоном. В данном прототипе для этого использовался двухканальный усилитель один канал усилителя работал с сильными сигналами, а другой со слабыми. То есть было 2 канала оцифровки, которые, к слову, никак не удавалось свести воедино. Регулируя усиление, прибор можно было настроить на какую-то глубину и в ней он что-то видел. Но не выше, не ниже заданного уровня он не видел, пока его не перенастроишь на иные параметры. Кроме того, у прототипа была очень большая чувствительность к помехам.



Итак, я начал проектировать собственное высокоскоростное, мощное устройство. Работа предстояла громадная, но были и хорошие новости: часть работы была уже проделана за основу я взял ту электронику, которую (та-дам!..:) я разработал для ДНК-секвенатора.

Работа над георадаром и комплексом программного обеспечения заняла около года. Я добился не только полной автоматизации, высоких скоростей, возможности управлять георадаром с любого устройства, мне удалось создать, по сути, универсальное ядро под любой высокотехнологичный проект: система на модуле (Zturn) система на кристале Zynq связка работы FPGA+CPU+Linux+Django+WebSockets+Javascript Впрочем, оставлю интригу. Чтобы не утомлять читателя, выношу техническую часть в отдельную статью.

Разработанный мною прибор представляет из себя 2 устройства Передатчик и Приёмник, у каждого своя антенна. Сейчас размеры устройств 22,2 х 14,6 х 5,5 см., в следующей партии планируется ещё уменьшить габариты. Для исследований на небольшой глубине применяются стандартные метровые антенны. В зависимости от длины и мощности антенн глубина зондирования на благоприятных почвах может достигать нескольких сотен метров. Для перемещения георадара может быть использована ручная переноска, платформа на колёсиках, складные/гибкие антенны глубинного зондирования с посадочными местами под Приёмник и Передатчик; под спец. задачи могут быть созданы иные средства перемещения (например, надувные для исследования подводных объектов).

Платформа с георадаром может быть прикреплена к транспортному средству, георадар способен производить съёмку в автоматическом режиме на скорости до 40 км/ч (то, что мы протестировали, вероятно и больше). Регистрация перемещения прибора осуществляется при помощи датчика колеса и GPS. Все процессы автоматизированы, прибор лёгок в управлении, не требует дополнительного оператора, со съёмкой справится и один человек. Приоритеты: мощность, компактность, лёгкость, пассивная охлаждаемость, возможность комплектации с различными средствами перемещения и антеннами разной мощности. Эти две небольшие коробочки могут стать сердцем георадарного комплекса любой мощности. Дальше скучное, спецификация:

Частотный диапазон (МГц) 1-300
Частота следования зондирующих импульсов до 1000 импульсов в секунду (при максимальной частоте импульсов при движении со скоростью 36 км/ч интервал измерений составит 1 мм)
Тип импульса: широкополосный видеоимпульс, без несущей частоты
Мощность в импульсе 5 МВт
Длительность импульса 1 нсек
Приёмник: частота дискретизации 1GSPS (1нсек)
Режим регистрации АЦП 16 бит
Усилитель логарифмический
Динамический диапазон изменяемый до 120 Дб
Временные развертки (количество отсчетов в одном сигнале зондирования) 16000
Автоматическая синхронизация с передатчиком
Встроенная память 64 Гбайт
Дополнительные датчики: барометр, акселерометр (измерение угла наклона прибора), GPS, энкодер колеса
Управление кнопками на приборе или через wifi с любого устройства. Есть web интерфейс, программа для PC
Время работы от встроенной батареи 4 часа, возможно присоединение внешнего аккумуляторного блока (+12 часов)
Встроенный генератор тестового сигнала
Габаритные размеры 222 х 146 х 55 mm (х 2)


image


Стандартные антенны имеют центральную частоту 100Mhz (1,5м), 200MHz (1м). Существуют также низкочастотные антенны для глубокого зондирования 10MHz (10м), 25Mhz (6м), 50Mhz (3м) (в скобках размер антенн, а не глубина зондирования). Увеличивая длину антенн увеличиваем глубину зондирования. Важно отметить, что глубина зондирования зависит также от мощности передатчика и от плотности/свойств грунтов. Таким образом, всегда можно увеличить антенны, поэтому глубина зондирования может достигать 500 м.

Георадар готов к работе сразу после нажатия кнопки ON георадар начнёт вести съёмку и записывать всё во внутреннюю память. Он может работать абсолютно автономно, без подключения к нему оператора.

Устройство раздаёт wi-fi, подключиться к нему можно с любого компьютера, планшета, телефона через web браузер. Вы попадаете в программу, позволяющую управлять устройством, менять настройки, видеть текущий трек. В верхней строке интерфейса различные статусы, позволяющие понять, что происходит в данный момент: напряжение аккумулятора, температура на процессоре, время, время работы, GPS данные. Ниже вкладки и кнопки настройки триггеров, offsetов, логарифмической шкалы, зума, режимы работы.

Прибор может производить одно измерение по нажатии кнопки, либо по времени, например, каждые 0,3 секунды, либо измерения могут быть синхронизированы с работой колеса, например на каждую четверть оборота колеса происходит срабатывание, то есть Передатчик посылает импульс, а Приёмник его принимает и записывает. Последний режим очень удобен, так как позволяет привязать платформу с георадаром к машине, например, и производить равномерные измерения независимо от скорости движения.

Первый запуск был произведён совместно с представителями ВНИИСМИ на строящейся асфальтной дороге на юго-западе Москвы. Вот так выглядела инсталляция моего прибора: 2 антенны на платформе максимально близко к земле, белая коробочка это Передатчик, чёрная Приёмник.

imageЗаранее было выбрано место, где в разных точках под землей лежали коммуникации. По плану было определить такие точки сначала прибором Лоза, а затем тоже самое место посмотреть моим прибором. Пока ВНИИСМИ свой прибор настраивали, я с пол оборота завёлся и проехался, и всё увиделось. Вот первая картинка, полученная с моего прибора. Мы подошли к рабочим, которые в это время занимались там строительством дороги и сверили результаты наших съёмок с картами коммуникаций, которые были у них всё совпало, там действительно были пластиковые трубы 200-го диаметра на глубине 2 м. Это был успех.

image

На картинке ниже кусок той же съёмки строящейся дороги. Справа вы видите извивающийся столбик это осциллограмма импульса. Маленькая пика вверху это тот импульс, который сгенерировал Передатчик, а всё что ниже, это то, что приходит на Приёмник со всех сторон, но в основном из под земли. Большая цветная картина рядом с осциллограммой импульса это набор таких столбцов, то есть каждый пиксельный столб это 1 импульс, только здесь амплитуда переведена в цвет. Красная полоса в верхней части цветной картинки это уровень земли. Дальше импульс уходит в землю, и ниже мы видим всё, что отразившись прилетело на Приёмник. Небольшие сине-красные бугорки чуть ниже второй красной полосы и есть пересечение диаграмм направленности антенн с точечными объектами в виде труб. То есть, в тех местах, где по ходу движения прибора мы пересекаем трубы, трубы дают вот такой вот радио-образ, причём в некоторых местах видно две трубы, если увеличить.

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

imageВот на этой картинке ручей, закопанный при строительстве дороги; видно русло на глубине порядка 2 метров. Ниже на этом участке тоже были большие параболы отражение от столбов освещения, вдоль дороги.

У меня рядом с домом есть река Незнайка с очень интересной геологией, где всё как в учебнике: есть старое русло, четвертичные отложения, аллювий, речные террасы. Там тоже мы ходили с ребятами из ВНИИСМИ испытывали, сравнивали работу разных георадаров. Здесь в конце видно старое русло реки, а большие зелёные параболы это отражение забора соседнего посёлка на металлических столбах.
imageКак видите, только человек уже знакомый с принципами интерпретации таких картинок, сходу определит, что это русло реки, а это отражение от забора. Таким образом, сама форма вывода данных существенно ограничивает круг людей, способных легко начать работать с георадарами. И здесь, по-моему, очень интересная задача стоит перевода этих данных в понимаемую любым пользователем картинку. Выполнена эта задача может быть с применением искусственного интеллекта (нейросети), который можно обучить маркировать, обозначать реальные контуры объектов и предполагать их назначение и свойства. Также решить проблему нечитабельности данных помогут объёмные построения, но об этом дальше.

На картинке ниже очень интересный рельеф (съёмка произведена 22.07.2020 в Ивановской области, д. Калинкино): видна какая-то продолговатая структура скорее всего это слой песка либо ещё какого-то менее плотного, чем превалирующие в наших краях суглинки, грунта. У меня есть мысль полностью автоматизировать идентификацию слоёв. Как мы определили выше, волна отражается от тех мест, где происходит изменение диэлектрической проницаемости среды. Мы можем померять скорость прохождения импульса на том или ином ярусе этой картинки и предположить состав среды и/или степень её влажности.



В конце лета мне удалось съездить в археологическую экспедицию под Рязанью близ села Терехово Шиловского района на место стоянки/поселения пятого века. Там на схождении двух рек есть холм, поросший деревьями остались многочисленные свидетельства, что в пятом веке там было Городище. Стали ходить там с георадаром на первом рисунке наш трек по GPS. В центре очень густо росли деревья, мы смогли пройти по периметру, дальше деревьев было меньше и мы всё довольно плотно исходили. У меня была возможность воспользоваться программой, которая анонсировалась, как программа позволяющая строить 3D модель из линейных проходок. Такого результата ждали от меня археологи построим и увидим наглядно что где зарыто. Лучше всего подобные построения делать из хорошо (плотно) исхоженных участков. Я загрузил в программу область с самыми плотными проходами и вошёл в режим 3D программа выстроила поле из цветовых столбиков, которые можно было вытягивать вниз и вверх, меняя настройки, больше никаких возможностей я не увидел проконсультировавшись выяснил, что программа строит лишь псевдо 3D.

image

Представим, что у нас есть некие столбцы, которые нам известны, где георадар сверху проехал, и мы знаем, что вертикально вниз в этой точке находится. Следующая задача, это достроить недостающие столбцы аппроксимацией, что и сделала программа. Но сделала она это банально решёткой по 2 осям. В результате такого построения возникают различные артефакты крестообразные, из горизонтальных и вертикальных линий. Если мы берём какой-то срез 3D плоскости, мы имеем цветную картинку из среза столбиков, которые программа просто вверх либо вниз вытягивает в зависимости от амплитуды сигнала.

image

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

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

Где может быть использован георадар?
Самые верхние под поверхностью земли слои называют культурным слоем это археология, прекрасно видны отдельные объекты, фундаменты строений, видны даже копаные и засыпанные когда-то ямы. Так, есть очень интересная задача связанная, например, с пустыней Сахара. Сахара расширяется, и многие археологи сходятся во мнении, что под её песками погребены города и строения. Во времена Египта там была саванна, текли реки, жили животные. Георадар может просканировать всё до грунта, и найти там строения и поселения. Это собственно единственный прибор, который способен это сделать.

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

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

image

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

В декабре 2020 года базовая версия прибора успешно прошла испытания и была сертифицирована. Георадар получил Название GEORA.

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

К широкому выпуску планируются 2 версии: для частного использования и для коммерческого. Есть возможность кастомизировать оборудование и софт под любые задачи.

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

Сейчас готовится к публикации вторая часть статьи с подробным описанием технической части. На сегодня всё, если интересно пишите, почта sokolov.labs@gmail.com. Всем пока, спасибо за внимание!
Подробнее..

Не только детекторы. Экскурс в прикладную физику нейтрино

17.01.2021 22:08:31 | Автор: admin

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

Немного истории

Нейтрино это совокупное название группы легчайших элементарных частиц, относящихся к фермионам. Существование нейтрино было предсказано Вольфгангом Паули в 1930 году, а экспериментально подтверждено в 1956 году Клайдом Коуэном и Фредериком Рейнесом. При этом Паули лишь неформально, в виде чистой гипотезы, предположил, что имеется возможность того, что в ядрах существуют электрически нейтральные частицы, которые я буду называть нейтронами и которые обладают спином . Масса нейтрона по порядку величины должна быть сравнимой с массой электрона и во всяком случае не более 0,01 массы протона. Таким образом он пытался объяснить наблюдаемую природу бета-распада. Он назвал такую неоткрытую частицу нейтроном. Только через два года, в 1932, Джеймс Чедвик открыл в атоме крупную элементарную частицу, сравнимую по массе с протоном, и назвал ее нейтрон, а неуловимый фермион Паули впоследствии удостоился названия нейтрино (нейтрончик) с легкой руки Энрико Ферми.

С тех самых пор нейтрино окружены ореолом загадочности в силу своих поразительных свойств. Они всерьез и надолго обосновались в научной фантастике так, Кельвин, главный герой Соляриса, предполагает, что именно из нейтрино разумный океан формировал своих фантомов, в том числе, фантом Хари, возлюбленной Кельвина. Вкратце напомню основные уникальные и парадоксальные аспекты нейтрино:

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

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

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

С этими и другими темами мне довелось познакомиться подробно и достаточно давно. В 2014-2015 году я работал переводчиком в издательстве Альпина Нон-Фикшн, и в этот период, в первые несколько месяцев 2015 года, перевел книгу Рэя Джаявардханы Охотники за нейтрино. На тот момент мне казалось, что эта тема слишком академична и сложна для широкой читательской аудитории, но книга оказалась настолько интересной, а я так надолго увлекся этой темой, когда она еще не была мейнстримом, что возвращаюсь к ней до сих пор. Разброс рассматриваемых тем, связанных с нейтрино, впрочем, меняется слабо. Пишут о новых детекторах, все более точных измерениях массы и скорости нейтрино, о космологической ценности этих исследований. Не буду здесь углубляться в эти аспекты, поскольку они хорошо и подробно рассмотрены на русском языке (и при этом богато иллюстрированы). Но позволю себе привести список статей с Хабра, которые вышли уже после книги Джаявардханы и кажутся мне наиболее интересными:

  1. BAIKAL-GVD. Охотники за нейтрино (2.02.2020). Пост автора @DNLP о кластере глубоководных детекторов нейтрино, устанавливаемых в озере Байкал. Отличный материал с большим количеством видео и техническими деталями. На самом деле, идея использовать Байкал в качестве естественного водного резервуара гениальна, а географическое расположение детектора позволяет в большом количестве отлавливать нейтрино, прилетающие к нам со стороны Южного полюса. Аналогичный, совсем свежий материал недавно вышел на сайте Медузы.

  2. Вещество и антивещество: что это такое, в чем разница и при чем тут нейтрино (24.09.2019). Автор @Bars21. Подробный разбор парадокса #3 из вышеприведенного списка.

  3. Поймай меня, если сможешь: радиоволны, каскад частиц и лед для поимки нейтрино (11.03.2020) Пост автора @Dmytro_Kikot, дающий представление о подледном лове нейтрино.

  4. Как неудавшийся ядерный эксперимент случайно породил нейтринную астрономию (15.07.2018) перевод с английского от редактора Хабра @SLY_G. Действительно, в тематических источниках зачастую умалчивается, что знаменитые нейтринные детекторы в шахтах Камиока в Японии (Super-Kamiokande, Hyper-Kamiokande) были оборудованы на месте совершенно другого физического мега-эксперимента, не увенчавшегося успехом, и переориентация этого комплекса на изучение нейтрино во многом счастливая случайность.

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

Нейтрино и телекоммуникация

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

В 2010 году было высказано предположение, что нейтрино могут применяться для однонаправленной связи с субмаринами, постоянно находящимися в подводном положении. Пучок нейтрино можно было бы направить в расположенный в условленной точке океана детектор. При попадании нейтрино в такой детектор рождались бы мюоны, испускающие излучение Черенкова, а подводная лодка могла бы считывать переданную информацию при помощи детекторов, проходя мимо при патрулировании зоны. В источнике такая точка именуется почтовый ящик; указано, что в силу почти полного отсутствия помех на такой глубине, ящиков в заданной зоне могло бы быть несколько, а скорость передачи информации в мюонное хранилище (muon storage) могла бы составлять более 100 бит/с.

Предыдущий пример подсказывает, что наибольшую проблему в данном случае составляет вычленение сигнала из шума, то есть, необходимость создания детектора, который различал бы информативную последовательность нейтрино на фоне нейтрино естественного происхождения. Кроме того, в примере с подводными лодками не удается уйти от пропускания пучков нейтрино через толщу воды, которой в данном случае служит сам океан. Но сохранились сведения о работах, предполагающих нейтринную коммуникацию без использования водного резервуара. В 2016 году исследовательская лаборатория Pirelli в Милане опубликовала материал о том, что с начала 2000-х ведутся работы по синтезу кристаллов кремния или кварца, используемых в качестве детекторов нейтрино. Кристалл охлаждается почти до абсолютного нуля, и при попадании пучка нейтрино температура этого кристалла немного возрастает, что можно зафиксировать при помощи приборов. Такие детекторы были бы не только гораздо компактнее водных резервуаров, но и обеспечивали бы когерентное рассеяние полученных нейтрино, и их осцилляции затем можно было бы считывать и расшифровывать. Правда, пока ни о каких практических результатах этих исследований не сообщается

Широкую известность получил эксперимент, поставленный в 2012 году в Национальной ускорительной лаборатории им. Энрико Ферми (Фермилабе) в Чикаго в ходе упоминаемого эксперимента ученым действительно удалось передать в пучке нейтрино информацию, а именно закодировать слово neutrino. Информация была передана с мощнейшего современного генератора нейтрино NuMI на расстояние более километра и зафиксирована детектором MINERvA. Вот как процесс проиллюстрирован на рис. 1 к упоминаемой научной статье:

Несомненно, перед нами лишь proof-of-concept (доказательство осуществимости), показанное почти 10 лет назад, но развитие нейтринной телекоммуникации ограничено техническими, а не фундаментальными сложностями. Таким образом, подобные технологии вполне могут стать реальностью.

Нейтрино и контроль над использованием ядерных реакторов

Одним из самых распространенных источников нейтрино на Земле являются ядерные реакторы. Предпринимаются попытки устанавливать уловители нейтрино (точнее антинейтрино) вблизи от ядерного реактора, на расстоянии в пределах 10 метров, для изучения свойств этих частиц. О таком эксперименте подробно рассказано в интервью д.ф.-м.н. Дмитрия Наумова, зам. директора по научной работе Лаборатории ядерных проблем Объединенного института ядерных исследований (ОИЯИ) в Дубне. При этом Наумов отмечает, что по свойствам улавливаемых антинейтрино можно достоверно определить, не используется ли реактор для производства оружейного плутония-239. Любая компания, занимающаяся эксплуатацией ядерного реактора, обязана предоставлять данные о том, сколько плутония получает в ходе работы, и нейтрино позволяют проверить эти данные. Впрочем, существовали гораздо более фантастичные проекты, связанные с нейтринным шпионажем. В 2010 году группа французских физиков во главе с Тьерри Ласьерром опубликовала статью о том, как, регистрируя потоки нейтрино, можно обнаруживать незадекларированные ядерные реакторы. Более того, авторы предлагали устанавливать нейтринные детекторы на кораблях и следить за развитием ядерных программ Ирана и Северной Кореи, заходя в прибрежные воды этих государств.

Нейтрино и геологические изыскания

Эта прикладная возможность отчасти смыкается с предыдущей предполагается, что детекторы нейтрино полезны при поиске месторождений урана и тория. Но гораздо более интересные возможности открываются при измерении свойств солнечных нейтрино, которые, как было указано выше, свободно пронизывают Землю насквозь. Изучая изменение осцилляций нейтрино при их прохождении через толщу пород можно было бы проводить томографию литосферы, находить в ней полости, анализировать плотность веществ, заполняющих эти полости. Такая технология открыла бы путь к обнаружению глубоких месторождений нефти. В частности, подобные идеи исследованы в статье перуанских ученых, опубликованной в 2015 году. Предполагается, что для такой цели могли бы использоваться не только солнечные нейтрино, но и направленные пучки нейтрино, сгенерированные искусственно. В таком случае геологоразведочные работы можно было бы существенно ускорить, затрачивая на поиск месторождений не годы, а месяцы. Впрочем, на момент публикации статьи еще не существовало столь мощных генераторов нейтрино, которые позволили бы воплотить эту технологию.

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

Заключение

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

Подробнее..

Перевод Как переделать планету атом за атомом

22.04.2021 10:12:30 | Автор: admin

Когда планета не планета? Когда идут гелиевые дожди? Как вода одновременно может быть в твердом и жидком состоянии? Чтобы ответить на эти вопросы, ученые берут вещества, из которых обычно состоят планеты, подвергают их экстремальному давлению и смотрят, что получится.

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

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

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

Это уловитель в Национальном комплексе лазерных термоядерных реакций, один из важнейших приборов для подобных экспериментов. Иллюстрация: Damien Jemison, Ливерморская национальная лаборатория им. ЛоуренсаЭто уловитель в Национальном комплексе лазерных термоядерных реакций, один из важнейших приборов для подобных экспериментов. Иллюстрация: Damien Jemison, Ливерморская национальная лаборатория им. Лоуренса

Мы стали размышлять, как реально можно создать в лаборатории такие условия, которые существуют в недрах планеты, сказалаАрианна Глисон, специалист по экспериментальной физике минералов в Национальной ускорительной лаборатории SLAC и в Стэнфордском университете в Менло-Парк, штат Калифорния. - берем самые бесхитростные минералы кварц, полевой шпат и задаемся вопросом: как они выглядят, каковы их свойства при экстремальных условиях?

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

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

Алмазы, пушки и лазеры

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

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

Научное сообщество уже почти 50 лет работает с давлениями порядка сотни гигапаскалей, то есть, 1 миллиона атмосфер, говорит Раймонд Джинлоз, ученый-планетолог из Калифорнийского университета в Беркли (давление на поверхности Земли составляет 1 атмосферу). сотня гигапаскалей важная отметка в нашей дисциплине, поскольку примерно соответствует давлению на границе ядра и мантии Земли. Вцентре ядра Землидавление примерно втрое выше, и такое давление вполне достижимо в новой, уменьшенной ячейке с алмазными наковальнями, где та же сила, что и в более ранней модели прикладывается к более мелкому образцу и, следовательно, давление удается увеличить.

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

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

Уикс, занимающаяся планетологией в Университете Джона Хопкинса, Балтимор, штат Мэриленд, ставит эксперименты по сжатию образцов лазером, изучая таким образом, как атомы и молекулы взаимодействуют глубоко внутри планет. В последние 20 лет на переднем крае этой дисциплины развивается сжатие с использованием оптических лазеров высоких энергий, как например в Национальной Ускорительной Лаборатории SLAC.

Фокусируем [лазер] на образце, и он превращает поверхность образца в плазму, а эта плазма расширяется и направляет равновеликую противоположную волну давления на образец, объясняет Уикс. Все это происходит за несколько миллиардных долей секунды.

В Z-машине для генерации высоких температур, высоких давлений и высокоэнергетических рентгеновских лучей используются электрические токи и магнитные поля. При помощи этой установки ученые изучают уравнения состояния тех веществ, из которых состоят планеты. Иллюстрация: Randy Montoya/Sandia National LaboratoriesВ Z-машине для генерации высоких температур, высоких давлений и высокоэнергетических рентгеновских лучей используются электрические токи и магнитные поля. При помощи этой установки ученые изучают уравнения состояния тех веществ, из которых состоят планеты. Иллюстрация: Randy Montoya/Sandia National Laboratories

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

Гелиевые дожди подсвечивают Сатурн

На Земле дождь состоит из жидкой воды, а на Сатурне из жидкого гелия. Это известно благодаря экспериментам с использованием лазера; такие эксперименты поставлены вНациональном комплексе лазерных термоядерных реакцийв Национальной Ливерморской Лаборатории им. Лоуренса. Они помогли определить, когда водород и гелий смешиваются друг с другом, а когда разделяются; это свойство называется смешиваемостью.

Водород самый распространенный элемент во Вселенной, он так или иначе присутствует в составе любого планетарного тела, как в чистом виде, так и в соединениях, например, в воде или метане, говоритТакуо Окучи, адъюнкт-профессор в Институте Планетарного Вещества при университете Окаямы в Японии. Он поясняет, что химическое состояние водорода может очень, очень отличаться в зависимости от условий окружающей среды, [в сущности] показателей температуры и давления.

Такие изменения на уровне атомов отражаются в планетарном масштабе. По словам Окучи, при давлениях, существующих в недрах Юпитера и Сатурна, водород становится металлическим, то есть, атомы водорода располагаются столь плотно, что их электроны перекрываются друг с другом. Жидкий металлический водород поддерживает магнитное поле внутри этих газовых гигантов (Окучи уточняет, что внутри Урана и Нептуна в металлическое состояние переходит вода.)

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

Джелноз сравнивает такую смесь с маслом и водой.

Сатурн примерно на 50% ярче, чем должен быть в своем возрасте и, как объясняет Стюарт, дело именно в гелиевых дождях: они идут на Сатурне, а на Юпитере их нет. Такаяидеясуществует уже некоторое время, но лишь сравнительно недавно появилась возможность создать такие условияв лаборатории.

На Нептуне небо в алмазах

На ледяных гигантах, Уране и Нептуне, выше доля метана (CH4), воды (H2O) и аммиака (NH3) чем на газовых гигантах, и эксперименты с динамическим сжатием показали, что там дожди еще более странные. Команда во главе сДомиником Краусомисследовала, что бывает с чистыми углеводородами, когда они попадают в условия, которые могут существовать в недрах такой планеты конкретно, Нептуна.

Мы зафиксировали формирование наноалмазов, говорит Краус. Давление, создаваемое при наносекундном лазерном сжатии, разрывает молекулярные связи между углеродом и водородом, в результате чего углерод спекается в нанометровые алмазы. Это открытие подтвердило давние теоретические выкладки.

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

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

Брать за основу результаты наносекундных экспериментов и экстраполировать их до планетарных масштабов это порой зыбкое дело. Определенно, какие-то вещи, происходящие за наносекунды, будут происходить и на протяжении миллионов лет, говорит Краус. Вопрос в том, что еще такого может произойти за это время, чего мы не замечаем в доступных нам микромасштабах?

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

Твердожидкая вода осложняет изучение газовых гигантов

Если мы сможем лучше понять, что происходит с веществом внутри Урана и Нептуна, это также поможет нам в изучении самого распространенного типаэкзопланет.

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

Эксперименты, опубликованные в 2018 году, что льдистость ледяных гигантов на самом деле, более сложное явление, чем считалось ранее. Мы обнаружили, что вода переходит в такое необычное суперионное состояние лишь при высоком давлении и при таких температурах, которые характерны для глубин Урана и Нептуна, говоритМариус Милло, физик из Национальной лаборатории им. Лоуренса в Ливерморе. Суперионный лед это новое состояние вещества.

Милло руководил исследовательской группой, открывшей это ранее не известное агрегатное состояние вещества. Сначала они попробовали использовать ячейку с алмазными наковальнями, а затем работали в лаборатории Omega Laser FacilityРочестерского университета, где смогли искусственно кристаллизовать воду в таком новом состоянии.

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

При давлении, устанавливающемся в мантии ледяного гиганта (около 200 миллионов атмосфер), суперионный лед тает при температуре около 4 700 C, гораздо выше, чем температура окружающей среды. В новом исследовании команда подтвердила свойства этого нового состояния льда. Вполне возможно, что суперионный лед не тает даже в самых глубинах Нептуна и Урана, считает Милло, поэтому ледяные гиганты могут оказаться вполне твердыми.

Атомытекучего водородаобладают электрическим зарядом, поэтому воздействуют с магнитным полем планеты, а, возможно, и влияют на него. Более того, структура и перенос энергии внутри такой планеты могут влиять и на другие наблюдаемые там явления, например, на погоду.

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

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

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

Мы как будто разобрали Землю и собрали ее заново

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

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

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

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

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

Многие планеты внутри рыхлые

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

Наши представления о слоистых планетах могут быть совершенно неверны. Планеты-суперземлимогут выглядеть подобно Земле железное ядро, дальше скалистые породы и атмосфера, более плотная, чем у нас или весьма отличаться от Земли в том, что камень и железо в их недрах могут смешиваться друг с другом в виде расплава, и четко выраженного металлического ядра у них не будет, говорит Стюарт.

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

Наше представление о слоистых планетах может оказаться полностью неверным. Мы пока не можем измерить эти величины экспериментально, но в течение ближайших 10 лет сможем, считает Уикс.

Импакты и хронология Солнечной системы

Примерно 4 миллиарда лет назад во внутренней части Солнечной системы шла настоящая каменная бомбардировка. Радиоизотопная датировка минералов, образовавшихся при таких столкновениях, позволяет оценить, когда это происходило, а следы ударных контактов в зернах минералов помогают понять силу удара. Циркониевый век или бадделеитовый век очень важны для нас при периодизации всей истории этой бомбардировки, - указывает Аи-Чэн Чжан, профессор минералогии в Нанкинском университете в Китае. Но для некоторых циркониевых руд не удается уверенно определить возраст.

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

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

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

Расширяющаяся дисциплина компрессионной физики

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

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

При помощи лазера удобно выравнивать аппаратуру и диагностические инструменты перед пуском Z-машины. В процессе сжатия рентгеновские лазеры зачастую используются для снятия данных. Иллюстрация: Sandia National LaboratoriesПри помощи лазера удобно выравнивать аппаратуру и диагностические инструменты перед пуском Z-машины. В процессе сжатия рентгеновские лазеры зачастую используются для снятия данных. Иллюстрация: Sandia National Laboratories

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

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

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

Благодарности

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

Подробнее..

От забоя до перегрузочного пункта. Пример интеграции GEOVIA Surpac и автоматизированной системы диспетчеризации ГТК

29.06.2020 12:20:00 | Автор: admin
image
Что добывают предприятия? Золото, железную руду, уголь, алмазы? Нет!
Каждое предприятие добывает деньги. Это и есть цель каждого предприятия. Если добытая тонна золота или железной руды не принесет вам доход или, хуже того, ваши затраты будут выше, чем прибыль от реализации продукции, какая ценность этой руды для предприятия?
Каждая тонна руды должна приносить максимальный доход или нести минимальные издержки в условиях безопасного производства и соблюдении технологии добычи. Т.е. распределение движения горной массы во времени должно приводить предприятие к цели. Для того чтобы достичь цели, необходимо создать хороший план, который будет моделировать производственный процесс с максимальным достижением объемно-качественных показателей. Каждый план необходимо обеспечивать верными, точными и актуальными данными. Особенно, если речь идет о краткосрочном или оперативном планировании.

Какие данные обеспечивают планирование горных работ? Это маркшейдерская и геологическая информация, данные проектирования и производственно-техническая информация (например, из ERP систем).
image

Все эти процессы несут в себе огромное количество графической, цифровой и текстовой информации, такой как облако точек лазерного сканирования, маркшейдерская база данных, оперативная съемка забоев, геологическая блочная модель, данные опробования скважин БВР, изменение контактов во взорванной горной массе, производственные показатели и их изменения, изменение динамики работы оборудования и т.д. Поток данных постоянен и бесконечен. И большая часть информации зависит друг от друга. Нужно не забывать, что все эти данные это первоисточник, информация, с которой начинается создание плана.
Поэтому для создания оптимального плана необходимо уметь получать максимально точный данные. Корректность исходной информации в геометрической прогрессии влияет на конечную цель.
image

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

Если говорить о краткосрочном планировании, то важно, чтобы эти данные были не только точными, но еще и актуальными. Необходимо иметь возможность получать информацию в любой момент времени, чтобы реагировать на изменения и оперативно редактировать сценарий производства. Соответственно, нужны такие системы и оборудование, которые позволят повысить эффективность процессов по получению и обработке информации. Лидарные сканеры позволяют оперативно получать данные с высокой точностью, технологии опробования горной массы дают картину положения рудного тела в массиве, системы позиционирования отслеживают положение и состояние оборудования в реальном времени, а системы проектирования и планирования GEOVIA Surpac и GEOVIA MineSched являются инструментами для создания проектов и сценариев развития горных работ. Для максимально быстрого достижения цели системы должны быть связаны в единую продуктивную цепочку. Представьте: вы получаете данные из разных систем и источников, но они доступны вам только по запросу, к тому же эти данные передает вам специалист, который в любой момент может изменить содержимое. Это ведет не только к снижению скорости получения данных, но и к потере точности или достоверности на одном из этапов передачи данных. Поэтому данные должны быть централизованными, храниться на одной платформе, в одной цифровой экосистеме и быть доступными в любой момент времени. Кроме того, важно обеспечить совместную работу всех подразделений, версионность, целостность и безопасность данных. С этой задачей справляется платформа 3DEXPEREINCE.

Информация, полученная из различных источников электронных систем, систем ГГИС (GEOVIA Surpac), ERP-систем, систем автоматизированного планирования горных работ (GEOVIA MineSched), систем управления горными работами (например, ВИСТ Групп) имеет разный формат данных.
Здесь встает вопрос об интеграции систем. Зачастую все решения в цепочке горного планирования и проектирования могут быть в большей или меньшей степени интегрированы между собой.
Но интенсивность потока данных, количество их типов, и изменчивость таковы, что человек не способен в относительно быстрое время произвести конвертацию из одной системы в другую. Будь то геолог или инженер по планированию, специалист должен тратить время не на импорт и экспорт файлов из одной системы в другую, он должен создавать ценность и двигать предприятие к цели. Поэтому процесс интеграции важно автоматизировать, настроить таким образом, чтобы количество манипуляций по обработке данных сводилось к минимуму.

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

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

Такой процесс реализован в маркшейдерской и геологической службах на Качканарском ГОКе компании ЕВРАЗ.

ЕВРАЗ КГОК входит в пятерку крупнейших в России горнорудных предприятий. Комбинат расположен в 140 км от ЕВРАЗ НТМК, в Свердловской области. ЕВРАЗ КГОК разрабатывает Гусевогорское месторождение титаномагнетитовых железных руд, содержащих примеси ванадия. Содержание ванадия позволяет выплавлять высокопрочные легированные сорта стали. Производственная мощность комбината составляет порядка 55 млн тонн железной руды в год. Основным потребителем продукции ЕВРАЗ КГОКа является ЕВРАЗ НТМК.

В настоящее время ЕВРАЗ КГОК добывает руду из четырех карьеров с дальнейшей ее переработкой в цехах дробления, обогащения, агломерации и окускования. Конечный продукт (агломерат и окатыши) загружается в железнодорожные вагоны и отправляется потребителям, в том числе за рубеж.

В 2018 году на ЕВРАЗ КГОКе добыто более 58,5 млн тонн руды, произведено 3,5 млн тонн агломерата, 6,5 млн тонн окатышей, около 2,5 млн тонн щебня.

Добыча руды производится в четырех карьерах: Главном, Западном, Северном, а также карьере Южная залежь. С нижних горизонтов руда доставляется БелАЗами, на дробильную фабрику горная масса перевозится железнодорожным транспортом. В карьерах используются мощные 130-тонные самосвалы, современные локомотивы НП-1, экскаваторы с объемом ковша 12 кубических метров.

Среднее содержание железа в руде составляет 15,6 %, содержание ванадия 0, 13 %.

Технология добычи железной руды на ЕВРАЗ КГОКе следующая: бурение взрывание экскавация транспортировка к месту переработки и вскрыши в отвалы. (Источник )

В 2019 году на Качканарском ГОКе была внедрена автоматизированная система диспетчеризации ВИСТ Групп. Внедрение данного решения позволило увеличить производственный контроль работы горнотранспортного оборудования, движения руды от забоев до перегрузочных пунктов, а также оперативно получать данные по объемно-качественным показателям в забоях и на перегрузочных пунктах. Была произведена двусторонняя интеграция систем АСД ВИСТ и GEOVIA Surpac, что позволило использовать получаемые данные (положение оборудования, степень отработки забоя, баланс горной массы на перегрузочных пунктов, распределение качества на перегрузочных пунктах и т.п.) для оперативного планирования и проектирования горных работ, а также контролировать производственный процесс на уровне линейного руководителя и машиниста экскаватора.
image

Благодаря разработкам ведущего геолога С.М. Некрасова и главного маркшейдера А.В. Безденежных, специалистами маркшейдерского и геологического отделов с использованием инструментов GEOVA Surpac было автоматизировано большинство процессов по обработке данных маркшейдерской съемки, проектирования, создания печатной документации, создания геологических блочных моделей, обновления геологической и маркшейдерской информации на сетевом ресурсе. Теперь специалистам не нужно выполнять повторяющиеся процессы ежедневно, будь то выгрузка/загрузка съемки с прибора/на прибор, поиск нужных данных для повседневной работы в огромном множестве папок. Макрокоманды GEOVIA Surpac делают это за них. Важно отметить, что эти данные доступны для всех причастных специалистов разных отделов. Например, чтобы открыть последнюю съемку карьера, обновленную блочную модель, блок БВР, коммуникации и т.п., специалисту по планированию не нужно искать это в большом числе маркшейдерских и геологических файлов. Все, что ему требуется для этого открыть в GEOVA Surpac соответствующее меню и выбрать данные, которые нужно загрузить в рабочее окно.
image

Инструменты автоматизации позволили легко произвести интеграцию GEOVIA Surpac и АСД ВИСТ Групп и сделать этот процесс максимально простым и быстрым.

Выбрав соответствующее меню в панели GEOVIA Surpac, геолог получает из АСД ВИСТ оперативные данные по отработке блока или данные на определенную дату и время. Эти данные используются для анализа текущей ситуации и обновления блочной модели.

image

После обновления блочной модели и контактов руда/вскрыша в GEOVIA Surpac геолог нажатием одной кнопки загружает эту информацию в систему АСД ВИСТ, после чего данные доступны для всех пользователей в обеих системах.

image

image

Благодаря объединению возможности инструментария позиционирования горнотранспортного оборудования в системе АСД ВИСТ Групп и инструментов GEOVIA Surpac были настроены процессы контроля движения горной массы от забоя до перегрузочного пункта, размещения горной массы в секторах перегрузочных пунктов, контроля баланса прихода/ухода горной массы по секторам и ведения мобильных остатков за период оперативного заполнения.
Для этого в GEOVIA Surpac были созданы блочные модели перегрузочных пунктов и разработана методика по их заполнению. По желанию геолога, процесс внесения в блочную модель (БМ) горной массы на виртуальный перегрузочный пункт, как и отгрузку с него, можно осуществлять как целиком за прошедший период, так и в оперативном режиме. Поставив на заполнение БМ с указанием времени окончания, макропрограмма сама производит запрос (через определенный интервал времени) на извлечение данных по экскаваторам, производящих черпание, а также извлекает информацию по движению и выгрузке автотранспорта на перегрузочный пункт.
Таким образом, по окончании работы макропрограммы формируется текущая информация по состоянию склада, наличию на данный период времени горной массы в трехмерном графическом виде и сводная таблица результатов оперативного изменения. Это позволило оперативно отслеживать движение руды, баланс и распределение горной массы на секторах перегрузочных пунктов, а также представлять графически эту информацию в обеих системах и обеспечить быстрый, свободный и безопасный доступ к информации для всех сотрудников. В частности, по словам ведущего геолога С.Н. Некрасова, такой процесс позволил повысить точность планирования отгрузки по качеству с перегрузочных пунктов на железнодорожный транспорт.
Также он отмечает, что если раньше можно было лишь предполагать, что было привезено на перегрузочные пункты и представлять только среднее значение качества по секторам, то сегодня известны показатели на каждом отдельно взятом участке сектора.

image

image

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

image

Подробнее о процессе и результатах интеграции, проведенной на Качканарском ГОКе, вы можете узнать из записи веб-семинара Новый подход к автоматизации планирования, ведения буровзрывных работ и управлению качеством на предприятии по ссылке

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

Подписывайтесь на новости Dassault Systmes и всегда будьте в курсе инноваций и современных технологий.

Dassault Systmes официальная страница

Facebook
Vkontakte
Linkedin
3DS Blog WordPress
3DS Blog on Render
3DS Blog on Habr
Подробнее..

Нефтянка для инженеров, программистов, математиков и широких масс трудящихся, часть 3

19.06.2020 08:12:45 | Автор: admin
image

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


Первую и вторую части серии можно прочесть здесь:
habr.com/ru/company/bashnipineft/blog/505300
habr.com/ru/company/bashnipineft/blog/506198

Осадконакопление, образование нефти, миграция

Теорий образования нефти есть несколько, я буду говорить только о той, согласно которой нефть образовалась из остатков живых существ зоопланктона и водорослей. Когда-то давным-давно, in a galaxy far far away, этот самый зоопланктон и водоросли размножался так бурно, а условия по наличию (точнее, отсутствию) кислорода и нужной температуры были такие, что он не успевал разложиться и падал на дно моря, будучи в дальнейшем засыпаем всякой осадочной породой (песочком и глиной). Сейчас даже деревья отказываются превращаться в каменный уголь, а просто, заразы, в присутствии кислорода гниют, а тогда они аналогично планктону подвергались захоронению и с течением времени погружались на всё большую и большую глубину, уплотнялись и обезвоживались. Это не означает, кстати, что вода куда-то пропадала вода всегда была и оставалась вокруг, пропитывая породу. Вода уходила из самих остатков клеток и органических соединений.

image

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

На пути к поверхности нефть встречает много разных слоёв, и они к этому времени совсем не обязательно будут горизонтальными. Река подточила глиняный утёс, он упал в воду и его смыло, вынесло в море, и там вся эта глина упала на дно в виде протяжённого, более или менее горизонтального пятна. Процессы движения литосферы, процессы горообразования вспучили земную твердь, и вот уже когда-то горизонтальный пласт низкопроницаемой глины выпятило в виде арки ловушка для нефти готова. Где-то нефть с газом попадёт в такие ловушки и застрянет на определённой глубине, а где-то просочится практически на поверхность. В последнем случае газ и лёгкие углеводороды, конечно, большей частью опередят всех остальных, достигнут поверхности и улетучатся, и останутся битуминозные пески Альберты и бассейна реки Ориноко, но в качестве противоположных примеров предлагаю вспомнить бакинские колодцы и озера с нефтью в 19 веке, а также горящие уже тысячи лет огни Химеры выходы метана на поверхность в районе турецкой горы Олимпос.

image

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

image

Динамика свойств и технологических показателей в процессе работы скважины

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

Какое было начальное пластовое давление? Если забыть про все эти ловушки для нефти и вообразить, будто вся земля представляет собой проницаемый песчаник, насыщенный водой, то давление жидкости на любой глубине будет определяться только плотностью этой воды и глубиной: это известная со школы формула ро-жэ-аш. Делайте как мы, забудьте начальную школу с её задачками на перевод из дециакров в кубофуты и считайте всё только в единицах СИ, и будет вам счастье. Получаем примерно 1000*10*2000 Паскалей или примерно 200 атмосфер.

image

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

В случае повышенного пластового давления воду из скважины будет выталкивать со всё повышающейся скоростью, и наконец из скважины забьёт фонтан нефти! Теперь представим, что у месторождения нормальное начальное пластовое давление, не повышенное, но скважина заполнена не водой, а нефтью. У нефти плотность меньше, чем у воды, поэтому ро-жэ-аш у такого столба нефти будет меньше, чем пластовое давление, и нефть из скважины будет выталкивать! Чтобы остановить такую скважину, нужно её заглушить, заменить нефть в стволе скважины на воду, а при аномально высоком пластовом давлении ещё и с растворенными солями, чтобы плотность была побольше и ро-жэ-аш как раз равнялось пластовому давлению (и чуть-чуть его превышало). Сами понимаете, при аномально высоком пластовом давлении ещё и при бурении будут проблемы нужно пробурить скважину, создав открытую связь месторождения под давлением с поверхностью, но при этом не допустив выбросов из-за повышенного давления.

image

Что происходит с физико-химическими свойствами месторождения при добыче? В первом приближении можно считать, что пористость, проницаемость и сжимаемость не меняются, потому что порода остаётся на месте, а вот все насыщающие её флюиды начинают мигрировать, поэтому меняются давление и насыщение. Пусть скважина была заполнена нефтью и при этом находилась в равновесии с пластовым давлением. Насос начинает добывать нефть из скважины, динамический уровень нефти в скважине начинает падать, забойное давление (давление столба жидкости по ро-жэ-аш на забой, то есть нижнюю часть скважины) тоже падает. Насос может быть настроен так, чтобы работать с заданным постоянным расходом, или (если в нём есть свой манометр) поддерживать заданное постоянное забойное давление. Перепад давления (называемый депрессией) между забойным (в стволе скважины) и пластовым (в прилегающей к скважине призабойной части месторождения) давлением вызывает приток нефти к скважине. Но и одновременно при этом понижает пластовое давление в призабойной зоне месторождения! Если проницаемость достаточная, то волна падения давления распространяется дальше, вызывая приток нефти из всё более дальних частей месторождения. Теоретически, при большой проницаемости можно всю нефть месторождения выкачать одной скважиной, просто это займёт очень много времени.

image

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

Итак, после начала добычи и закачки на месторождении начинается невиданная за тысячи лет движуха перераспределяется пластовое давление, вызывая фильтрацию (напоминаю, слово фильтрация можно смело заменять на перетоки) содержимого от мест с большим давлением к местам с меньшим, меняя, таким образом, во всех точках нефтенасыщенность, водонасыщенность и газонасыщенность. Если вдруг добычу остановить, то через некоторое время пластовое давление выровняется, на месторождении останется только одна сила тяжести и снова начнутся те же медленные процессы восстановления равновесия: более тяжёлая вода вниз, более лёгкий газ вверх, нефть между ними. Коктейль снова расслаивается.

Физико-химические свойства нефти, воды и газа

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

Но оказывается, когда они все находятся в одном месте, в одних и тех же порах, они друг другу начинают мешать, и скорость фильтрации их всех зависит не только от абсолютной проницаемости породы и собственной вязкости, но и от количества остальных флюидов в наличии. Так, например, опыты показывают, что когда в одних и тех же порах находятся нефть и вода, скорость фильтрации нефти при соотношении 10:1 может оказаться во много раз больше, чем скорость фильтрации нефти при соотношении 1:10.

image

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

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

image

В итоге, есть абсолютная проницаемость породы, которая зависит только от породы, и есть относительная проницаемость для воды, нефти, газа, которая зависит и от их пропорции, и это всё отображается на графиках, подобных тому, что я привёл выше. Как такую картинку читать? Водонасыщенности меньше 20% на этом гипотетическом месторождении нет, то есть воды всегда как минимум 20%. При таком содержании воды её относительная фазовая проницаемость равна нулю, то есть вода есть, её 20%, но она не течёт никуда, течёт только нефть. С другой стороны, когда воды становится 75% (а нефти остаётся, соответственно, только 25%), относительная фазовая проницаемость нефти падает до нуля, и значит уже нефть никуда больше не течёт, и эти оставшиеся 25% из пор никак не выковырять.

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

Ещё одна особенность нефти заключается в её способности растворять углеводородный газ (чаще всего метан). В одном кубометре нефти могут быть растворены десятки и сотни кубометров газа, в зависимости от давления. Отношение объёма газа, выделившегося из нефти при добыче её на белый свет и в нормальные условия, к объёму той нефти, из которой он выделился, называется газосодержанием. Если вы добываете нефть с газосодержанием, равным 100, то при добыче 1 кубометра нефти вы получите попутно ещё и 100 кубометров газа, которые в ней были растворены и которые из неё при добыче выделятся.

image

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

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

image

В случае месторождения всё происходит похожим образом, но в разных направлениях в зависимости от того, что было в начале и что происходит. Если у месторождения была газовая шапка при давлении в 200 атмосфер, это значит, что нефть уже приняла весь газ, который могла растворить, и находится в насыщенном состоянии. Если начать закачивать в месторождение воду и повышать пластовое давление, то растворимость газа (способность нефти вмещать газ) начинает расти, и газ тут же пользуется этой возможностью и в нефти начинает растворяться (правда, происходит это не быстро). Так что повышая давление можно теоретически добиться, чтобы вся газовая шапка на месторождении растворилась в нефти. Допустим, что это произойдёт при 250 атмосферах, а мы продолжаем нагнетать, что будет происходить, когда пластовое давление вырастет, скажем, до 300 атмосфера? Да ничего, газа как не было, так и не будет. Но нефть как бы запомнила, при каком давлении в ней растворился последний кубометр газа, это давление и называется давлением насыщения.
Если затем начать добывать нефть и снижать давление с 300 до 290, 280, 270 атмосфер, то ничего происходить не будет, но как только давление упадёт до давления насыщения 250 атмосфер, газ начнёт снова выделяться, причём не в шапку, а по всему объёму, и быстро.
Если у месторождения с пластовым давлением в 200 атмосфер с самого начала не было газовой шапки, означает ли это, что в этой нефти нет растворённого газа? Нет, не означает растворённый газ всегда в нефти есть, просто из-за условий её образования. Для такой нефти важно знать, какое у неё газосодержание и давление насыщения, потому что если вдруг в процессе разработки пластовое давление упадёт ниже давления насыщения, то по всему объёму месторождения из нефти начнёт выделяться газ, и этот газ будет очень сильно мешать нормальной фильтрации нефти к скважинам.

image

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

Кроме этого, важно знать, что при растворении в нефти газа, объём этой самой нефти увеличивается. И наоборот, при добыче 1 кубометра нефти из пласта, когда мы её поднимем на поверхность и из неё выделится растворённый газ, её объем уменьшится в соответствии с так называемым объёмным коэффициентом. Когда нефть поднимают на поверхность, происходит сразу несколько явлений, в разном направлении влияющих на её объём: падает давление (нефть увеличивает объём в соответствии со своим коэффициентом сжимаемости), падает температура (нефть уменьшает объём в соответствии со своим коэффициентом температурного расширения), улетучивается растворённый газ (объём и масса нефти уменьшаются). Всё это вместе объединяют в единый объёмный коэффициент. Если он равен, например, 1.2, то чтобы получить 1 кубометр на поверхности, нужно забрать 1.2 кубометра из пласта.

И вязкость, и объёмный коэффициент, конечно, зависят от давления и газосодержания: чем больше газосодержание (при одном и том же давлении), тем больше объёмный коэффициент и тем меньше вязкость. Чем больше давление (при одном и том же газосодержании), тем больше вязкость и меньше объёмный коэффициент.

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

Продолжение следует
Подробнее..

Нефтянка для инженеров, программистов, математиков и широких масс трудящихся, часть 4

02.07.2020 08:05:44 | Автор: admin
image

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

Сегодня мы расскажем о том, зачем нужны модели месторождений, и как их строить. Модель это тот самый план действий, который обязательно должен быть и предполагаемый результат этих действий.


Моделирование, прогноз, неопределённость


Все перечисленные в предыдущих статьях (раз, два, три) физические эффекты важно понимать не просто для того, чтобы знать, как устроен мир. Скорее всего их придётся учитывать при построении модели, которая умеет правильно предсказывать будущее. Зачем нам уметь предсказывать будущее в нефтедобыче, если цену нефти и коронавирус всё равно не предсказать? Да затем же, зачем и везде: чтобы принимать правильные решения.

image

В случае с месторождением мы не можем непосредственно наблюдать, что происходит под землёй между скважинами. Почти всё, что нам доступно, привязано к скважинам, то есть к редким точкам на необъятных просторах болот (все что мы можем измерить заключено примерно в 0.5% породы, о свойствах остальных 99.5% мы можем только догадываться). Это проведённые измерения на скважинах тогда, когда скважина строилась. Это показания приборов, которые на скважинах установлены (давление на забое, доля нефти воды и газа в получаемой продукции). И это измеряемые и задаваемые параметры работы скважин когда включать, когда выключать, с какой скоростью качать.

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

image

Если модель не может воспроизвести историю, необходимо менять её входные данные, но какие? Фактические данные менять нельзя: это результат наблюдения и измерения реальности данные с приборов. Приборы, конечно, имеют свою погрешность, да и используются приборы людьми, которые тоже могут и накосячить, и приврать, но неопределённость фактических данных в модели, как правило, мала. Менять можно и нужно то, что имеет наибольшую неопределённость: наши предположения о том, что происходит между скважинами. В этом смысле построение модели это попытка уменьшить неопределённость в наших знаниях о реальности (в математике этот процесс известен как решение обратной задачи, и обратных задач в нашей области как велосипедов в Пекине!).

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

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

Модели в жизни месторождения


Для того, чтобы принимать решения в процессе разработки месторождения, нужна целостная модель всего месторождения. Более того, сейчас без такой модели разрабатывать месторождение вообще нельзя: такую модель требуют государственные органы РФ.

image

Всё начинается с сейсмической модели, которая создаётся по результатам сейсморазведки. Такая модель позволяет увидеть под землей трехмерные поверхности специфические слои, от которых хорошо отражаются сейсмические волны. Она не даёт почти никакой информации о нужных нам свойствах (пористости, проницаемости, насыщении, и т. д.), но зато показывает, как изгибаются в пространстве некоторые слои. Если вы сделали многослойный бутерброд, а потом как-то его изогнули (ну или кто-то на него сел), то у вас есть все основания считать, что все слои изогнулись примерно одинаково. Поэтому мы можем понять, как изогнулся слоёный пирог из различных нападавших на дно океана осадков, даже если на сейсмической модели увидим только один из слоёв, по счастливой случайности хорошо отражающий сейсмические волны. На этом месте дата-сайнс инженеры оживились, потому что автоматическое выделение таких отражающих горизонтов в кубе, чем и занимались участники одного из наших хакатонов, классическая задача распознавания образов.

image

Затем начинается разведочное бурение, и по мере бурения скважин в них спускают на кабеле приборы, измеряющие всякие разные показатели вдоль ствола скважины, то есть проводят ГИС (геофизические исследования скважин). Результат такого исследования каротаж ГИС, то есть кривая определённой физической величины, измеренная с определённым шагом вдоль всего ствола скважины. Разные приборы измеряют разные величины, а обученные инженеры затем проводят интерпретацию этих кривых, получая значимую информацию. Один прибор измеряет естественную гамма-радиоактивность породы. Глины фонят сильнее, песчаник фонит слабее это знает любой инженер-интерпретатор и выделяет на каротажной кривой: тут глины, тут слой песчаника, тут что-то среднее. Другой прибор измеряет естественный электрический потенциал между соседними точками, возникающий при проникновении в породу бурового раствора. Высокий потенциал показывает наличие фильтрационной связи между точками пласта, знает инженер и подтверждает наличие проницаемой породы. Третий прибор измеряет сопротивление насыщающего породу флюида: солёная вода ток пропускает, нефть ток не пропускает и позволяет отделить нефтенасыщенные породы от водонасыщенных и так далее.
На этом месте дата-сайнс инженеры снова оживились, потому что входные данные у этой задачи это простые численные кривые, а заменить инженера-интерпретатора какой-нибудь ML-моделью, умеющей вместо инженера по форме кривой сделать вывод о свойствах породы это значит решить классическую задачу классификации. Это только потом у дата-сайнс инженеров начинает дергаться глаз, когда выясняется, что часть этих накопленных кривых со старых скважин есть только в виде длинных бумажных портянок.

image

Кроме этого, при бурении из скважины достают керн образцы более или менее целой (если повезло) и неразрушенной при бурении породы. Эти образцы отправляют в лабораторию, где определят их пористость, проницаемость, насыщение и всякие разные механические свойства. Если известно (а при правильном проведении это должно быть известно), с какой глубины был поднят конкретный образец керна, то когда придут данные из лаборатории, можно будет сопоставить, какие значения на этой глубине показывали все геофизические приборы, и какие значения пористости, проницаемости и насыщенности имела порода на этой глубине по данным лабораторных исследований керна. Таким образом можно пристрелять показания геофизических приборов и затем только по их данным, не имея керна, делать вывод о таких нужных нам для построения модели свойствах породы. Весь дьявол в деталях: приборы замеряют не совсем то, что определяют в лаборатории, но это уже совсем другая история.

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

image

На скважинах мы точно знаем, какие свойства имеет порода на какой глубине, но не знаем, как слои породы, наблюдаемые на скважинах, распространяются и изгибаются между ними. Сейсмическая модель не позволяет точно определить, какой слой на какой глубине расположен, но зато уверенно показывает характер распространения и изгиба всех сразу слоёв, характер напластования. Тогда инженеры отмечают на скважинах определённые характерные точки, ставя на определённой глубине маркеры: на этой скважине на этой глубине кровля пласта, на этой глубине подошва. А поверхность кровли и подошвы между скважинами, грубо говоря, рисуют параллельно той поверхности, которую видят в сейсмической модели. В итоге получается набор трёхмерных поверхностей, которые охватывают в пространстве интересующие нас, а нас интересуют, конечно же, пласты, содержащие нефть. То, что получилось, называется структурной моделью, потому что она описывает структуру пласта, но не его внутреннее содержание. О пористости и проницаемости, насыщении и давлении внутри пласта структурная модель не говорит ничего.

image

Затем наступает этап дискретизации, при котором область пространства, занимаемая месторождением, разбивается на такой изогнутый в соответствии с залеганием слоёв (характер которого виден ещё на сейсмической модели!) параллелепипед из ячеек. Каждая ячейка этого изогнутого параллелепипеда однозначно определяется тремя номерами, I, J и K. Все слои этого изогнутого параллелепипеда лежат согласно распространению слоёв, а количество слоёв по K и количество ячеек по I и J определяется детальностью, которую мы можем себе позволить.
Насколько детальная информация о породе у нас есть вдоль ствола скважины, то есть по вертикали? Настолько детальная, насколько часто делал замеры своей величины геофизический прибор при движении по стволу скважины, то есть, как правило, каждые 20-40 см, поэтому каждый слой может быть и 40 см, и 1 м.

Насколько детальная у нас есть информация по латерали, то есть в сторону от скважины? Ни насколько: в сторону от скважины у нас информации нет, поэтому смысла разбивать на очень маленькие ячейки по I и J, как правило, нет, и чаще всего они бывают по 50 или 100 м по обеим координатам. Выбор размера этих ячеек является одной из важных инженерных задач.

image

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

image

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

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

image

Затем за работу берутся специалисты-гидродинамики. Они могут укрупнить геологическую модель, объединив несколько слоёв по вертикали и пересчитав свойства породы (это называется апскейлинг, и представляет собой отдельную непростую задачу). Потом они добавляют остальные нужные свойства для того, чтобы гидродинамический симулятор мог моделировать, что куда будет перетекать: кроме пористости, проницаемости, нефте-, водо-, газонасыщенности, это будут давление, газосодержание и так далее. Они добавят в модель скважины и внесут по ним информацию о том, когда и в каком режиме они работали. Вы ещё не забыли, что мы пытаемся воспроизвести историю, чтобы иметь надежду на корректный прогноз? Гидродинамики возьмут отчёты из лаборатории и добавят в модель физико-химические свойства нефти, воды, газа и породы, всяческие зависимости их (чаще всего от давления) и всё, что получилось, а это будет гидродинамическая модель, отправят в гидродинамический симулятор. Тот честно рассчитает, из какой ячейки в какую всё будет в какой момент времени перетекать, выдаст графики технологических показателей на каждой скважине и скрупулёзно сравнит их с реальными историческими данными. Гидродинамик вздохнёт, глядя на их расхождение, и пойдёт изменять всё неопределённые параметры, которые он пытается угадать так, чтобы при следующем запуске симулятора получить что-то близкое к реально наблюдавшимся данным. А может при следующем запуске. А может при следующем и так далее.

image

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

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

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

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

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

Продолжение следует
Подробнее..

Нефтянка для инженеров, программистов, математиков и широких масс трудящихся, часть 5

06.08.2020 08:23:05 | Автор: admin

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

Симулятор гидроразрыва пласта





Одной из самых важных на сегодня операций, проводимых на скважинах, является гидравлический разрыв пласта (ГРП). Суть гидроразрыва пласта в следующем. В скважину под большим давлением (до тысячи атмосфер!) закачивают специальную жидкость, похожую на желе (собственно, это и есть желе при его производстве используют пищевой гелеобразователь). Давление разрывает пласт, раздвигая слои породы. На той глубине, где обычно производится ГРП, породу проще раздвинуть в стороны, чем вверх, поэтому трещина получается почти плоская и вертикальная, при этом ширина её составляет считанные миллиметры или сантиметры. Затем вместе с жидкостью начинает подаваться пропант (иногда также пишут проппант, с двумя п и так, и так правильно) похожая на песок смесь крепких керамических гранул диаметром от долей миллиметров до миллиметров. Цель ГРП закачать нужное количество пропанта в пласт (а, точнее, в коллектор в ту часть пласта, где находится нефть) так, чтобы образовалась проницаемая область, соединённая со скважиной. Жидкость, конечно, утечёт в пласт, а пропант останется там, куда успел дойти. Низкопроницаемые месторождения (а сейчас почти все новые месторождения именно такие, все большие хорошие традиционные месторождения давно уже нашли и разбурили) бесполезно бурить обычными скважинами: проницаемость мала, и к маленькому стволу скважины нефть еле сочится. Можно, конечно, теоретически и скважину пробурить очень толстую, чтобы площадь стенок у неё была большая, но очевидно, что на практике это сделать нереально. А вот если сделать на скважине ГРП, то вокруг скважины образуется обширная, хорошо соединённая со скважиной трещина с большой площадью стенок. Нефть сочится так же медленно, как и раньше, но теперь скважина собирает её потихоньку с большой площади, и суммарный дебит получается хороший, и скважина экономически оправдана. Чтобы ещё увеличить площадь, с которой собирается одной скважиной нефть, скважины бурят горизонтально, и проводят ГРП несколько раз в разных местах горизонтального ствола.



Что тут может пойти не так, и зачем нужно моделирование? Во-первых, трещина может пойти не туда, куда нужно, и пропант может распределиться по ней не так, как нам бы хотелось. Месторождение могло разрабатываться уже давно, и часть слоёв породы уже обводнилась, насытилась водой, а, возможно, и с самого начала на месторождении нижние слои были водонасыщены. Если вдруг трещина пройдёт от скважины к водонасыщенным или газонасыщенным пластам и заполнится пропантом, то вместо части нефти скважина будет добывать воду или газ.

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

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



Как можно смоделировать то, что в скважине будет происходить при закачке под большим давлением жидкости и пропанта? Этим занимается специализированный софт под названием симулятор ГРП, который использует численную физико-математическую модель развития трещины. Разработка такого софта, как и разработка любого другого софта для численного моделирования физических процессов, требует одновременного участия физиков, математиков, программистов высокопроизводительных параллельных вычислений и программистов, умеющих создавать удобный интерфейс, потому что в конечном счёте симулятором ГРП будут пользоваться инженеры в болотах Западной Сибири, нещадно кусаемые комарами!

Программы, которые позволяют рассчитать ГРП, буквально по пальцам сосчитать, и все они, до недавнего времени, были сделаны только в одной стране мира и легко запрещались к продаже в России. Да, мы гордимся, что разработали и продолжаем развивать первый в нашей стране промышленный симулятор гидроразрыва пласта, РН-ГРИД. До нашей разработки не только мы, но и вся отечественная индустрия была вынуждена пользоваться только американским ПО, потому что другого, собственно, не было, а с 2014 года с покупкой этого ПО возникли проблемы (как говорится, не было бы счастья, да несчастье помогло). Сейчас РН-ГРИД не только полностью заменил всё аналогичное ПО внутри Роснефти, но и успешно продаётся сторонним компаниям.



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



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



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

В итоге получится самый быстрый и самый навороченный в мире консольный exe-шник. Чтобы подать ему на вход все нужные данные, и красиво отобразить на выходе результаты расчёта, а потом еще и выгрузить их в отчет по нужной форме, понадобится ещё и такое удобное приложение, как пре-пост-процессор. Полевые инженеры, напомню, довольно суровые ребята, которые сидят ночами, интерпретируя графики и пересчитывая дизайны ГРП. Они всегда готовы высказать всё, что думают о режущем глаза шрифте или о кнопке в диалоге, которая уезжает за пределы экрана на старых ноутбучных мониторах 1024 на 768, которые ты, программист, развалясь в своём уютном офисе за парой 4К мониторов, забыл предусмотреть. А так как работа идёт и днём, и ночью, они тоже очень любят тёмную тему:



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

Симулятор операций с гибкой трубой



В хирургии с незапамятных времён привыкли к операциям, выполняемым через вену, а медицинские эндоскопы тоже давно и успешно используются для исследования и манипуляций со всякими внутренними органами человека. В нефтянке есть свой аналог этой операции это гибкая насосно-компрессорная труба (ГНКТ) меньшего диаметра, которая спускается внутрь обычной насосно-компрессорной трубы и позволяет выполнять разного рода работы, которые иначе выполнять сложно.



Как используют гибкую трубу? Основных вариантов два. Это может быть просто гибкая труба с открытым концом. Её спускают в обычную трубу, по которой до того осуществлялась закачка или добыча, и начинают через неё закачивать жидкость и газ. Но при этом хотят, чтобы они не уходили в пласт, а поднимались по трубе обратно. Посмотрите на схему, и вам станет понятна идея: в трубу 1 засунута труба 2, по трубе 2 подаётся жидкость и газ, а между трубой 1 и трубой 2 закачанные жидкость и газ поднимаются обратно на поверхность. Зачем это нужно? Например, если основная труба 1 была в результате предыдущей неудачной операции забита до половины пропантом, её таким образом можно промыть: пропант будет с обратным током жидкости и газа выноситься обратно на поверхность это называется промывка ствола. Например, если соединение скважины с пластом забилось и плохо пропускает нефть, можно закачивать вместе с жидкостью больше азота, тогда пластовое давление жидкости само прочистит призабойную область. Управляя давлением на входе в трубу 2 и на выходе из трубы 1 можно добиваться, чтобы циркулирующая жидкость не уходила в пласт, потому что нам нужно промыть скважину а можно, чтобы, наоборот, уходила, если мы закачиваем кислоту, чтобы подрастворить породу или нежелательные загрязнения вокруг скважины.

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



Мы в нашем институте сейчас разрабатываем наше ПО для проведения расчётов различных операций, производимых через гибкую трубу РН-ВЕКТОР. Угадаете, где сделан массово применяемый на просторах РФ аналог, который сейчас приходится использовать инженерам?

Что здесь нужно моделировать? Во-первых, нагрузки на гибкую трубу. Её будут спускать в другую трубу, она будет там тереться о стенки, и сначала её нужно будет заталкивать с дополнительным усилием, а потом наоборот, удерживать. Но даже когда в скважину будет спущено 2-3 километра трубы, все эти 2-3 километра будут находиться в очень разных условиях: верхняя часть будет висеть и под собственным весом растягиваться, а нижняя часть, ушедшая в горизонтальный ствол, будет лежать и стараться застрять. Соответственно, важно рассчитать, чтобы нигде и не застряло, и не порвалось. Расчёт разнообразных нагрузок в соответствии с весом гибкой трубы и её трением о стенки основной трубы, гулянием трубы по спирали и механическими свойствами стали это первая часть любого такого продукта.

Во-вторых, гидравлика. На входе в гибкую трубу одно давление, которое теряется сначала при движении жидкости по намотанной на бобину трубе, потом теряется при движении по гибкой трубе вниз, превращается в какое-то давление внизу. Если вдруг внизу окажется давление больше, чем в пласте, жидкость будет уходить в пласт, теряться, вместо того, чтобы уходить наверх между трубами 1 и 2 на рисунке выше. Уважаемый мастер бригады Иван Иванович, скажет рабочий бригады ГНКТ, несомненно, мы наблюдаем потерю циркуляции, как Вы считаете?. Если наоборот, окажется меньше, то из пласта будет сквозить нефть, а мы же ремонт на скважине делаем, а не добычей занимаемся. Так что симулятор этой технологической операции должен уметь правильно рассчитывать все перепады давления. Кстати, а давление в гибкой трубе обратно влияет на её нагруженное состояние.



В-третьих, вынос частиц. На рисунке выше представьте себе взбаламученный песок на дне скважины: будет ли он потоком жидкости подниматься наверх, или скорости потока не хватит, и он будет падать обратно вниз? А если закачивать жидкость вместе с азотом азотный пузырь быстрее воду наверх выталкивает, особенно расширяясь на пути снизу вверх, но при этом сам песок выносить не может. Как много взбаламученного песка можно вынести наверх и при этом не уронить пачку? Нужен расчёт, причём с разными режимами течения.

В-четвертых, нужен учёт усталости трубы. Гибкая труба распрямляется с барабана, снова сгибается на гузнеке, потом снова распрямляется на входе в скважину а потом всё в обратном порядке. Трубу то спускают, то поднимают, соответственно, какие-то её части сгибаются больше, какие-то меньше. Сталь трубы рассчитана на определённую накопленную усталость, и нужно от работы к работе рассчитывать очередную добавку усталости каждого метра трубы и хранить этот профиль в базе данных. Вот как выглядит этот профиль, после определённого числа работ:



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

Графики, графики, графики!



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



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



Но мы отклонились от темы. Вот, например, как может выглядеть в том же самом софте записанная операция ГРП:



Смотрит грамотный инженер на эти графики, и вся картина произошедшего разворачивается перед ним, как линии судьбы на ладони. В точке 1 началась закачка, и забойное давление (давление внизу скважины) начинает резко расти с точки 2 до точки 3, пока, наконец, при давлении 380 атмосфер не открывается трещина гидроразрыва. Обратите внимание, что давление внизу у скважины так и останется почти постоянным, пока трещина будет расти вплоть до точки 11. Кстати, когда трещина открывается при 380 атмосферах, на поверхности манометр показывает более 500 атмосфер в точке 4. Давление в трещине остаётся почти постоянным, а на поверхности показания манометра падают из точки 4 в точку 5. Инженер и глазом не поведёт: он отлично знает, что это почти вода в скважине замещается на тот самый гель, и потери давления на трение в трубе падают именно на разницу между точками 4 и 5. Пытливый инженер даже измерит наклон линии от точки 4 к точке 5, и получит таким образом отношение коэффициентов трения у жидкости, которая была в скважине и которая туда поступает.

В точке 6 начинает подаваться пропант, и смотрите как сразу в точке 7 начинают падать показания манометра на поверхности это столб жидкости в скважине из-за пропанта становится всё тяжелее и тяжелее. В точке 8 до точки 9 в скважину подбавили понизитель трения, чтобы пропант не слишком сильно тормозился о стенки трубы. В точке 10 перестали подавать пропант, он перестал под своим весом проваливаться в скважину, и поэтому сильнее приходится давлением его продавливать в трещину, устьевое давление возрастает до точки 11. А там насосы выключают, и давление мгновенно падает к точке 12, и там уже начинает медленно снижаться по мере того, как из трещины гидроразрыва утекает жизнь жидкость, и она постепенно закрывается.

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

Перевод Учёные раскрыли универсальную геометрию геологии, и оказалось, что мир состоит из кубов

18.12.2020 14:16:13 | Автор: admin

Упражнения в чистой математике привели к созданию масштабной теории об устройстве мира




Где-то в середине лета 2016 года венгерский математик Габор Домокош взошёл на крыльцо дома Дугласа Джерольмака, геофизика из Филадельфии. С собой у Домокоша были дорожные чемоданы, сильная простуда и жгучая тайна.

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

Сколько граней у каждого из этих камушков? спросил он. Затем он ухмыльнулся. Что, если я скажу вам, что их количество обычно равно шести? А затем он задал ещё более общий вопрос, который, как он надеялся, надолго поселится в мозге его коллеги. Что, если мир состоит из кубов?

Джерольмак сначала возразил: может, дома и строятся из кирпичей, но Земля состоит из камней. А форма у камней, очевидно, разная. Слюда крошится на чешуйки, кристаллы ломаются по жёстко определённым осям. Однако Домокош утверждал, что из одной лишь чистой математики следует, что любые камни, ломающиеся случайным образом, будут порождать формы, имеющие в среднем по шесть граней и восемь вершин. Если взять среднее по всем им, оно будет стремиться к некоему идеальному кубу. Домокош сказал, что доказал это математически. Теперь ему нужно было, чтобы Джерольмак помог ему показать, что это происходит и в природе.

Это было чёткое геометрическое предсказание, порождённое природой, причём без какой-либо физики, сказал Джерольмак, профессор из Пенсильванского университета. Как, чёрт возьми, природа это вообще допустила?

В последовавшие несколько лет парочка изучала свою геометрическую идею, исследуя всё, от микроскопических фрагментов камней до обнажения геологических пород, поверхностей планет и даже диалога Платона "Тимей". Всё это покрывало проект налётом мистицизма. Один из величайших философов примерно в 360 году до н.э. сопоставил пять платоновых тел с пятью элементами мироздания: землёй, воздухом, огнём, водой и звёздной материей. По удачному совпадению и/или предвидению Платон сопоставил кубы, которые лучше всего складываются в штабеля, с землёй. И я подумал ладно, вот мы уже слегка зашли и на территорию метафизики, сказал Джерольмак.


Габор Домокош и Дуглас Джерольмак

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

Несколько геофизиков, с которыми связалась редакция журнала, говорят, что ту же математическую платформу можно использовать и для других задач, вроде изучения эрозии разломов скал, или предотвращения опасных оползней. Это очень интересно, сказал геоморфолог Микаэль Аттал из Эдинбургского университета, один из двух рецензентов этой работы. Другой рецензент, геофизик Дэвид Фёрбиш из университета Вандербильта, сказал: Подобные работы заставляют меня задуматься: не могу ли я как-то воспользоваться этими идеями?

Все возможные разломы


Задолго до визита в Филадельфию у Домокоша возник более безобидный математический вопрос.

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

Сначала Домокошу стало интересно, можно ли только средствами геометрии предсказать, из каких фигур в среднем будет состоять подобная мозаика. Затем он захотел научиться описывать все остальные возможные наборы таких фигур.

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

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

Эту же задачу можно рассматривать и в трёх измерениях. Лет 50 назад русский физик-ядерщик, лауреат нобелевской премии мира, и позднее диссидент, Андрей Дмитриевич Сахаров задумался над такой же проблемой, когда со своей женой резал кочаны капусты. Сколько вершин в среднем будет у каждого из полученных кусочков? Сахаров передал эту задачу легендарному советскому математику Владимиру Игоревичу Арнольду и его студенту. Однако полного решения они не нашли, и их попытки по большей части были забыты.


Валуны Моераки в Новой Зеландии

Домокош, не знавший об их работе, написал доказательство, ответом которого стали кубы. Но он захотел проверить его правильность. Он решил, что если ответ к этой задаче уже существует, то он должен быть спрятан в непостижимом труде немецких математиков Вольфганга Вайля и Рольфа Шнайдера 80-летнего титана из области геометрии [в оригинале не указано название видимо, имеется в виду книга "Стохастическая и интегральная геометрия" / прим. пер.]. Домокош профессиональный математик, но текст книги оказался тяжеловат даже для него.

Я нашёл человека, согласившегося прочесть мне нужную часть книги, и перевести её обратно на человеческий, сказал Домокош. Он нашёл там теорему для любого количества измерений. Она подтвердила, что в трёх измерениях в ответе действительно фигурируют кубы.

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

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


Геометрия Марса. Для анализа поверхности в данном случае, поверхности кратера Марса, похожей на соты исследователи размечают все вершины и стороны. Они подсчитывают количество вершин для каждой из ячеек и количество ячеек, для которых каждая из вершин является общей.

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

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

Опять-таки, этот диапазон достаточно легко найти в двух измерениях, но гораздо труднее в трёх. В трёхмерном пространстве кубы очень хорошо складываются вместе, но есть и другие виды фигур в том числе, формирующие трёхмерные версии диаграммы Вороного. Чтобы не переусложнять задачу, Домокош ограничился мозаикой из правильных выпуклых ячеек с общими вершинами. В итоге они с математиком Жолтом Ланги вывели новую гипотезу, набросав кривую, в которую укладываются все возможные трёхмерные мозаики. Они опубликовали работу в журнале Experimental Mathematics, а потом я отправил всё это Рольфу Шнайдеру, нашему божеству, сказал Домокош.


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


Я спросил его, нужно ли объяснять, как я пришёл к такой гипотезе, но он сказал, что он знает об этом, смеётся Домокош. Для меня это было в сто раз более важным, чем принятие статьи любым журналом мира.

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

Но чтобы всё это имело значение для кого-то, кроме небольшой группки математиков, Домокош должен был доказать, что этим правилам подчиняется и реальный мир.

От геометрии к геологии


К тому времени, как в 2016 году Домокош оказался в Филадельфии, в решении задачи применительно к реальному миру он уже кое-чего достиг. Они с коллегами из Будапештского технологическо-экономического университета собрали осколки доломита, отколовшиеся от скалы Хармашатархег, находящейся в Будапеште. Несколько дней сотрудник лаборатории без всякого предубеждения насчёт кубов усердно подсчитывал количество граней и вершин сотен кусочков. Какие средние показатели он получил? Шесть граней, восемь вершин. Домокош совместно с Яношом Тёроком, специалистом по компьютерным симуляциям, и Ференцем Куном, экспертом по фрагментационной физике, обнаружили, что средние кубоиды появлялись и в породах другого типа например, в гипсе и известняке.

Вооружившись математикой и ранними физическими свидетельствами, Домокош выдал свою идею поражённому Джерольмаку. Он будто загипнотизировал меня, и всё остальное на какое-то время просто исчезло, сказал Джерольмак.

Их альянс был не нов. Много лет назад Домокош обрёл известность, доказав существование гёмбёца забавной трёхмерной фигуры, упорно переворачивающейся в определённую позицию равновесия. Чтобы узнать, может ли гёмбёц существовать в реальности, он привлёк Джерольмака, помогшего применить эту концепцию для объяснения круглой формы гальки на Земле и Марсе [Владимир Арнольд приложил руку и тут, впервые поставив вопрос о существовании подобных тел / прим. пер.]. Теперь Домокош снова просил помочь превратить некие теоретические математические концепции в осязаемый камень.


Гёмбёц выпуклая трёхмерная однородная фигура, имеющая ровно одну точку устойчивого равновесия и одну неустойчивого

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

Сначала всё вроде бы работало, сказал Джерольмак. Математика Домокоша предсказывала, что осколки камней в среднем должны быть кубами. Всё большее число реальных осколков вроде как укладывалось в эту теорию. Однако вскоре Джерольмак понял, что для доказательства теории необходимо разобраться и с исключениями из правил.

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

Возможно, эти примеры полностью опровергли бы теорию кубического мира. А возможно, что было бы интереснее, они появлялись бы только в особых случаях, из которых геологи смогли бы извлечь новые уроки. Я сказал, что знаю, что это работает не везде, и мне нужно знать, почему, сказал Джерольмак.

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


Вселенная плиток. Все возможные выпуклые плитки, полностью закрывающие плоскость, можно нанести на график соответствия среднего числа вершин у плитки (ось y) и среднего количества клеток, делящих одну вершину (ось х). Примеры из реального мира:
6 мостовая гигантов, 7 вечная мерзлота на Аляске, 8 высохшая грязь, 9 поверхность гранита.


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

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


Диаграмма Вороного делит плоскость на отдельные участки, каждый из которых состоит из всех точек, ближе всего расположенных к начальной.

Подсчёт корки


В процессе разработки команда встретилась в Будапеште и три безумных дня лихорадочно старалась включить в модель больше примеров из реальной жизни. Вскоре Джерольмак вывел на экран компьютера новый узор: мозаику тектонических плит Земли. Плиты располагаются на литосфере, почти двумерной коже на поверхности планеты. Узор выглядел знакомым, и Джерольмак подозвал других, чтобы полюбоваться им. Мы все были потрясены, сказал он.

На первый взгляд кажется, что плоскостные рисунки тяготеют к диаграмме Вороного, а не к квадратной сетки. А потом команда произвела подсчёты. В идеальной мозаике Вороного из шестиугольников на плоскости у каждой ячейки должно быть по шесть вершин. У реальных тектонических плит в среднем оказалось по 5,77 вершин.

На этом месте гефизик уже мог бы праздновать победу. Но математика это не устраивало. Настроение Дага поднималось. Он работал, как заведённый, сказал Домокош. А на следующий день я расстроился, поскольку размышлял об этом разрыве.

Вечером Домокош отправился домой, всё ещё снедаемый этой разницей. Он вновь записал все числа. И вдруг на него снизошло откровение. Мозаика из шестиугольников может замостить плоскость. Но Земля не плоская по крайней мере, за пределами некоторых одиозных уголков на YouTube. Представьте себе футбольный мяч, состоящий из пятиугольников и шестиугольников. Домокош обработал данные с учётом сферической поверхности и обнаружил, что на шаре у ячеек мозаики Вороного в среднем должно быть по 5,77 вершин.

Эта идея помогла исследователям решить один из важных и открытых вопросов в геофизике: как формируются тектонические плиты Земли? Одни считают, что эти плиты побочный продукт движущихся глубоко в мантии конвекционных потоков. Их противники полагают, что земная кора является отдельной системой. Она расширялась, стала хрупкой и сломалась. Соответствие плит диаграмме Вороного, напоминающей грязевую корку, может поддержать вторую теорию, сказал Джерольмак. А ещё это дало мне понять, насколько важной была та работа, сказал Аттал. Феноменально.

Переломный момент


В трёх измерениях же исключений из правила кубов встречалось достаточно мало. И их тоже можно было объяснить при помощи симуляции необычных сил, давящих наружу. Одна явно некубическая формация находится на побережье Северной Ирландии, где волны бьются о десятки тысяч базальтовых колонн. На ирландском языке это называется Clochn na bhFomhrach, дорога из камней для сверхъестественных существ. По-английски это называется "мостовой гигантов".

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


Мостовая гигантов в Северной Ирландии

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

Они должным образом проверили свою модель на соответствие реальности, сказала Марта-Кэри Эппс, специалист по естественным наукам из университета Северной Каролины. Мой изначальный скептицизм угас.

Математика говорит нам, что если мы будем дробить камни, как угодно, случайно или специально, у нас всё равно окажется ограниченный набор возможностей, сказал Фёрбиш. Разве это не умно?

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

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

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

Это буквально самый очевидный из примеров, которые можно придумать. Статистическое среднее всех этих наблюдений представляет собой куб, сказал Джерольмак. Но такого куба не найти".
Подробнее..

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

21.12.2020 14:15:23 | Автор: admin

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


Рений (Re) 75-й элемент таблицы Менделеева. Был открыт в 1925 году Вальтером и Идой Ноддак, является хронологически последним открытым элементом, имеющим стабильные изотопы. Рений является чрезвычайно рассеянным редким металлом, содержится в основном в медно-молибденовых рудах. Практическая ценность и даже незаменимость рения заключается в его крайней жаропрочности. Температура плавления рения составляет 3186 градусов Цельсия, поэтому применение рения для легирования стали имеет принципиальное значение при производстве деталей для реактивных двигателей и для аэрокосмической промышленности в целом. Так, рений добавляют в лопатки турбин реактивных двигателей, что позволяет повысить их рабочую температуру на 50-80 градусов. Сплавы рения с вольфрамом находят применение в электронике, а сплавы рения с платиной в нефтехимической промышленности, в качестве катализаторов, участвующих в удалении серы. Феноменальная прочность рения сегодня изучается, в частности, в Национальном исследовательском технологическом университете МИСиС; подробнее об этом можно почитать в блоге института на Хабре в замечательной статье Материал, который выжил. В качестве небольшого спойлера отметим, что речь идет о модификации рения при помощи азота, а именно о пернитриде рения. Показатели статьи:

Основные промышленно разрабатываемые залежи рения находятся в Чили, а также в Казахстане, Узбекистане и Армении. Разведанные запасы металла также имеются в США. Объяснимый интерес США к приобретению рения для развития программы Шаттл привел к тесному сотрудничеству США и Казахстана в этой сфере в 1990-е годы. В тот период Российская Федерация, также крайне нуждающаяся в 75 элементе, была вынуждена извлекать его из промышленных отходов. Но еще в 1992 году начались подвижнические исследования Генриха Семеновича Штейнберга и Феликса Исааковича Шадермана, открывших уникальное месторождение рения на фумарольных полях вулкана Кудрявый.

Подробный и очень увлекательный рассказ об этих исследованиях вы найдете в документальном фильме "Эффект рения" (2014), повествующем об очередной экспедиции на Итуруп. Фильм содержит комментарии Генриха Семеновича, а также других специалистов, завораживающие пейзажи российского Дальнего Востока и сцены героической, по-настоящему адовой работы геологов.

История проблемы

Генрих Семенович Штейнберг еще в начале 90-х был осведомлен о необычайно разнообразном минеральном составе фумарол на Южных Курилах. Вулканы выбрасывают в атмосферу дым, насыщенный индием, золотом, германием и другими редкими элементами. В поисках образцов рения Штейнберг обратил внимание на Кудрявый, так как дым в его фумаролах имеет необычно высокую температуру: на Главном поле она достигает 870 градусов Цельсия, а на куполе вулкана около 750 градусов. Геохимические исследования показали не только высокое содержание рения в фумаролах, но и позволили открыть ранее не известный минерал, сульфид рения ReS2, встречающийся на склонах вулкана и на его фумарольных полях. Минерал был назван рениитом.

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

Остров Итуруп сравнительно невелик, и возведение горно-обогатительного комбината на территории острова, а также утилизация пустой породы в настоящее время не представляются осуществимыми. Кроме того, вулкан Кудрявый отличается значительной геологической активностью. Последнее серьезное извержение зафиксировано в 1883 году, но в 1946 и 1999 годах наблюдались фреатические взрывы. Более того, в середине 2019 года прогнозировалось близкое полноценное извержение Кудрявого, так как открытые минералы на его склонах раскалились до 1600 градусов по Цельсию.

В вышеупомянутом фильме Эффект рения констатируется, что в настоящее время стоит сосредоточиться на извлечении минерала из корки, то есть, из рениита. Наряду с сульфидом рения в районе Кудрявого также встречается соль перринат аммония, содержание рения в которой составляет порядка 70%.

Как можно добывать фумарольный рений

Пока не существует промышленной технологии изготовления ловушек, которые позволили бы стабильно собирать фумарольный дым и извлекать из него рений. Основные сложности связаны с быстрой изнашиваемостью конструкций во время эксплуатации при таких температурах (возможно, конечно, помогло бы легирование рением). Тем не менее, опытный образец ловушки для рения был сконструирован под руководством Штейнберга еще в 2015 году и запатентован. В ходе испытаний устройство позволило собрать 17 килограммов концентрата рения. Но лишь в апреле 2020 года появилась новость о том, что ВНИИ химической технологии (ВНИИХТ) и Институт вулканологии и геодинамики (ИВиГ) разработали оборудование по улавливанию рения из вулканических газовых испарений. Конструкция представляет собой систему куполов, которые предполагается развернуть прямо в кратере вулкана и при помощи катализаторов извлекать из паров рений.

При этом еще в 1999 году началось развитие технологий, связанных с пропусканием рения через фильтрующий слой. Требуется вещество, которое могло бы впитывать рениевые пары как губка. В таком качестве могут использоваться, в частности, минеральная вата, активированный уголь и, в особенности, алюминиесодержащие минералы из группы цеолитов. Необходимо отметить, что крупное месторождение цеолитов (Лютогорское) имеется на Сахалине. При использовании такой технологии целесообразно охлаждать обрабатываемый газ до температуры около 150 градусов путем разбрызгивания воды в дымососе. В таком случае из газа выделяется элементная сера, которая служит дополнительным уловителем рения.

Наряду с абсорбционными технологиями улавливания рения разрабатываются электрохимические; правда, пока электрохимические технологии далеки от практического применения ввиду дороговизны оборудования. Тем не менее, активное исследование абсорбционных и электрохимических технологий на Кудрявом продолжается и в 2020 году. Специализированная экспедиция была организована силами Томского государственного университета и специалистов РАН.

В любом случае, после 2014 года в описанном проекте достигнуты серьезные успехи, а сам проект промышленной добычи рения, индия и германия на Кудрявом решается на государственном уровне. Согласно государственной программе Социально-экономическое развитие Курильских островов (Сахалинская область) на 2016-2025 годы, строительство инфраструктуры для проекта предполагается провести в 2020-2023 годах. Правда, из-за труднодоступности Кудрявого работа предприятия в настоящее время возможна только в летний период.

Кудрявый продолжает покуривать свою драгоценную магму и никуда не спешит. В том числе, не спешит делиться с нами своими сокровищами. Но перспективы этого уникального российского месторождения кажутся оптимистичными.

Кстати, Генриху Семеновичу Штейнбергу 13 февраля 2021 года исполнится 86. Надеюсь, что он успеет лично оценить самые интересные события на Кудрявом и плоды, которые принесет труд всей его жизни. В этом была бы экзистенциальная и подлинная справедливость.

Подробнее..

Подвижный элемент

23.03.2021 08:12:20 | Автор: admin

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

Олимпиадинское месторождение находится на полтысячи километров севернее города Красноярска. Начиная с 80-х годов 20 века на месторождении добыли более пятисот тонн золота и, по расчётам специалистов, еще около тысячи тонн осталось в недрах.
Первичные руды Олимпиады это плотные породы с вкраплением сульфидных минералов (соединений металлов и серы) и очень мелкого самородного золота (рис. 1).


Содержание золота в таких рудах составляет 3-5 граммов на тонну.

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

Максимальное содержание золота в окисленных рудах Олимпиадинского месторождения достигало 450 граммов на тонну. Ранее считалось, что в таких условиях золото ведет себя как инертный металл не подвергается коррозии и окислению, и, в отличие от большинства неблагородных металлов, не реагирует с щелочами и кислотами (кроме царской водки). Однако исследователи СФУ получили доказательства его подвижности и активности золота в гипергенных условиях, при воздействии на него особых геологических и климатических факторов.

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

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

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

Сейчас известно, что подобные процессы могут происходить даже с участием бактерий, которые способны восстанавливать наночастицы металлического золота из его соединений с другими элементами.
Исследователи отмечают, что процесс формирования окисленных руд, а также миграция и новообразование золота Олимпиадинского месторождения хорошо иллюстрируются в разрезе толщи окисленных руд. В нижней части содержится много реликтового золота первичных руд, а также большое количество губчатого золота. Выше по разрезу доля губчатого и реликтового золота снижается, но увеличивается количество новообразованных глобуль и микрокристаллов. А вот наибольшее количество глобуль и микрокристаллов, как и максимальное суммарное содержание золота (до 60 грамм на тонну в изученных пробах), наблюдается в верхнем слое изученного разреза. Эта зона примечательна еще тем, что здесь найден минерал церианит (оксид церия, CeO2), выступающий индикатором резкого изменения условий окружающей среды. Именно с этим изменением условий связано новообразование золота, считают авторы.
Хотя окисленные руды Олимпиадинского месторождения были отработаны еще в 2007 году, изучение форм нахождения, поведения и распределения золота в них важно по ряду причин. В первую очередь, окисленные руды ценны, поскольку не требуют сложных процессов обогащения и рентабельны даже при низких содержаниях золота в них. С другой стороны, зоны окисления золоторудных месторождений представляют и фундаментальный интерес, связанный с экзогенной геохимией и металлогенией золота, отметил ведущий инженер кафедры геологии, минералогии и петрографии СФУ Сергей Сильянов.
Соавтор исследования, инженер аналитической лаборатории Научно-технологического R&D центра ГМК Норильский никель Борис Лобастов, в свою очередь сообщил, что зона окисления исследованного месторождения была самой настоящей лабораторией под открытым небом.
Нами было изучено немало техногенных объектов там процессы переноса вещества и образования новых минералов происходят особенно быстро. Удивительно было увидеть, насколько похожи микрокристаллы золота и многих других минералов из окисленных руд Олимпиадинского месторождения на новообразованные минералы других объектов. Такое внешнее сходство не всегда списывается на конвергенцию механизмы формирования новых минералов в столь разных обстановках схожи, а изучение одних процессов в данном случае помогает понять все остальные, сказал учёный.
Выполненное сибирскими учёными исследование показало, что формирование окисленных руд Олимпиадинского месторождения происходило при сложных геохимических условиях, в которых золото вело себя как подвижный элемент, способный и к перераспределению, и к переотложению на геохимических барьерах (участках изменения условий окружающей среды).
Источник
Подробнее..

Перевод Реальная история легендарной денежной ямы Острова Оук

05.06.2021 18:10:05 | Автор: admin

Золото тамплиеров, пиратское хранилище, природная карстовая воронка или один гигантский обман? Проклятие острова Оук рассказывает нам, что семь человек должны умирать до того, как остров раскроет своё легендарное сокровище. Шесть человек погибли в поисках миллиардов в золоте, но опасность только подогревает исследования и спекуляции.

Разочаровывающий, увлекательный, манящий вы можете поставить любое прилагательное перед словами остров Оук, и будете правы, рассказывает Чарльз Баркхаус, историк шоу канала History Channel The Curse of Oak Island, в котором на протяжении восьми сезонов (с некоторыми результатами) рассказывается о продолжающихся поисках сокровищ.


Если вы не готовы к эмоциональной нагрузке, то можете собирать свои игрушки и уходить.

Охота за сокровищами денежной ямой острова Оук, то есть яма в 100 футов (около 30 с половиной метров) на острове в Новой Шотландии, где якобы хранится что угодно от пиратских сокровищ до Ковчега Завета, началась ещё в 1795 году. Хотя сокровища так и не были найдены, сопутствующие открытия очевидные подсказки, возможные ловушки и геологические диковинки заставляют искателей продолжать поиски, даже когда историки оспаривают более сенсационные заявления, связанные с кладом. Был ли остров Оук сокровищницей рыцарей-тамплиеров, секретным промышленным центром Британии или злополучной природной воронкой? Чтобы ответить на этот вопрос, конечно же, нужно копать.

Поиски начинаются

Остров Оук впервые привлёк к себе внимание вскоре после золотого века пиратства (примерно в 16501730 годах), когда Эдвард Лоу и Бартоло Мью Робертс патрулировали моря к северо-востоку от Америки. В 1795 году подросток из Новой Шотландии со своего дома на материке увидел парящие над островом странные огни.

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

Мальчики начали выкапывать то, что впоследствии стали называть денежной ямой. На глубине двух футов (60 сантиметров) они обнаружили круг из камней, окаймляющих окружность ямы, а на глубине 10 футов (3 метров) платформу из подогнанных в стенки ямы обрезков брёвен. Вторая платформа лежала на 20 футов (около 6 метров) ниже, но на этом рассказ о первом поиске заканчивается.

История возобновляется в начале 1800-х годов, когда компания Онслоу отправилась в первую официальную экспедицию для раскопок. Они продолжили раскопки с того места, где остановились в первый раз, обнаружив новые платформы через каждые 10 футов (около 3 метров) (приблизительно три метра), иногда со слоями замазки, древесного угля или кокосовых волокон. Кокосы не растут в радиусе 900 миль (приблизительно полтора километра) от Новой Шотландии, но в истории утверждается, что экипаж сделал ещё более грандиозное открытие на высоте 90 футов (более 27 метров): прямоугольный камень, исписанный странными знаками.

Остров Оук. Пиковая высота острова Оук составляет всего 36 футов (почти 11 метров) над уровнем моря. Архивы Новой ШотландииОстров Оук. Пиковая высота острова Оук составляет всего 36 футов (почти 11 метров) над уровнем моря. Архивы Новой Шотландии

Исследователи и охотники за сокровищами сочли, что эти отметки были сделаны по ошибке инструментами экскаватора, но другие были уверены, что это секретный код, ведущий к зарытым сокровищам. В 1860-х годах профессор лингвистики из университета Далхаузи в Новой Шотландии изучил камень и определил, что код представляет собой подстановочный шифр: Сорок футов [кооло 12 метров] ниже погребены два миллиона фунтов. Но другая попытка перевода в 1970-х годах интерпретировала код как христианское предупреждение коптов не забывать о своём долге перед Господом.

Компания Онслоу продолжала копать, и на глубине 98 футов (кооло 30 метров) они обнаружили нечто, по звуку удара напоминавшее полый контейнер предположительно, хранилище сокровищ. Бригада прекратила работу на вечер, но когда они вернулись на следующее утро, то обнаружили, что котлован заполнен водой на 60 футов (около 18 метров). Предполагалось, что раскопки привели к срабатыванию мины-ловушки. И похоже, что наводнение, похоже, положило конец усилиям Онслоу; в 1805 году компания была распущена.

Удар проклятия

Ещё одна экспедиция на Оук была запущена из близлежащего города Труро в 1849 году. Команда могла откачать воду из ямы и укрепить стены ямы перед бурением в хранилище. Сверло проникало последовательно расположенные слои древесины и свободного металла и это предполагает сундук с сокровищами, а, согласно опубликованной годы спустя новости в газете, всплыли три небольших звена золотой цепочки. Но до того, как команда смогла получить доступ к хранилищу, нижняя часть денежной ямы рухнула и снова была затоплена, унеся за собой предполагаемое сокровище.

Ничуть не смутившись, люди подумали, что нашли затопленный тоннель, соединявший яму и рукотворную бухту Смита примерно на 500 футов (около 150 метров) к востоку от площадки, где они копали.

Остров Оук. К моменту съёмки в 1947 году охота унесла две жизни. Архивы Новой Шотландии.Остров Оук. К моменту съёмки в 1947 году охота унесла две жизни. Архивы Новой Шотландии.

По мере распространения новостей об опасной, драматической охоте, на остров Оук прибывало всё больше экскурсий, которые открывали все новые ямы, но это лишь скрывало истину о сокровищах. В 1897 году поисковики обнаружили крошечный кусочек пергамента с буквами vi.

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

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

Не дали прорыва и современные технологии. Съёмочная группа 1971 года привезла камеры и монитор, чтобы исследовать примерно 235-футовую (чуть менее 72 метров) шахту (называемую скважиной 10X и находящуюся примерно в 180 футах (около 55 метров) от денежной ямы), и хотя они утверждали, что во время исследования на дне ямы был обнаружен деревянный сундук и отрубленная человеческая рука, этот инцидент так и не был зарегистрирован.

Рик и Марти Лагина, входящие в группу владельцев острова Оук, ведут сериал на канале History Channel.Рик и Марти Лагина, входящие в группу владельцев острова Оук, ведут сериал на канале History Channel.

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

Несмотря на это, рассказывает Баркхаус, количество доказательств, которые мы находим [на острове Оук], заставляет вас поверить, что какая-то группа или какой-то человек сделали там нечто особенное.

Дыры в поиске

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

Есть другие пробелы: многие считают, что звенья золотой цепи от группы 1849 года были подброшены самой группой для поощрения будущих экспедиций, а камень с надписью, найденный в начале 1800-х годов, был зарегистрирован как найденный только в 1862 году. Этот камень вообще не упоминался в инвестиционном проспекте компании Oak Island Treasure Company за 1893 год: ни сам камень, ни его маркировка не были зарисованы или сфотографированы, а существующее сегодня изображение камня датируется 1949 годом. Именно с этого времени начинаются современные переводы.

Остров Оук. Ведущая на материк дамба позволяет искателям сокровищ доставлять на остров Оук тяжёлую технику. Архивы Новой ШотландииОстров Оук. Ведущая на материк дамба позволяет искателям сокровищ доставлять на остров Оук тяжёлую технику. Архивы Новой Шотландии

Можно привести аргумент существует общая структура [историй] о зарытых сокровищах, рассказывает Даунс. Сокровища каким-то образом теряются, истории о них рассказываются так, будто они правда, а тот факт, что сокровища так и не найдены подпадает под категорию подтверждающих формул, добавляемых в легенду, чтобы она стала достовернее.

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

Небылицы о денежной яме

Существует множество спекуляций по поводу того, что именно представляет собой денежная яма острова Оук: общий банк для пиратов, хранилище сокровищ, награбленных британскими военными, или сберегательный счёт для финансирования Американской революции в натуральном виде.

Согласно одной из твёрдых теорий, яма была построена рыцарями ордена тамплиерами, предшественниками масонов. Исследователи интерпретировали некоторые гравюры на острове, а также большой и неясный масонский каменный треугольник как свидетельство того, что остров был местом, где прятались религиозные артефакты, такие как Святой Грааль и Ковчег Завета.

Даунс считает,эти теории движимы теми же механизмами, что и и вера в теории заговора. Людям нравится верить, что в окружающем нас мире существует некий порядок, считает она.

Реальность такова, что они не всегда [достаточно знают об ордене]. Чтобы придать смысл миру, который нам трудно понять, мы часто придумываем повествования.

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

Сцена из сериала Проклятие острова Оук. Часто люди, предполагающие существование сокровищ на острове Оук, имеют долю его землиСцена из сериала Проклятие острова Оук. Часто люди, предполагающие существование сокровищ на острове Оук, имеют долю его земли

Тем не менее правда об острове Оук сохраняет элементы интриги. Исследование, проведённое историком Джой А. Стил и отставной морской геолог Гордон Фейдер доказывают, что на острове Оук располагался секретный британский промышленный центр.

Изучив деловую документацию и современную переписку, пара пришла к выводу, что в 1720 году корона совместно с британскими военными зафрахтовала частные компании, чтобы вести дела на острове Оук, включая производство сосновой смолы, изготовление латуни и волочение проволоки, чтобы помочь погасить задолженность. В то время это была крупнейшая промышленная разработка в Канаде, рассказывает Фейдер: Был миллион причин отправиться на остров Оук он ближе всего к пресной воде, ближе всего к берегу, безопасен; остров самый большой в заливе это хорошая стоянка.

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

По словам Фейдера, В те времена сосновая смола имела такое же значение, как сегодня нефть: Корабль не выходил в море, если не был пропитан сосновой смолой. Именно этим люди и занимались на острове. Все артефакты, которые мы видим, точно соответствуют этой теории.

Канцлер казначейства Англии (по сути, министр финансов) и другие высокопоставленные банковские чиновники того времени часто ссылались на Секрет в своей переписке, рассказывает Стил, который [несомненно] проект Оук-Айленд.

Я удивлена, что историки никогда не обращали на это внимания и не осмеливались выяснить, в чём заключался секрет, говорит она. [Остров Оук] воплотил то, что могло бы стать очень прибыльной товарной схемой.

Контрсвидетельства природы

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

Фейдер составил карту подземных условий вдоль побережья Новой Шотландии, работая морским геологом в канадском правительстве, и он рассказывает, что карстовые воронки в этом районе обычное явление.

Можно наблюдать, как развивается воронка и как сосна в 60 футов (примерно 18 метров) падает [в неё] за две секунды, говорит он, иллюстрируя, как в яме могли разместиться предметы размером с мачты британских кораблей.

На острове Оук обнаружено по меньшей мере две карстовые воронки, а на прилегающем материке их гораздо больше. Стивен Эйткен, доктор наук, геофизик с более чем 25-летним опытом изучения района острова Оук, считает: природные данные указывают на то, что Денежная яма сама по себе является карстовой воронкой. В соответствии с утверждениями Фейдера Эйткен рассказывает, что горная порода под этой стороной острова локально растворилась, образовав карстовую систему, и некоторые пещеры этой системы обрушились и образовали карстовые воронки, включая Денежную яму.

Карстовые ямы это мусорные баки геологического мира, рассказывает Эйткен. Они часто заполнены обрушившейся брекчией [разновидностью осадочной породы], перекрытой смесью органических обломков и отложений.

Эйткен пишет, что отложения над брекчией в Денежной яме утолщены до 37 футов (чуть более 11 метров), а это указывает на позднее заполнение впадины. По словам Эйткена, первоначальная впадина, обнаруженная в 1795 году, также соответствует типичному проявлению карстовой воронки на поверхности.

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

Идея о том, что пираты вырыли сокровищницу кирками в скальных породах, просто смешна, также добавляет он.

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

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

Куда бы вы ни шли по этому острову, вы идёте сквозь историю, утверждает Баркхаус. Нельзя взять и свалить всё найденное на одну теорию; вот что так странно на острове. Я никогда не перестану верить, что на острове Оук зарыт клад, иэтот клад всё ещё там. Я чувствую это.

Притягательность настоящих сокровищ, будь то Святой Грааль, Ковчег Завета или огромный пиратский банк, затмевает собой все остальные находки на острове Оук. Охота может никогда не прекратиться. Навязчивая идея богатства, как указал один из зарегистрированных искателей сокровищ ещё в 1862 году, может стоить борьбы.

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

Как формируются карстовые воронки?

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

В зависимости от геологических условий и вышележащих материалов карстовые воронки могут образовываться в течение десятилетий и разрушаться в считаные секунды. Обрушение может быть ускорено внезапным притоком воды в недра или замерзанием и оттаиванием. Фейдер рассказывает, что при благоприятных обстоятельствах карстовая воронка может быть от 130 до 165 футов (чуть более 50 метров) в поперечнике. Он уверен, что именно этот феномен истинная причина происхождения денежной ямы острова Оук.

Если вам интересны не только исследования и поиски сокровищ и вы хотите исследовать данные и извлекать пользу из них, то вы можете присмотреться к нашему курсу Data Science, итог которого эквивалентен двум-трём годам активных самостоятельных поисков в науке о данных.

Узнайте, как прокачаться и в других специальностях или освоить их с нуля:

Другие профессии и курсы
Подробнее..

AB-тест в инженерно-геологических изысканиях на языке Python

05.06.2021 14:06:41 | Автор: admin

1. Введение

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

При данной постановке задачи можно применить методику A/B-тестирования со следующими параметрами:

  1. Измеряемой метрикой будет среднее значение плотности скелета грунта (pd, г/см3), характеризующее сложение проб. Данная величина имеет нормальный закон распределения;

  2. Критерием проверки гипотезы будет служить t-критерий (критерий Стьюдента):для двух независимых выборок, если сопоставляемые полевые (до транспортировки) и лабораторные (после транспортировки) данные проводились на разных пробах грунта;для двух зависимых выборок, если исследования выполнены на одних и тех же пробах.

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

2. Генерация выборок

2.1 Оценка объема выборок

В рамках дизайна эксперимента, перед генерацией выборок плотностей, прикинем их необходимый объем при заданномразмере эффекта (ES - effect size),мощности (power)идопустимой ошибке I рода ()(определения данных терминов приведено ниже). Расчет произведем с привлечением пакетаstatsmodels.

Размер эффекта (стандартизированный) величина, характеризующая различие, которое мы хотим выявить, равная отношению разности средних значений по выборкам к взвешенному стандартному отклонению. В нашем случае:

{ES = \frac{{(\bar{X}_1 - \bar{X}_2)}_{obs}}{{S}_{pooled}}}{ES = \frac{{(\bar{X}_1 - \bar{X}_2)}_{obs}}{{S}_{pooled}}}

Взвешенное стандартное отклонениеSpooledдля выборок одинакового размера можно расcчитать по формуле:

{S}_{pooled}= {\sqrt{\frac{{S}_{1}^2+{S}_{2}^2}{2}}}

Существует условная классификация размера эффекта (Cohen, 1988) ES = 0.2 - маленький; 0.5 - средний; 0.8 - большой.

Мощность вероятность не совершить ошибку II рода (обычно принимается равной 80%).

Пояснения по ошибкам I и II рода приведены в таблице ниже:

H0верна

H1верна

H0принимается

H0верно принята

Ошибка II рода ()

H0отвергается

Ошибка I рода ()

H0верно отвергнута (power = 1-)

Для описанных выше величин примем следующие значения:

  • = 0.05 (вероятность выявить различия между средними при их отсутствии)

  • ES = 0.5 (размер эффекта составит половину от дисперсии измеряемых величин плотности).

  • Power = 0.8 (вероятность выявления установленного различия между средними значениями).

Теперь к коду:

#Импорт библиотекimport numpy as npfrom statsmodels.stats.power import TTestIndPowerfrom matplotlib.pyplot import figureimport matplotlib.pyplot as pltimport scipyfrom statsmodels.stats.weightstats import *
#Задаем параметрыeffect = 0.5alpha = 0.05power = 0.8analysis = TTestIndPower()#Оценка размера выборкиsize = analysis.solve_power(effect, power=power, alpha=alpha)print(f'Размер выборки, шт.: {int(size)}')

Размер выборки, шт.: 63

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

Давайте построим график зависимости необходимого размера выборок от размера эффекта при заданной мощности и уровне значимости.

plt.figure(figsize=(10, 7), dpi=80)results = dict((i/10, analysis.solve_power(i/10, power=power, alpha=alpha))                for i in range(2, 16, 1))plt.plot(list(results.keys()), list(results.values()), 'bo-')plt.grid()plt.title('График зависимости необходимого объема выборки \n от размера эффекта')plt.ylabel('Размер выборки n, шт.')plt.xlabel('Размер эффекта ES, д.е.')for x,y in zip(list(results.keys()),list(results.values())):    label = "{:.0f}".format(y)    plt.annotate(label,                  (x,y),                  textcoords="offset points",                  xytext=(0,10),                  ha='center')plt.show()

Данный график позволяет увидеть, как быстро изменяется необходимый объем выборок при уменьшении фиксируемого размера эффекта ES. Например: при выявлении различия в плотности проб грунта до и после их транспортировки в 0,03г/см3при стандартном отклонении в 0,1г/cм3(ES = 0,03г/см3/ 0,1г/см3= 0,3 д.е.), необходимый объем проб по каждой выборке должен составить не менее 175 проб для заданной мощности и уровня значимости (power=0.80,=0.05).

2.2 Генерация выборок

Теперь зная необходимый минимальный размер выборок, сгенерируем их с помощью библиотекиnumpy.

Измеряемая физическая характеристика грунта (плотность скелета) имеет нормальный закон распределения. В рамках данного примера зададим генератору следующие значения среднего (X) и стандартного отклонения (S):

  • для первой выборки X1= 1,65г/см3,S1= 0.15г/см3;

  • для второй X2= 1,60г/см3,S2= 0.15г/см3.

loc_1 = 1.65sigma_1 = 0.15loc_2 = 1.60sigma_2 = 0.15sample_size = 65#Генерируем выборки с заданными параметрамиsample_1 = np.random.normal(loc=loc_1, scale=sigma_1, size=sample_size)sample_2 = np.random.normal(loc=loc_2, scale=sigma_2, size=sample_size)

Постоим гистограммы и "ящик с усами" по полученным выборкам.

fig, axes = plt.subplots(ncols=2, figsize=(18, 5))max_y = np.max(np.hstack([sample_1,sample_2]))#Гистрограмма по выборке 1count_1, bins_1, ignored_1 = axes[0].hist(sample_1, 10, density=True,                                           label="Выборка 1", edgecolor='black',                                          linewidth=1.2)axes[0].plot(bins_1, 1/(sigma_1 * np.sqrt(2 * np.pi)) *               np.exp( - (bins_1 - loc_1)2 / (2 * sigma_12)),         linewidth=2, color='r', label='плотность вероятности')axes[0].legend()axes[0].set_xlabel(u'Длина сессии, с')axes[0].set_ylabel(u'Количество сессий, шт.')axes[0].set_ylim([0, 5])axes[0].set_xlim([1.1, 2.2])#Гистрограмма по выборке 2count_2, bins_2, ignored_2 = axes[1].hist(sample_2, 10, density=True,                                           label="Выборка 2", edgecolor='black',                                           linewidth=1.2, color="green")axes[1].plot(bins_2, 1/(sigma_2 * np.sqrt(2 * np.pi)) *               np.exp( - (bins_2 - loc_2)2 / (2 * sigma_22)),         linewidth=2, color='r', label='плотность вероятности')axes[1].legend()axes[1].set_xlabel(u'Длина сессии, с')axes[1].set_ylabel(u'Количество сессий, шт.')axes[1].set_ylim([0, 5])axes[1].set_xlim([1.1, 2.2])plt.show()
#Ящик с усамиfig, ax = plt.subplots(figsize=(8, 8))axis = ax.boxplot([sample_1, sample_2], labels=['Выборка 1', 'Выборка 2'])data = np.array([sample_1, sample_2])means = np.mean(data, axis = 1)stds = np.std(data, axis = 1)for i, line in enumerate(axis['medians']):    x, y = line.get_xydata()[1]    text = ' ={:.2f}\n ={:.2f}'.format(means[i], stds[i])    ax.annotate(text, xy=(x, y))plt.ylabel('Плотность скелета грунта, г/см3')plt.show()

3. Формулировка гипотез

Пришло время для формулировки гипотез. У нас могут быть два случая:

  • Случай 1. Сопоставляемые полевые и лабораторные данные по определению плотности скелета грунта относятся к разным пробам, тогда t-критерий будет рассчитываться для двух независимых выборок;

  • Случай 2. Исследования в поле и лаборатории выполнены на одних и тех же пробах, тогда t-критерий будет рассчитываться для двух зависимых выборок.

Начнем с первого варианта.

Вариант 1. Для двух независимых выборок

С помощью двухвыборочного критерия Стьюдента проверим гипотезу о равенстве средних выборок.

Нулевая гипотезаH0:средние значения равны1=2.

Альтернативная гипотезаH1:средние не равны12.

Статистика:

T({{X_1}^{n_1}},{{X_2}^{n_2}}) = \frac{\bar{X_1}-\bar{X_2}} {\sqrt{\frac{S_1^2}{n_1} + \frac{S_2^2}{n_2}}}

Нулевое распределение:T(X1n1,X2n2)~St(), где степень свободывычисляется по следующей формуле

{\nu = \frac{ ({\frac{S_1^2}{n_1} + \frac{S_2^2}{n_2}})^2 } {\frac{S_1^4}{n_1^2(n_1-1)}+ \frac{S_2^4}{n_2^2(n_2-1)} } }

Для расчета достигаемого уровня значимости воспользуемся методомttest_indмодуляstats.

t_st, p_val = scipy.stats.ttest_ind(sample_1, sample_2, equal_var = False)print(f't-критерий составил {round(t_st, 2)}')print(f'Рассчитанный t-критерий дает достигаемый \уровень значимости (p-value) равный {round(p_val, 3)}')

t-критерий составил 2.92

Рассчитанный t-критерий дает достигаемый уровень значимости (p-value) равный 0.004

Вывод для варианта 1

Нулевая гипотезаH0о том, что средняя плотность скелета грунта не изменилась после транспортировки,отвергаетсяна уровне значимости 0,05 (достигаемый уровень значимостиp-valueдля сгенерированных выборок составил 0.004) в пользу альтернативной.

Давайте интервально оценим разность средних по данным выборкам.

c_m = CompareMeans(DescrStatsW(sample_1), DescrStatsW(sample_2))print("95%% доверительный интервал: \[%.4f, %.4f]" % c_m.tconfint_diff(usevar='unequal'))

95% доверительный интервал: [0.0235, 0.1228]

Так как ноль не попадает в рассматриваемый 95% доверительный интервал, мы можем сделать вывод, что средние значения рассматриваемых выборок отличаются на уровне значимости в 5%.

Вариант 2. Для двух связанных выборок

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

Нулевая гипотезаH0:средние значения равны 1=2.

Альтернативная гипотезаH1:средние не равны12.

Статистика:

T({{X_1}^{n}},{{X_2}^{n}}) = \frac{\bar{X_1}-\bar{X_2}} {\frac{S}{\sqrt{n}}}S^2 = \frac{1}{n-1} \sum_{i=1}^n (D_i - \bar{D})^2, D_i = X_{1i} - X_{2i}

Нулевое распределение: T(X1n, X2n)~St(n-1)

Для расчета достигаемого уровня значимости воспользуемся методомttest_relмодуляstats.

t_st, p_val = stats.ttest_rel(sample_1, sample_2)print(f't-критерий составил {round(t_st, 2)}')print(f'Рассчитанный t-критерий дает достигаемый \уровень значимости (p-value) равный {round(p_val, 3)}')

t-критерий составил 2.79

Рассчитанный t-критерий дает достигаемый уровень значимости (p-value) равный 0.007

Вывод для варианта 2

Нулевая гипотезаH0о том, что средняя плотность скелета грунта не изменилась после транспортировки,отвергаетсяна уровне значимости 0,05 (достигаемый уровень значимостиp-valueдля сгенерированных выборок составил 0.007).

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

print("95%% confidence interval: [%.4f, %.4f]"      % DescrStatsW(sample_1 - sample_2).tconfint_mean())

95% confidence interval: [0.0208, 0.1255]

Так как ноль не попадает в рассматриваемый 95% доверительный интервал, мы можем сделать вывод, что средние значения рассматриваемых выборок отличаются.

5. Итог

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

Подробнее..

Категории

Последние комментарии

  • Имя: Макс
    24.08.2022 | 11:28
    Я разраб в IT компании, работаю на арбитражную команду. Мы работаем с приламы и сайтами, при работе замечаются постоянные баны и лаги. Пацаны посоветовали сервис по анализу исходного кода,https://app Подробнее..
  • Имя: 9055410337
    20.08.2022 | 17:41
    поможем пишите в телеграм Подробнее..
  • Имя: sabbat
    17.08.2022 | 20:42
    Охренеть.. это просто шикарная статья, феноменально круто. Большое спасибо за разбор! Надеюсь как-нибудь с тобой связаться для обсуждений чего-либо) Подробнее..
  • Имя: Мария
    09.08.2022 | 14:44
    Добрый день. Если обладаете такой информацией, то подскажите, пожалуйста, где можно найти много-много материала по Yggdrasil и его уязвимостях для написания диплома? Благодарю. Подробнее..
© 2006-2024, personeltest.ru