Компьютерное моделирование структуры и порового пространства терригенных осадков тема диссертации и автореферата по ВАК РФ 04.00.01, кандидат геолого-минералогических наук Амосова, Ольга Евгеньевна

• Амосова, Ольга Евгеньевна
• кандидат геолого-минералогических наук
• 2000

Оглавление диссертации кандидат геолого-минералогических наук Амосова, Ольга Евгеньевна

ВВЕДЕНИЕ

Глава 1. СТРУКТУРА ТЕРРИГЕННЫХ ГОРНЫХ ПОРОД

1.1. Общие положения

1.2. Размеры зерен кластического осадка

1.3. Форма зерен как элемент структуры кластического осадка

1.4. Взаимное расположение зерен

1.4.1. Упаковка и координационные числа осадка

1.4.2. Ориентировка зерен как важный показатель структуры осадка

Глава 2. МАТЕМАТИЧЕСКИЕ (ГЕОМЕТРИЧЕСКИЕ) ОСНОВЫ

ГРАВИТАЦИОННОЙ СЕДИМЕНТАЦИИ ЧАСТИЦ

2.1. Феноменологическое описание моделируемого явления

2.2. Терминология и краткие комментарии к ней

2.3. Виды элементов траектории спуска

2.4. Система обозначений

2.5. Условия смены видов элементов траектории спуска

2.6. Условия остановки

2.7. Вопросы, связанные с ограниченностью бассейна седиментации

Глава 3. КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ОСАДКА

3.1. Общий алгоритм моделирования кластического осадка

3.2. Алгоритмы отдельных процедур

3.2.1. Первое касание

3.2.2. Спуск по меридиану частицы

3.2.3. Спуск по ложбине между двумя частицами

3.2.4. Условия остановки или выбора вида элемента траектории спуска в момент касания трех частиц

3.2.5. Спуск по ложбине между частицей и гранью контейнера седиментации

3.2.6. Условия остановки или выбора вида элемента траектории движения в момент касания двух частиц и грани контейнера седиментации

Глава 4. РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ ОСАДКА

4.1. Коллекторские и некоторые другие структурные свойства осадка

4.2. Сферическое распределение точек касания

4.3. Координационные числа осадка

Глава 5. СТРУКТУРА ПОРОВОГО ПРОСТРАНСТВА

5.1. Анализ структуры порового пространства

5.2. Характеристики порового пространства правильных упаковок шаров

Введение диссертации (часть автореферата) на тему «Компьютерное моделирование структуры и порового пространства терригенных осадков»

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

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

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

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

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

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

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

Первая задача заключалась в создании алгоритмов и программ моделирования гравитационного осадка с заданной гранулометрической характеристикой.

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

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

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

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

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

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

Впервые получены эталоны случайного гравитационного осадка, как для моноразмерного (пористость 41.8%, координационное число 6, средний поперечник пор 69.5% от диаметра слагающих частиц), так и для полиразмерного.

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

Предложен принцип и алгоритм определения принадлежности частиц осадка к каркасу или матриксу.

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

Введено понятие расстояния для порового пространства на основе трех вариантов метрик порового пространства.

Новой и перспективной является схема описания структуры горной породы (осадка), в основу которой положена заимствованная из программирования информационная структура - запись.

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

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

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

Апробация и публикации. Результаты исследований опубликованы в 11 научных работах. Основные положения диссертационной работы докладывались на конференциях: всероссийских - XIII Геологический съезд Республики Коми (Сыктывкар, 1999); Юбилейная конференция к 70-летию ВНИГРИ (Санкт-Петербург, 1999); "Моделирование геологических систем и процессов" (Пермь, 1996); региональных - IV, V, VI, VII, VIII научных конференциях Института геологии Коми НЦ УрО РАН (Сыктывкар, 1995-1999); XIII Коми республиканская молодежная научная конференция (Сыктывкар, 1996); X конференция, посвященная памяти К.О.Кратца (Апатиты, 1999).

Объем работы. Диссертация состоит из введения, пяти глав и заключения, содержит 126 страниц текста, 53 рисунка, 6 таблиц, список использованной литературы из 97 наименований и 4 приложения.

