ФИЗИЧЕСКАЯ ЭКВИВАЛЕНТНОСТЬ ГОРНЫХ ПОРОД И ПЛАСТОВЫХ ФЛЮИДОВ ПО НЕЙТРОННЫМ СВОЙСТВАМ
Результаты теоретических расчетов нейтронных характеристик горных пород можно непосредственно использовать для определения в реальной системе прибор — скважина — пласт изменений показаний нейтронных методов, обусловленных изменениями минерального состава скелета пласта, его плотности, глинистости, концентрации элементов-поглотителей, водо-, нефте-и газонасыщенности. Возможность такого эффективного использования нейтронных характеристик обусловлена принципом физической эквивалентности горных пород и пластовых флюидов.
Принцип физической эквивалентности можно сформулировать следующим образом: при фиксированных технических условиях измерений в скважине породы, имеющие различный состав и петрофизические параметры ( , , , , и т. д.), отмечаются одинаковыми показателями, если совпадают их нейтронные характеристики.
Поля излучения в сравниваемых системах могут оставаться физически эквивалентными и при изменении нейтронных характеристик этих систем. В таких случаях в соответствии с принципом подобия трансформируются по известному правилу геометрические размеры и параметры всех зон [15].
Способ моделирования, основанный на физической эквивалентности пород, впервые был предложен и применен Е. И. Крапивским и А. А. Бремом [49]. Эти исследователи показали, что при экспериментальных исследованиях можно добиться совпадения нейтронных свойств породы (Ls, Pa, D, t) с натурной моделью, добавляя в насыпную среду-наполнитель (кварцевый песок, мраморную муку) вещества, хорошо рассеивающие (вода, парафин, полиэтилен, бериллий, углерод) и поглощающие (бор, литий, ртуть и др.) нейтроны, а также изменяя плотность модельной среды добавлением специальных утяжелителей (свинцовой или чугунной дроби, барита, галенита и т. д.). Соотношение компонентов модельной смеси подбирается таким образом, чтобы ее нейтронные характеристики: полный возраст нейтронов и вероятность избежать поглощения при замедлении, сечения рассеяния и поглощения тепловых нейтронов совпадали с характеристиками изучаемой среды (породы). При одинаковой весовой влажности и плотности сред параметры Ls, , D оказываются близкими. Добавление поглотителей нейтронов позволяет добиться совпадения параметров Ра,t, характеризующих поглощающие свойства.
Изменение водородосодержания (пористости) горных пород существенно влияет на величину нейтронных характеристик, как правило, значительно сильнее, чем изменение их химического состава.
Длина замедления быстрых нейтронов зависит от водородосодержания через замедляющую способность (и средний косинус угла рассеяния).
Однако погрешность определения т по величине , обусловленная неконтролируемым изменением минералогического состава пласта, может быть значительной. Так, =12 см соответствует песчанику с т = 21 % или доломиту с т =16 %, или ангидриту и известняку с . Соответствующая максимальная относительная погрешность определения объемного водонасыщения по длине замедления достигает 30 % отн.
Коэффициент диффузии с изменением минералогического состава изменяется почти так же, как и длина замедления. Время жизни тепловых нейтронов в первую очередь определяется концентрацией сильно поглощающих элементов. При малых концентрациях аномальных поглотителей на величину t оказывает заметное влияние также водородосодержание породы.
Из анализа нейтронных характеристик и закономерностей поля излучения в системе скважина — пласт следует, что при количественном определении водородосодержания (водо-, нефте- и газонасыщенности) стационарными методами нейтронометрии скважин критерием физической эквивалентности служит совпадение длин замедления нейтронов в пласте (для ННМнт) или длин миграции тепловых нейтронов (для ННМт). Обоснованность этого вывода можно проиллюстрировать результатами прямых расчетов методом Монте-Карло, выполненных Ф. X. Еникеевой и Б. К. Журавлевым для различных литотипов при реальной скважинной геометрии (рис. 74), а также экспериментальными данными зарубежных исследователей (рис.75) для одозондового (а) и двухзондового (б) приборов.
В системе скважина-пласт длина переноса захватного гамма-излучения для нейтронного гамма-метода физического смысла не имеет, ине может служить даже приближенным критерием физической эквивалентности пластов различного состава. При одинаковых { , , }разные зоны в системе прибор — скважина—пласт могут сильно различаться по величине излучающей способности и соотношению вкладов в результирующую зависимость . Если основной вклад дает излучение скважины (например, при регистрации жесткого гамма-излучения с кадмиевым фильтром, усиливающим поглощение тепловых нейтронов в корпусе прибора), то длина миграции тепловых нейтронов может служить приближенным критерием физической эквивалентности и для НГМ.
Рис.10.7. Обратные величины показаний ННМнт и ННМт (в единицах показаний в пресной воде) в породах различного состава как функции водонасыщенной пористости, выраженной через длину замедления (для ННМнт, зонд Z= 47,5 см; кривая 1) и длину миграции М (для ННМт, зонд Z = 50 см; кривая 2) (по данным Ф.X.Еникеевой и Б.К.Журавлева, расчет методом Монте-Карло).
Скважина необсаженная, dc = 200 мм, прибор ДРСТ прижат к стенке скважины.
Шифр кривых: 1 — известняк; 2 — кварцевый песчаник; 3 — доломит; 4 — полимиктовый песчаник.
Еще один тип инверсионного эффекта для нейтронных методов, впервые отмеченный Ю.А.Гулиным (рис.76), можно назвать «эффектом литологической инверсии», поскольку он связан с влиянием минерального состава.
Зависимость показаний нейтронных (ННМнт и ННМт) и нейтронного гамма-метода можно описать экспонентами, декременты пространственного затухания в которых определяются обратными длинами замедления надтепловых или миграции тепловых нейтронов (λ):
Рис.10.8 а,б. Показания однозондового прибора ННМнт (по Г.Эдмундсону и Л.Раймеру) а): 1 – кварцевый песчаник; 2 – доломит; 3 – известняк; 4 – вода; и интерпретационного параметра двухзондового прибора ННМнт (по Г.Скотту, С.Флауму и Г.Шерману) б): 1 –песчаник, 2 – известняк, 3 – доломит, 4 – насыщение пресной водой, 5 – насыщение соленой водой; 6 – газ при различных соотношениях газ/вода. как функция длины замедления нейтронов для различных литотипов. Длины замедления рассчитаны теоретически. Необсаженная заполненная водой скважина диаметром 20 см.
N(z) ~ exp(-z/λ). (10.271)
Пусть в изучаемом геологическом разрезе присутствуют пласты различных литотипов, нейтронные характеристики которых λ изменяются в неперекрывающихся интервалах. Перенумеруем литотипы пород в последовательности изменения величин λi , выберем пласт «опорной литологии» λоп и сформируем интерпретационный параметр в виде отношения
Ni(z)/ Nоп(z) ~ exp[z(1/λоп-1/λi)] ~ exp(zDi) , (10.272)
где через Di обозначена разность
Di = 1/λоп-1/λi. (10.273)
При изменении длины зонда (это управляемый параметр) величина интерпретационного параметра ( ) на диаграмме рис.9.76 будет возрастающей при Di > 0 или убывающей при Di < 0, причем скорость роста или убывания определяется величиной модуля |Di|.
Размер зонда z, см
Рис.10.9.«Литологическая инверсия» показаний ННМ и НГМ (по Ю.А.Гулину). Зависимость относительных показаний ННМнт и НГМ от длины зонда в пластах разного минерального состава. 1 и 2 – ангидрит; 3 и 4 – каменная соль; 5 и 6 – пористые известняки; 1,3,5 – НГМ; 2,4,6 – ННМнт. За единицу приняты показания в плотных известняках.
Необходимость использования физической эквивалентности пластовых флюидов (воды и газа) впервые возникла в связи с задачей определения газонасыщенности коллекторов по данным нейтронометрии скважин. В случае двухфазного флюида (вода, газ) эффективное водонасыщение неглинистого пласта определяется соотношением [6]
(10.274)
где р — давление газа; р0 — давление, при котором плотности ядер водорода в воде и в сжатом газе совпадают (при t = 20 °C для воды без растворенного газа р0 = 62,3 МПа). Это соотношение устанавливает эквивалентность флюидов по плотности ядер водорода. Однако совпадение не обеспечивает совпадения параметров замедления , и длины диффузии тепловых нейтронов, т. е. не является условием физической эквивалентности пластовых флюидов (исключение составляет время жизни тепловых нейтронов, поскольку микросечение поглощения у водорода на два порядка больше, чем у углерода, и на три порядка больше, чем у кислорода). Физическая эквивалентность водородного газа и воды в действительности выполняется при давлении , слабо зависящем от начальной энергии нейтронов.
Поскольку нейтронные и термодинамические характеристики углеводородного и водородного газов различаются, физическая эквивалентность воды (нефти) и природного газа будет иметь место при давлении р0:
(10.275)
где q — атомная единица массы; R — универсальная газовая постоянная; v — числа атомов в молекуле газа; Т — температура; Z — коэффициент сжимаемости; п — ядерная плотность газа; и — соответственно возраст (или площадь миграции) нейтронов в газе и воде (нефти) в пластовых условиях.
При одинаковом водородосодержании (точнее, при физической эквивалентности флюидов) газоносные и водоносные коллекторы не удовлетворяют условию физической эквивалентности, поскольку они отличаются объемом твердой фазы (плотностью). Этот плотностной эффект называют эффектом замещения или экскавационным эффектом [107] .
В задачах рудной геологии Е. И. Крапивский предложил при количественном определении концентраций нейтронных поглотителей пользоваться понятием эквивалентной концентрации. Эквивалентными называются такие составы и концентрации различных элементов-поглотителей, которым соответствуют одинаковые макросечения поглощения.
Принцип физической эквивалентности и принцип подобия взаимно дополняют друг друга как способы экспериментального и теоретического исследования в задачах ядерной геофизики. При ограниченности палеточного фонда и натурных моделей пластов принцип физической эквивалентности позволяет распространить имеющиеся эталонировочные зависимости, получаемые в результате математического и физического моделирования, на большое число разнообразных геологических ситуаций, существенно повышая эффективность геологической интерпретации результатов измерений нейтронными методами.
[1] Восстановление истинных закономерностей изменения полей по интегральным характеристикам. Изв. АН СССР, сер. Геофизич. 1962, №7.
[2] Не вдаваясь в петрохимические детали (существенные при интерпретации данных нейтронометрии), под этим термином здесь и далее понимается водонасыщенная пористость гранулярного коллектора.
[3] Д.А.Кожевников. Атомная энергия. 1964, № ,
[4]расcчитано А. И. Пшеничнюком и Б. Е. Лухминским.
[5] В связи с тем, что неправомерно распространение закона рассеяния надтепловых нейтронов в область тепловых энергий.
[6] Для максвелловского спектра .
[7] В Башкирском государственном университете (задача поставлена Ю. А. Гулиным).
[8] Предполагается устойчивость модальных или средних значений распределений пористости по данным межскважинной корреляции.
[9] на которые нормируются результаты расчетов методом Монте-Карло.
[10] Геофизика, 2005, № 4, с. 31 – 37.
[11] Подчеркнем, что приведенные материалы не претендуют на роль инструкции по интерпретации.
[12] Мартьянов И. А., Старцев А. А., Шиканов А. Е., Федына Е. А., Рудов И. В., «Двухкомпонентная модель показаний импульсного нейтронного каротажа», Геоинформатика, М., 1998, № 2, 31 – 34.
Дата добавления: 2020-08-31; просмотров: 484;