Заключение диссертации по теме «Общая и региональная геология», Амосова, Ольга Евгеньевна

ЗАКЛЮЧЕНИЕ

Данная работа посвящена созданию формальной математической теории структуры осадочной терригенной породы. Работа состоит из экспериментальной и исследовательской частей. Экспериментальная часть заключается в компьютерном моделировании осадка. Исследовательская - в его изучении.

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

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

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

Для изучения коллекторских и других структурных свойств было проведено компьютерное моделирование моноразмерных и двухразмерных моделей осадка, значение отношения размеров крупных и мелких частиц в которых изменялось от 2 до 10. При этом исследовались осадки как насыщенные, так и недосыщенные и пересыщенные в 1.5-2 раза мелкой фракцией. Для моделирования двухразмерных осадков, с заданным значением насыщенности мелкой фракцией, вычислялись количества крупных и мелких частиц, необходимые для заполнения контейнера определенного объема, используя среднее значение пористости моноразмерного осадка, следующим образом. Количество крупных частиц определяется как частное от деления объема части контейнера (заполненного моноразмерным осадком), приходящейся на частицы моноразмерного осадка, на объем крупной частицы. Количество мелких частиц - частное от деления объема части пустотного пространства контейнера (заполненного моноразмерным осадком), приходящейся на частицы моноразмерного осадка, на объем мелкой частицы.

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

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

Для изучения влияния силы гравитации на строение осадка, в частности, на его анизотропию, исследовался средний поперечник поровых каналов в трех направлениях: вертикальном (направление силы гравитации) и двух перпендикулярных ему и друг другу. В результате исследований установлено, что с увеличением показателя насыщенности мелкой фракцией и значения отношения диаметров крупных и мелких частиц двухразмерного осадка заметно уменьшается поперечник поровых каналов. В горизонтальном сечении он на 15-50% меньше, чем в двух других перпендикулярных ему. Это особенно характерно для недосыщенных осадков с большим значением отношения диаметров крупных и мелких частиц. Этот эффект хотя и слабо, проявился даже в моноразмерном осадке (69.38±0.28% в горизонтальных сечениях против 69.68±0.75% в вертикальных).

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

118 1

Список литературы диссертационного исследования кандидат геолого-минералогических наук Амосова, Ольга Евгеньевна, 2000 год

1. Абрикосов И.Х., Гутман И.С. Общая, нефтяная и нефтепромысловая геология. -Москва: Недра, 1974. -360 с.

2. Амосова O.E. Координационное число и упаковка изометричных частиц гравитационного осадка // Геология и минеральные ресурсы европейского северо-востока России: Материалы XIII Геологического съезда Республики Коми. -Сыктывкар, 1999. Т. IV. - С. 227-228.

3. Амосова O.E. Структура порового пространства осадка // Геология и полезные ископаемые северо-запада и центра России: Материалы X конференции, посвященной памяти К.О. Кратца. -Апатиты, 1999. -С. 184-186.

4. Амосова O.E. Модель гравитационной седиментации терригенного осадка // Структура, вещество, история литосферы Тимано-Североуральского сегмента: Материалы 6-ой научной конференции Института геологии Коми НЦ УрО РАН. -Сыктывкар: Геопринт, 1997. -С. 9-11.

5. Амосова O.E. Основные моменты моделирования пористости осадка с заданной гранулометрической кривой // Материалы XIII Коми республиканской молодежной научной конференции. -Сыктывкар,1997. С. 57.

6. Амосова O.E., Ткачев Ю.А. Компьютерное моделирование гравитационной седиментации изометричных частиц с заданной гранулометрической кривой // Моделирование геологических систем и процессов: Информ. материалы. Пермь: ПТУ, 1996.-С. 31-32.

7. П.Белкин В.И., Рязанов И.В. О понятиях сортированности и однородности гранулометрического состава осадочных пород // Литология и полезные ископаемые. -Москва, 1974. -N2. -С. 133-139.

8. Белощенко Д.К.// Расплавы.-1987 -Т.1, N3. -С.45.

9. Бурлаков И.А. Новые фильтрационно-емкостные параметры гранулярных горных пород // Геология нефти и газа. -Москва, 1974. NX. -С. 60-64.

10. Войтеховский Ю.Л. Минеральный агрегат: топологии, метрики, ковариации, оценивание // Препр.: Геол. Ин-т Кольск. Научн. Центра РАН. — Апатиты, 1998. 40 с.

11. Выгодский М.Я. Справочник по высшей математике. -Москва, 1963. -872 с.

12. Гольдберг В.М., Скворцов Н.П. Проницаемость и фильтрация в глинах. -Москва: Недра, 1986. -160 с.

13. Гриффите Д. Научные методы исследования осадочных пород. -М.: Мир, 1971.-410 с.

14. Гришин Ф.А. Промышленная оценка месторождений нефти и газа. -Москва: Недра, 1975. -304 с.

15. Гудок Н.С. Изучение физических свойств пористых сред. -Москва: Недра, 1970. -208 с.

16. Дзисько В.А., Карнаухов А.П., Тарасова Д.В. Физико-химические основы синтеза окисных катализаторов. -Новосибирск: Наука, 1978. -384с.

17. Жданов М.А., Гординский Е.В., Ованесов М.Г. Основы промысловой геологии газа и нефти. -Москва: Недра, 1975. -295 с.

18. Жданов М.А. Нефтегазово-промысловая геология и подсчет запасов нефти и газа. -Москва: Недра, 1981. -453 с.

19. Заграфская Р.Ф., Карнаухов А.П.,Фенелонов В.Б. Глобулярная модель пористых тел корпускулярного строения. П1. Исследование случайных и частично упорядоченных упаковок шаров "Кинетика и катализ", 1975. Т.16, N6,с. 1583-1590.

20. Иванов В.А., Храмова В.Г., Дияров Д.О. Структура порового пространства коллекторов нефти и газа. -Москва: Недра, 1974. -96 с.

21. Исследование структуры пор в компьютерных моделях плотных и рыхлых упаковок сферических частиц / Волошин В.П., Медведев H.H., Фенелонов В.Б., Пармон В.Н. // Журнал структурной химии. -1999. -Т.40, N 4. -С. 681-692.

22. Кноринг Л.Д. Условия формирования пористости песчано-алевритовых пород и их выявление методом факторного анализа // Литология и полезные ископаемые. Москва, 1974. -N2. -С. 66-77.

23. Коваленко Э.К., Халимов Э.М. К вопросу о взаимосвязи проницаемости с пустотностью и удельной поверхностью при учете структуры породы. -Труды УфНИИ, 1967,вып. ХУЛ. -С. 60-65.

24. Компьютерная порометрия монодисперсных силикагелей / Волошин В.П., Медведев H.H., Фенелонов В.Б., Пармон В.Н. // Доклады АН. -1999. -Т.364, N3. -С. 337-341 .

25. Конюхов И.А. О зависимости петрофизических свойств терригенных пород от их литологии // Вестник московского университета. -Москва, 1961. -N3. -С. 57-70.

26. Котяхов Ф.И. Количественная оценка структуры поровых каналов пород нефтяных и газовых залежей // Новости нефтяной и газовой техники. Нефтепромысловое дело, 1962. N1. -С. 35-38.

27. Котяхов Ф.И. Оценка однородности и ее связи с коллекторскими свойствами горных пород // Геология нефти и газа, 1971.7V4, С. 26-29.

28. Котяхов Ф.И. Результаты анализа методов определения среднего радиуса пор пород нефтяных и газовых залежей // Новости нефтяной и газовой техники. Нефтепромысловое дело, 1962. N2. -С. 37-42.

29. Котяхов Ф.И. Физика нефтяных и газовых коллекторов. -М.: Недра, 1977.-287 с

30. Куренков А.И. О некоторых методах определения удельной поверхности и проницаемости гранулярных коллекторов // Литологические исследования пород-коллекторов в связи с разведкой и разработкой нефтяных месторождений.-Москва:Наука, 1970. -С.52-57.

31. Леворсен А. Геология нефти и газа. -М.: Мир, 1970, -638 с.

32. Лидер М.Р. Седиментология. -М.: Мир, 1986. -439 с.

33. Математическая статистика. Под ред. A.M.Длина. -Москва: Высшая школа, 1975. -398 с.

34. Медведев H.H. Метод Вороного-Делоне в исследованиях структуры неупорядоченных систем. Докторская диссертация, 1996, Новосибирск, ВЦ СО РАН.

35. Медведев H.H. Приложение метода Вороного-Делоне к описанию структуры межшарового пространства в полидисперсных системах // Доклады АН. -1994. -Т.337,7V6. -С.767-771."

36. Медведев H.H., Волошин В.П., Наберухин Ю.И. Структура простых жидкостей как перколяционная проблема на сетке Вороного // Журнал структурной химии. -1989. -Т.30,7V2. -С. 98.

37. Медведев H.H., Наберухин Ю.И. Исследование структуры простых жидкостей и аморфных тел методами статистической геометрии // Журнал структурной химии. -1987. -Т.28, N 3. -С. 117.

38. Методика исследования фильтрационно-емкостных параметров горных пород на образцах больших размеров/ Гороян Г.Ф., Рабиц Э.Г., Петерсилье В.И., Белов Ю.Я. // Методика разведки и подсчета запасов месторождений нефти и газа-Москва, 1977. -С. 146-154.

39. Могилевская Т.М., Куренков А.И. Определение извилистости поровых каналов // Цитологические исследования пород-коллекторов в связи с разведкой и разработкой нефтяных месторождений. -Москва, 1975. -С. 5-9.

40. Наберухин Ю.И., Волошин В.П., Медведев H.H. Исследование межатомного пространства в моделях одноатомных систем с помощью методов Вороного-Делоне // Журнал физической химии. -1992. -Т. 66, N1. -С. 155-162.

41. Наберухин Ю.И., Волошин В.П., Медведев H.H. // Расплавы. -1987-Т.1, N2. -С.71.

42. Осовецкий Б.М. Дробная гранулометрия аллювия. Пермь: Изд-во Перм. ун-та, 1993. - 343 с.

43. Прошляков Б.К. Зависимость коллекторских свойств от глубины залегания и литологического состава пород // Геология нефти и газа, 1960. N12, -С. 24-29.

44. Прошляков Б.К., Кузнецов В.Г. Литология. -М.: Недра, 1991. -444 с.

45. Петтиджон Ф.Дж. Осадочные породы. -М.: Недра, 1981. -751 с.

46. Рожков Г.Ф. Геологическая интерпретация гранулометрических параметров по данным дробного ситового анализа // Гранулометрический анализ в геологии. -Москва, 1976.-С. 5-25.

47. Рухин Л.Б. Основы литологии. (Учение об осадочных породах). 3-е изд., перераб. и доп. -Л.: Недра, 1969. -703 с.

48. Седиментология / Градзиньский Р., Костецкая А., Радомский А., Унруг Р. -М.: Недра, 1980-Пер. Изд.ПНР, 1976. -640 с.

49. Словарь по геологии нефти и газа. Под ред. К.А.Черникова и др. -Ленинград: Недра, 1988. -679 с.

50. Справочник по теории вероятностей и математической статистике. Под ред. Королзок B.C., Портенко Н.И. Скороход A.B., Турбин А.Ф. -Москва: Наука, 1985. -640 с.

51. Ткачев Ю.А., Амосова O.E. Терригенный осадок.с изнанки // Вестник Института геологии. -1999. -N2- С. 3-5.

52. Ханин A.A., Бурова Е.Г. Зависимость проницаемости горных пород от структуры порового пространства // Разведка и охрана недр. -Москва: Недра, 1964. N12. - С. 22-27.

53. Ханин A.A. Основы учения о породах-коллекторах нефти и газа. -Москва: Недра, 1965.-360 с.

54. Ханин A.A. Петрофизика нефтяных и газовых пластов.-М.: Недра, 1976. -295 с.

55. Энгельгарт В. Поровое пространство осадочных пород. -Москва: Недра, 1964. -232 с.

56. Allen М.Р:, Tildesley D.J. Computer simulation of liquids. Oxford: Clarendon Press, 1987. P. 386.

57. Bieshaar R., Geiger A., Medvedev N.N. // Mol. Simulation 1995. - 15. - P. 189-196.

58. Bryant S., Blunt M. // Phys.Rev. 1992. - 46, N 4. - P. 2004-2011.

59. Coogan A.H., Manus R.W. Compaction and diagenesis of carbonate sands // Compaction of coarse-grained sediment. Vol. 1. Amsterdam-Oxford-New York. 1975. P. 79-166.

60. Dickinson, W.R., 1970, Interpreting detrital modes of graywacke and arkose: Jour. Sed. Petrology, v. 40, pp. 695-707.

61. Dinger D.R., Funk J.E.: Particle Packing I Fundamentals of Particle Packing, Monodisperse Spheres. Interceram41 (1992) 1. 10-14.

62. Dinger D.R., Funk J.E.: Particle Packing II Review of Packing of Polydisperse particle Systems. Interceram 41 (1992) 2. 95-97 and [3] 176-179.

63. Finney J.L., Wallace J. // J. Non-Crystalline Solids. 1981. V. 43. P. 165.

64. Ghidaglia C., Guazzelli, Oger L. // J. Phys. D: Appl. Phys. 1991. V. 24. P. 2111-2114.

65. Graton L.C:, Fraser H.J. Sistematic Packing of Spheres with Particular Relation to Porosity and Permeability, J. Geol., 43,785-909,1935.

66. Greenfield M. D., Theodorou D.N. // Macromolecules. 1993. - 26 - P. 961-996.

67. Greqq S.J., Sing K.S.W. Adsorption, surface area and porosity. L.: Acad. Press, 1982.

68. Guyon E., Oger L., Piona TJ. // J. Phys. D. 1987. V. 20. P. 1637.

69. Jullien R., Meakin P., Pavlovitch A. Disorder and Granular Media. Elsevier Science Publ., 1993.

70. Liu H., Lin Zhang, Seaton N.A. // Chem. Eng. Sci. 1992. V. 47. P. 4393.

71. Luchnikov B.A., Medvedev N.N., Naberukhin Yu.L., Novikov V.N. // Phys. Rev. B. -1995.-51.-P. 15569-15572.

72. Mason G. // J.Coll.Interface Sci.-l 972.-42, N2. -P.208-227.

73. Mason G., Mellor D.W. // J. Colloid and Interface Sci.-1995.-176(l).-P.214-225.

74. Meakin P., Jullien R. // J. Phys. France. 1990. V. 51. P. 2673-2680.

75. Medvedev N.N. // Ibid. 1986. V. 67. P. 223.

76. Medvedev N.N., Voloshin V.P., Naberukhin Yu.I. Physical Chemistry of Colloids and Interfaces in Oil Production. Paris, 1992.

77. Medvedev N.N., Naberukhin Yu.i. // J. Phys. A: Math. Gen. 1988. V. 21. P. L247-L252.

78. Oger L., Gervois A., Troadec J.P., Rivier N. // Philosophical Magazine B. 1996, - 74 (2). -P. 177-197.

79. Rydgway K., Turbuk K.J. // British Chem. Eng.-l967.-12, N3.-P.384-388.

80. Sastry S., Corti D.S., Debenedetti P.G., Stillinger F.H. // Phys. Rev. E. 1997. - 56(5). - P. 5524-5532.

81. TanakaH. // Chem. Phys. Lett. 1997. - 282. - P. 133-138.

82. TanemuraM., OgawaT., OgitaN. III. Comput.Phys. 1983. V. 51. P. 191.

83. Thompson K.E., Fogler H.S. // AIChE Journal. 1997. V. 43(6). P. 1377-1389.

84. Trask P.D., Suggested Classification of Detrital Particles, Natl. Res. Council Rept. Comm. Sedimentation 1941-1942, appendix G, p.3, 1942.

85. VoloshinV.P., Naberukhin Y.I., Medvedev N.N. // Molecular Simulation. -1989. V.4. -P. 209.

86. Voronoi G.F. // J. Reine Angew. Math. 1909. V.136. P.67-181.

87. Zhdanov V.P., Fenelonov V.B., Efremov D.K. // J. Colloid. Interface Sci. 1987. - 120. -P. 218-231.