Text
                    С. С. ИТЕНБЕРГ, Т. Д. ДАХКИЛЬГОВ
ГЕОФИЗИЧЕСКИЕ
ИССЛЕДОВАНИЯ
В СКВАЖИНАХ
ДОПУЩЕНО
Министерством высшего и среднего
специального образования СССР
в качестве учебника для студентов
вузов, обучающихся по специальности
«Геология и разведка
нефтяных и газовых месторождений»
МОСКВА «НЕДР А» 1982


УДК 550.832 @75) Итенберг С. С, Дахкильгов Т. Д. Геофизические исследования в скважинах. М., Недра, 1982. 351 с. Описаны физические свойства горных пород и связи между ними. Изложены теория, техника и методика электрических, радиоактивных, аку- акустических, термических и других геофизических методов, изучение технического состояния нефтя- нефтяных и газовых скважин, контроль процессов раз- разработки месторождений, автоматические средства преобразования геофизической информации, тех- техника безопасности и охрана окружающей среды. Для студентов нефтяных и геологоразведочных институтов по курсу «Геология и разведка нефтя- нефтяных и газовых месторождений». Табл. 16, ил. 193, список лит.— 12 назв. Рецензенты: инж. Л. Л. My хер (МНП), кафедра геофизики АзИнНефтехим им. М. Азизбекова 1904050000-262 043 @1) — 82 © Издательство «Недра», 1982
ВВЕДЕНИЕ Поиски и разведка полезных ископаемых сопровождаются бу- бурением скважин. В нефтяной и газовой промышленности последние служат не только лля поисков и разведки месторождений нефти и газа, но и для их разработки. В целях изучения геологического разреза скважин, их технического состояния и контроля за режи- режимом разработки месторождений в них проводят геофизические ис- исследования (ГИС). Комплекс ГИС назван промысловой геофизикой. Изучение геологического разреза скважины заключается в опре- определении последовательности и глубины залегания пластов, их ли- толого-петрографических свойств, оценке наличия и количествен- количественного содержания в недрах полезных ископаемых. Изучение разреза возможно путем отбора керна. Однако он связан с применением специальных долот, что приводит к ограничению и замедлению бурения, а в конечном счете — к его удорожанию. Кроме того, керн не всегда удается извлечь из нужного интервала (неполный вынос керна), а при его отборе и выносе на поверхность свойства породы и насыщающей ее жидкости заметно изменяются, поэтому результаты анализа керна и шлама не дают полного представления о геологическом разрезе. Вместе с тем некоторые физико-химиче- физико-химические свойства пород (электропроводность, электрохимическая ак- активность, радиоактивность, температуропроводность, упругость и др.) поддаются изучению непосредственно в скважине в условиях их естественного залегания путем проведения в ней соответствую- соответствующих геофизических или геохимических исследований. Такие ис- исследования, заменяющие частично или полностью отбор керна, названы каротажем; результаты их изображаются в виде диаграммы изменения физических свойств пород вдоль скважин — каротажных диаграмм. В зависимости от изучаемых свойств горных пород известны следующие виды каротажа: электрический, радио- радиоактивный, термический, акустический и др. Результаты каротажа дают возможность определять последо- последовательность и глубину залегания пластов, вскрываемых скважи- скважиной, их литологические свойства и содержание в них полезных ис- ископаемых (нефть, газ, уголь, руды, нерудное сырье). Получен- Полученные данные являются исходными для изучения геологического строения месторождения и региона в целом, а также для подсчета запасов и проектирования рациональной системы разработки за- залежи. Геофизические данные являются в настоящее время основ- основными и служат для оценки коллекторских свойств пород и степени их насыщения нефтью, газом или водой. Отбор керна в таких сква-
жинах доводится до оптимального минимума, а в тех случаях, когда разрез месторождения хорошо изучен, бурение ведется без отбора керна. Однако полностью отказаться от него, особенно в разведочных скважинах, нерационально. Данные о пористости, проницаемости, глинистости, нефтегазонасыщенности и других свойствах, полученные при анализе керна, являются часто исход- исходными для корректировки результатов интерпретации материалов геофизических исследований. Контроль за разработкой нефтяных и газовых месторождений включает комплекс геофизических исследований в действующих и контрольных (одиночных) скважинах, размещенных в пределах эксплуатируемой залежи для изучения процесса вытеснения нефти в пласте и закономерностей перемещения водонефтяного, газонеф- газонефтяного и газоводяного контактов. Изучение технического состояния скажин проводится на всех этапах их действия: бурения, перед вводом в эксплуатацию, в про- процессе эксплуатации. Во время бурения инклинометром определяют искривления ствола скважины, каверномером — ее диаметр, ре- зистивиметром и электрическим термометром — места поступле- поступления жидкости из пласта в скважину и поглощения промывочной жидкости. Перед вводом скважины в эксплуатацию проводится изучение технического состояния колонны на герметичность и ка- качество цементирования. В эксплуатационных скважинах контроль их технического состояния заключается в выявлении мест наруше- нарушения герметичности цементного кольца, нарушений сцепления це- цемента с колонной и породой, вызывающих возникновение затруб- ной циркуляции жидкости. К ГИС принято также относить прострелочно-взрывные работы, опробование пластов приборами на кабеле, отбор керна боковыми грунтоносами, перфорацию колонн при вскрытии пластов, обса- обсаженных трубами, торпедирование. Связь этих работ с геофизиче- геофизическими вызвана тем, что для их выполнения применяется то же оборудование, что и при ГИС. В процессе бурения производят отбор образцов пород из стенки скважины боковым грунтоносом для дополнительного изучения вскрываемого разреза, отбор жидкостей и газов приборами на ка- кабеле (ОПК) для оперативного опробования пластов в необсаженных скважинах, торпедирование бурильного инструмента и металличе- металлических предметов, оставленных в скважине при ликвидации аварий. При вводе скважины в эксплуатацию для сообщений скважины с пластом осуществляют перфорацию колонны против продуктив- продуктивных горизонтов. В эксплуатационных и нагнетательных скважинах с открытым забоем с помощью пороховых генераторов давления и торпедиро- торпедирования производят разрыв пластов и тем самым повышают их от- отдачу или приемистость. Геофизические исследования являются в настоящее время неотъемлемой частью геологических, буровых и эксплуатационных работ, проводимых при разведке и разработке
нефтяных, газовых и других месторождений полезных ископаемых. Первые геофизические исследования в виде геотермических из- измерений в скважинах были выполнены в 1906 г. Д. В. Голубятни- ковым на месторождениях Биби-Эйбат, Сураханы (Баку). Широ- Широкое применение ГИС началось после внедрения электрического каротажа методом сопротивления, предложенного в 1928 г. фран- французским ученым Шлюмберже. В СССР электрический каротаж (замер кажущихся сопротивлений КС) был впервые проведен в 1929 г. в Грозном по инициативе И. М. Губкина и Д. В. Голубят- никова; с 1930 г. эти исследования начали проводиться в Баку. В 1931 г. одновременно с замером КС стали измерять и потенциал естественного поля в скважине (ПС). Развитию электрического каротажа и превращению его в современный метод изучения раз- разрезов скважин способствовала плодотворная работа коллектива советских геологов (В. В. Денисевича, В. А. Долицкого, В. А. Сель- Сельского, А. Н. Снарского, Г. М. Сухарева, М. А. Жданова и др.) под руководством И. М. Губкина и Д. В. Голубятникова. Разработка электрического каротажа-, его совершенствование и внедрение про- проводились весьма активно советскими учеными-геофизиками Л, М. Алыпиным, В. Н. Дахновым, С. Г. Комаровым, С. Я. Литви- Литвиновым и др. В результате электрический каротаж уже в 1933 г. начал применяться на всех нефтепромыслах СССР. Наряду с электрическим каротажем шло развитие и других ме- методов ГИС. В 1933 г. советским» геофизиками В. А. Шпаком, Г. В. Горш- Горшковым, А. Г. Граммаковым и Л. М. Курбатовым была начата раз- разработка гамма-каротажа. В 1937 г. ими были получены первые кривые. В 1941 г. известным физиком Б. М. Понтекорво был пред- предложен нейтронный каротаж. Радиоактивный каротаж в СССР по- получил широкое развитие с 1951 г. в связи с освоением аппаратуры и методики, разработанных Б. Б. Лапуком, Г. Н. Флеровым, Л. С. Полаком, Д. Ф. Беспаловым и другими в МИНХиГП им. И. М. Губкина. С 1959 г. получило развитие новое направление, разработанное под руководством Г. Н. Флерова, основанное на использовании импульсных источников нейтронов (ИННК). Опытные исследования по газовому каротажу были впервые проведены в 1933—1934 гг. М. В. Абрамовичем и М. И. Бальзамо- Бальзамовым. Промышленное внедрение этого метода началось с 1956 г. в связи с работами, выполненными Ю. М. Юровским, Л. М. Поме- Померанцем и др. Первые попытки измерения скорости распространения упругих волн в скважинах на малых базах были осуществлены в 1937 г. в СССР В. С. Воюцким. С 1955 по 1969 г. разработка акустического каротажа проводилась в ИФЗ АН СССР под руководством Е. В. Ка- руса, а в 1958 г. был создан первый макет трехэлементной уста- установки для дискретного ультразвукового каротажа. Специализиро- Специализированная установка ЛАК-1 была выполнена в конце 50-х годов, (Г. В. Дахнов, ВНИИГеофизика; А. А. Перельман, ВИРГ). В те- 5
чение 1970—1975 гг. была создана аппаратура акустического ка- каротажа для исследования как необсаженных, так и обсаженных скважин, а также для оценки качества их цементирования (Е. В. Ка- рус, О. Л. Кузнецов и др.). Наряду с развитием и внедрением промысловой геофизики в нефтяную промышленность, геофизические методы находят все более широкое применение при разведке угольных и рудных место- месторождений. Одновременно с развитием геофизических методов исследова- исследований разрезов скважин разрабатывались методы изучения их тех- технического состояния. В 1932 г. начаты работы по определению ази- азимута и угла искривления скважин при помощи инклинометров. В 1933—1934 гг. в нефтяных районах СССР начали широко при- применять термометрию для определений интервалов затрубного дви- движения жидкости, высоты подъема цемента за трубами и мест при- притока воды в скважину как термометром, так и резистивиметром. В последующие годы для решения этих задач начали использовать радиоактивные и акустические методы каротажа. В 1936 г. для измерения диаметра скважины был предложен каверномер. В это же время в цикл промыслово-геофизических работ были включены отбор керна боковым грунтоносом из стенок скважин, перфорация колонн при вскрытии пластов, торпедирование. В 1956—1959 гг. под руководством Н. Г. Григоряна были разработаны и внедрены кумулятивные (беспулевые) перфораторы. С 1970 г. под руководст- руководством П. А. Бродского разработан и внедрен метод опробования пла- пластов приборами на кабеле (ОПК). Начиная с 1950 г., запись каротажных диаграмм стала осущест- осуществляться автоматическими каротажными станциями, что улучшило качество замеров и повысило в несколько раз по сравнению с полу- полуавтоматической регистрацией скорость измерений. Этому способст- способствовало создание автоматических каротажных станций АКС (С. Г. Комаров и др.), ОКС (В. М. Запорожец и др.). В 60-е годы в СССР интенсивно разрабатывалась и внедрялась аппаратура новых методов исследования скважин: индукционного (С. М. Аксельрод, М. И. Плюснин) и бокового (В. Т. Чукия, A. Г. Мельников), для комплексирования измерений (А. Г. Бар- минский) и пластовых наклономеров (П. А. Зельцман, С. И, Кри- воносов и др.). К числу геофизических методов исследования скважин, полу- получивших развитие в СССР в 70-е годы, относятся высокочастотный электромагнитный каротаж, разработанный Д. С. Даевым и СБ. Де- Денисовым, ядерный магнитный каротаж — С. М. Аксельродом, B. М. Запорожцем и др.** * Для обработки и интерпретации результатов измерений с по- помощью ЭВМ необходим перевод геофизических исследований сква- скважин на цифровую запись. Для этого разработаны и внедряются цифровые каротажные преобразователи-регистраторы («Триас» — В. IX Логвинов, ПЛК-4 — А. Г. Мельников, Н. Н. Сохранов);
создана автоматическая газокаротажная станция с цифровой за- записью АГКС-4А (Л. М. Померанц и др.). Основоположником теории электрического каротажа по методу сопротивления является академик В. А. Фок, давший уже в 1932 г. теоретическое решение задачи для пласта неограниченной мощ- мощности» пересеченного скважиной. Позднее, в 1938 г., Л. М. Аль- пиным и С. Г. Комаровым были разработаны основные положения теории электрического каротажа по методу сопротивлений, пред- предложен метод бокового каротажного зондирования (БКЗ) и разра- разработана методика его интерпретации. Большая роль в разработке интерпретации данных электрического каротажа принадлежит В. Н. Дахнову. Им дана методика интерпретации кривых кажуще- кажущегося сопротивления в пластах малой мощности и определения коллекторских свойств пород по данным удельных сопротивле- сопротивлений. Важное значение в разработке электрического каротажа имели экспериментальные работы, выполненные в Баку и Грозном в 1933— 1936 гг. (С. Г. Комаров, И. М. Коган, Л. А. Горбенко, Г. С. Мо- Морозов и др.), обосновавшие влияние нефтеводонасыщенности на величину удельного сопротивления горных пород, доказавшие роль диффузионно-адсорбционных потенциалов в возникновении естественных потенциалов в скважине. В это же время Л. М. Аль- пиным и С. М. Шейнманном была создана теория электрического поля ПС. Исследования естественных потенциалов в скважине по- получили дальнейшее развитие в работах А. М. Нечая A952 г.) и Б. Ю. Вендельштейна A957—1959 гг.). Разработка физико-математических основ термических иссле- исследований в скважине была начата А. И. Заборовским A932 г.) и за- затем продолжена В. Н. Дахновым, Д. И. Дьяконовым и др. Первые исследования по теории нейтронного каротажа проводились А. И. Заборовским и позднее Ю. П. Булашевичем и С. А. Канто- Кантором, иа В области изучения петрофизических зависимостей горных по- пород, в том числе с моделированием естественных пластовых усло- условий, большой объем работ выполнен на кафедре промысловой гео- геофизики МИНХиГП им. И. М. Губкина В. Н. Дахновым, В. М. Доб- Добрыниным, В. Н. Кобрановой (I960—1980 гг.). Значительные ра- работы в области теории геофизических исследований скважин вы- выполнены зарубежными учеными Г. Доллем, Г. Гюйо, Р. Дебран- дом и др. Настоящий учебник составлен по программе курса «Геофизи- «Геофизические исследования в скважинах», утвержденной Учебно-методи- Учебно-методическим управлением по высшему образованию 8.IV.1976 г. для выс- высших учебных заведений по специальности 0103 — «Геология и раз- разведка нефтяных и газовых месторождений». В учебнике введение, главы III, кроме § 15, и IV, кроме § 22, главы V, VI, VIII, XI* XII, кроме § 56, и заключение составлены С. С.Итенбергом> глава I, §7,8главы II, § 15главы III, §22главы IV,
главы VII, X, §56 главы XII — Т. Д. Дахкильговым; § 5, 6, 9 главы II, глава IX—С. С. Итенбергом совместно с Т. Д. Дах- Дахкильговым, Авторы выражают искреннюю признательность сотрудникам кафедры промысловой и разведочной геофизики ГНИ им. акад. М. Д. Миллионщикова А. И. Соловьевой и Л. Г. Ишхановой за помощь в подготовке и оформлении материалов учебника.
ЧАСТЬ ПЕРВАЯ ПЕТРОФИЗИКА Горные породы в зависимости от условий их образования и рас- распространения обладают присущими им структурными и текстур- текстурными признаками. Они характеризуются определенным комплек- комплексом физических свойств — пористостью, проницаемостью, плот- плотностью, упругостью, удельным электрическим сопротивлением, радиоактивностью и др. Наука, занимающаяся изучением физиче- физических свойств горных пород и установлением численной взаимо- взаимосвязи различных параметров между собой, и называется п е т р о - ф и з и к о й. Глава I КОЛЛЕКТОРСКИЕ СВОЙСТВА И НЕФТЕНАСЫЩЕННОСТЬ ПОРОД § 1. ПОРИСТОСТЬ Под пористостью горной породы понимается совокупность пустот (пор) между частицами ее твердой фазы в абсолютно сухом состоянии. Пористость в основном определяет содержание в поро- породах жидкостей и газов и является одним из основных параметров, характеризующих их коллекторские свойства. Поры горных по- пород могут быть различны по происхождению, форме, размерам и внутренней взаимосвязи. Форма и размеры пор Формы и размеры отдельных пор, характер их распределения в породах и соотношение объема пор различных размеров, их взаи- взаимосвязь, извилистость и удельная внутренняя поверхность поро- вых каналов определяют структуру порового пространства породы. По своей форме поры пород крайне разнообразны и могут быть близкими к ромбоидальным и тетраэдрическим, щелевидным, ка- верновидным, трещиновидным, каналовидным, ячеистым, пузыр- пузырчатым и др. Примеры наиболее часто встречающихся типов пори- пористости пород представлены на рис. 1. По размерам поры условно подразделяются на сверхкапилляр- сверхкапиллярные, капиллярные и субкапиллярные. К сверхкапилляр- н ы м относятся поры радиусом 0,1 мм и выше. Движение жид- жидкостей и газов в них происходит под действием силы тяжести или напора по обычным для трубчатого канала законам гидродинамики.
Капиллярные поры имеют просвет от 0,0002 до 0,1 мм. Движение жидкостей и газов в них осуществляется при участии капиллярных сил и возможно лишь тогда, когда силы тяжести и напора значительно превосходят силы молекулярного взаимодей- взаимодействия между фильтрующимся флюидом и поверхностью канала фильтрации. Размер субкапиллярных пор составляет меньше 0,0002 мм. Вследствие малого расстояния между стенками, субкапиллярных каналов жидкость в них находится в сфере дейст- действия повеохностных молекулярных сил и при градиентах давлений. Рис. 1. Породы с различными типами пористости Уежзерновая пористость: а — с хорошо отсортированными зернами, б — с хорошо от- отсортированными зернами и цементирующим веществом в промежутках между ними, в — глинистый песчаник с рассеянным глинистым материалом; пористость* г — трещинно- каверновая, д — трещинная; / — зерна; 2 — глинистые частицы, 3 — цементирующий материал, 4 — блоковая часть породы наблюдаемых в природе, перемещаться не может. Породы, имею- имеющие в основном субкапиллярные поровые каналы, не являются коллекторами (глины, глинистые сланцы, сильноглинистые разно- разности терригенных пород и др.). Трещинная пористость классифицируется по раскрытое™ тре- трещин аналогично рассмотренной выше схеме. Трещины сверхкапил- сверхкапиллярны, если ширина их раскрытости больше 0,25 мм, капиллярны при ширине от 0,0001 до 0,25 мм, субкапиллярны при ширине мень- меньшей 0,0001 мм. Иногда выделяют макротрещины и микротрещины. Первые имеют раскрытость больше 0,1 мм, а вторые от 0,01 до 0,1 мм. Одна и та же порода может содержать поры различных разме- размеров. Их соотношение и распределение по объему в кластических породах зависит от отсортированное™, плотности укладки и формы породообразующих частиц. В практической работе поры горных пород по своей морфологии (форме и происхождению) принято де- делить на межзерновые (гранулярные), трещинные и каверновые. 10
Коэффициент пористости Поры горных пород могут быть взаимосвязанными и изолирован- изолированными. Первые соответствуют открытой, а вторые — закрытой ча- частям порового пространства породы. Общая (абсолютная) пористость породы равна суммарному объему открытых и закрытых пор. Количественно пористость оценивается коэффициентом, чис- численно равным отношению объема пор УПОр к объему образца по- породы Уобр 1 и выражается в долях единицы или процентах. Для оценки внутренней взаимосвязи пор используют: 1) коэффициент общей пористости kn = Упор, общ/Уобр» A.1) гяе Упор, общ — общий объем пор в образце породы; 2) коэффициент открытой пористости ?п.о= Упор. о/Уобр, A.2) где Fnop. о — объем открытых пор в образце породы; 3) коэффициент закрытой пористости &П.З — Упор. з/Уобр! A.3) где Упор, з — объем закрытых пор в образце породы. Статическая полезная емкость пород-коллекторов определяется объемом пор, которые могут быть заняты нефтью или газом. Вели- Величина этой пористости характеризуется коэффициентом эффектив- эффективной статической пористости Ъ л.— ^пор. о Упор, ев %. эф A.4) Уб ИЛИ kn. эф == A — ?в. ев) &п. о» A5) где УПор. ев — объем порового пространства, занятый связанной водой; Ав.св — коэффициент связанной водонасыщенности. При подсчете запасов нефти и газа используется средневзвешен- средневзвешенное значение коэффициента пористости, рассчитываемое по формуле k' (L6) где hi, St и knt — соответственно мощность, площадь и коэффи- коэффициент пористости отдельных участков пласта-коллектора. Для пород с жесткой связью частиц наряду с межзерновой по- пористостью (первичной) (kn. M) значительную роль играют поры вто- вторичного происхождения (Ап. вт) — трещины, каверны и другие пустоты выщелачивания. Суммарный объем пор в таких породах kn = ku. и + A — kn, M) ku. BT, (I 7) 11
откуда коэффициент вторичной пористости ъ ь 1 кп кп. м В зависимости от преобладающего типа пористости выделяют породы с межзерновой пористостью (поровые), трещиноватые, ка- кавернозные или порово-кавернозно-трещиноватые (смешанные). Влияние различных факторов на коэффициент пористости Величина пористости обломочных пород зависит от формы и раз- размера породообразующих частиц, степени их отсортированности, сце- ментированности и уплотненности. Породы с низкой пористостью (меньше 5 %) при отсутствии трещин и каверн обычно не являются промышленными коллекторами. Породы с пониженной пористостью характе- характеризуются kn = 5—10 %, со средней пористо- k 1015 % йk h -»L стью — kn = 10—15 %, с повышенной—kn = *•«$•* — 15—20%. Высокопористыми считаются по- *$\ .• роды с &п>20 %. в#*^ *в Увеличение глинистого и другого по ве- V»V#\ щественному составу цементирующего матери- # «"^* . ала снижает пористость • vVj#i пород (рис. 2). При от- • ^ в % сутствии в породе раство- ^^ •• римого компонента (кар- * • * e^^i *° бонатного и лимонитного) •*в »^^^^^._^ содержание глинистого ма- в •• •""• териала учитывается с по- помощью следующих основ - —'—fQ '—^ '—-6q А'гл'м,% ных параметров: 1) массовой гли- РисГ 2. ! Зависимость коэффициента пори- НИСТОСТИ кгл м ВЫ- стости^ от глинистости нижнемеловых тер- ражаЮщейся отношением ригенных пород Прикумского нефтегазо- г ^ ^ ^ш^пп^ш носного района Ставрополья массы ГЛИНИСТОЙ фракции (частиц размером меньше 0,01 мм) Мгл к полной массе жесткого скелета породы М: Ъ — АЛ IАЛ' /т Q\ «гл. м — ^гл/Л<1, (l.yj 2) об ъ е?н ой глинистости &гл, характеризующейся отношением объема глинистого материала (объема глинистых ча- частиц с прочно связанной с ними водой) к объему всей породы: &гл - VrJlI(VCK + УГл + Упор), A.9) где Угл, FCK, ^пор — объемы соответственно глинистой фракции, скелета и пор; 3) относительной глинистости т|гл, показы- показывающей степень заполнения глинистым материалом пор неактив- 12
ного скелета: ть = *гл/(*гл + ад- A.10) Если считать, что минералогические плотности вещества ске- скелетной и глинистой фракций равны, то справедливы выражения: kr3l — kTn. м A — ku) ИЛИ kn = &п.ск — &гл*> n, (I 11> 1 — /ггл. м где ?п.ск — пористость скелета породы; kn — открытая пористость. Если порода не содержит глинистого материала (см. рис. 1, а), ТО kn = &п. ск • В естественных условиях горные породы под действием всесто- всестороннего давления находятся в состоянии напряжения. В осадочных породах при определении среднего нормального напряжения до- допускают, что оно равно напряжению, созданному весом вышележа- вышележащих пород (горным или геостатическим давлением /?). Разность между горным давлением р и пластовым (поровым) риор называется эффективным давлением (напряжением) /?эф, дейст- действующим на скелет породы: РэЬ = Р—Рпор. A-12) Величина геостатического давления с учетом изменяющейся плотности пород определяется как h р= ЁяМп* A.13) или P = gHdn. ср, где Snt — плотность литологически однородного пласта мощностью Ы\ 6П> Ср — средневзвешенное значение плотности вышележащих пород до глубины Я определения геостатического давления; g — ускорение свободного падения. Величина АюР = ?#бв.ср, A.14) где бв. Ср — среднее значение плотности воды по разрезу. Соответст- Соответственно формула A.12) может принять вид Рэф = ?# (бп. ср — бв. ср). A.15) Под действием эффективного давления породы претерпевают как упругие, так и неупругие деформации, которые сопровождаются уменьшением пористости. Согласно экспериментальным данным под действием давления вышележащих пород глины, уплотняясь, отдают связанную воду и их пористость с глубиной уменьшается по экспоненциальному закону: *n^noexP(-P/4)- <1Л6> Здесь kn и kUo — пористость глин соответственно на заданной глубине и поверхности; $ — постоянная величина, характеризую- 13
щая степень уплотнения глин с глубиной, отвечает коэффициенту необратимого уплотнения породы; /?эф — эффективное давление (напряжение). Определение коэффициента пористости В лабораторных условиях обычно определяют коэффициенты открытой &п. о и общей kn пористости пород. Для определения ко- коэффициента открытой пористости используются различные (весо- (весовой, объемно-весовой, объемный или газометрический) способы [10]. Весовой способ (способ И. 'А. Преображенского) является наи- наиболее распространенным. Величина kut o этим способом опреде- определяется в следующей последовательности. 1. Экстрагированный от нефти образец породы помещают в бюксу и сушат в сушильном шкафу при температуре 105 °С до постоянной массы. 2. Определяют массу сухого образца взвешиванием в воздухе (Мг). Взвешивание производят на технических весах 1-го класса чувствительностью 10 мг. 3. Образец насыщают рабочей жидкостью под вакуумом, созда- создаваемым с помощью специальной вакуумной установки. Методика насыщения должна обеспечивать полное заполнение всех открытых (взаимосвязанных) пор образца рабочей жидкостью. 4. Определяют массу образца гидростатическим взвешиванием его в рабочей жидкости (М2). Для этого образец помещают в спе- специальную корзинку, подвешенную к серьге коромысла весов и по- погружающуюся в стакан с рабочей жидкостью, устанавливаемый на мостике. Массу корзинки Мк определяют гидростатическим взве- взвешиванием отдельно и вычитают из результатов взвешивания. 5. Определяют массу насыщенного рабочей жидкостью образца взвешиванием в воздухе (М3). Перед этим его обкатывают несколько раз на предметном стекле или слегка обтирают тряпочкой, чтобы снять с поверхности излишние капли жидкости. А По результатам взвешивания вычисляют коэффициент открытой пористости k ^ Упор. о_ (М3- МХ)/6Ж ^ М9-Мг П'° Уобр (Мв-М2+Мк)дж М3-М2+Мк ' ' где ]/обр — объем [образца; Упор. 0 — объем открытых пор образца; бж — плотность рабочей жидкости; (Мд—Мг) — масса жидкости в поровом пространстве образца; (М3—М2 + Мк) — масса объема жидкости, равного объему образца (масса жидкости в объеме об- образца). Рабочая жидкость, используемая для определения пористости, не должна вступать с веществом породы в химические взаимодейст- взаимодействия; вызывать набухания, отслаивания и деформации частиц по- породы; она должна иметь низкую упругость насыщающих паров и быть термостабильной. Если породы не содержат растворимой или набухающей в воде компоненты, то в качестве рабочей жидко- 14
сти рекомендуется использовать пластовую воду или ее модель. Если же для насыщения пород нельзя использовать воду, в каче- качестве рабочей жидкости выбирают очищенный от смол керосин, плотность которого периодически должна контролироваться. При определении коэффициента общей пористости kn исходят из очевидного соотношения Ум /т 1Я\ A.16) кобр "обр где у^ — объем твердой фазы образца (минералов). Так как Vo6v = Мс/бс; Vu = Мм/бм иМм^ Мс, то AЛ9) где Мс и Мы — массы соответственно сухой породы и ее твердой фазы; бс — объемная плотность породы в воздушно-сухом состоя- состоянии; бм — минералогическая плотность ее твердой фазы (скелета). Отсюда следует, что коэффициент общей пористости породы определяется отношением плотностей сухой породы и твердой фазы. Значение этих плотностей определяют весовым способом, предложенным Мелчером: бс — парафинированием образца с по- последующим гидростатическим взвешиванием, бм — пикнометри- чески (см. § 3). Зная коэффициенты общей и открытой пористости, оценивают величину закрытой пористости § 2. ВОДОНАСЫЩЕННОСТЬ И НЕФТЕГАЗОНАСЫЩЕННОСТЬ Содержание воды в горных породах называется их влаж- влажностью, а способность пород удерживать в себе то или иное количество воды в определенных условиях — влагоем- влагоемко с т ь ю. В естественных условиях поровое пространство пород-коллекто- пород-коллекторов может быть заполнено водой полностью или частично. В по- последнем случае оставшаяся его часть может быть заполнена нефтью или газом. Количественное содержание воды (нефти, газа) и ее состояние в породах существенно влияет на формирование ряда физических свойств и имеет большое значение при электрических, нейтронных и других геофизических методах исследования скважин. Водонасыщенность Вода, содержащаяся в породах, в зависимости от характера ее взаимодействия с твердыми частицами имеет различное состояние и подразделяется на две основные категории: связанную и свободную. 15
Содержание связанной воды в породах обусловлено химическими и физико-химическими (адсорбционными) процессами. В связи с этим различают химически связанную и адсорбционно связанную воду. К химически связанной относится вода кристаллиза- кристаллизационная и конституционн а я. Кристаллизацион- Кристаллизационная вода входит наряду с другими молекулами и ионами в кристал- кристаллическую решетку некоторых минералов в виде Н2О (гипс CaSO4-2H2O, эпсомит MgSO4-7H2O, мирабилит Na2SO4-10 H2O, сода Na2CO3 10H2O и др.) и, как правило, удаляется из них при температурах более 110°С. Исключение составляют сода и мира- мирабилит, из которых кристаллизационная вода удаляется при обычных температурах воздуха. Конституционная вода образуется при нагревании ряда мине- минералов из входящих в их кристаллическую решетку гщфоксильных ионов ОН~, Н+, Н3СИ~ и часто называется гидроксильной. Среди минералов, в которых синтезируется конституционная вода, можно отметить гидраргиллит А1 (ОНK, тальк Mg3 (OHJ (Si4O12), диаспор А1 (ОН) О, аллофан Al2O3-nSiO2, каолинит А12 (ОНLХ X(Si2O5), галлуазит Al2 (OHL(Si2O5)-H2O, монтмориллонит (AlMgJ(OHJ(Si4O10).H2O и др. В осадочных породах конституционная вода преимущественно образуется за счет минералов глин, особенно таких, как аллофан, каолинит и монтмориллонит. Синтез гидроксильной воды и ее уда- удаление из различных минералов обычно происходят при температу- температурах от ПО до 1300 °С. Химически связанная вода (кристаллиза- (кристаллизационная и конституционная) практически не оказывает влияния на физические свойства пород, за исключением нейтронных. Слои воды, находящиеся на разном расстоянии от поверхности твердых частиц или поровых каналов, отличаются по своему фи- физическому состоянию и химическому составу (рис. 3). Некоторая часть воды в пористых средах подвержена влиянию физико-хими- физико-химических (адсорбционных) сил, которые по своей природе являются электрическими и проявляются на поверхности гидрофильных ча- частиц. Эту часть воды принято называть адсорбционной. (Породы, твердая фаза которых смачивается водой, считаются гидрофильны ми, если твердая фаза не смачивается во- водой — гидрофобными.) Действие электрического поля вокруг твердой дисперсной ча- частицы приводит к образованию ориентированного и прочно адсор- адсорбированного слоя связанной воды. В зависимости от заряда твер- твердой поверхности молекулы воды ориентируются отрицательными или положительными зарядами к поверхности, образуя двойной электрический слой Гельмгольца (см. § 6). Адсорбированная вода подразделяется на прочно связанную и рыхло связанную. Прочно связанная вода удерживается на по- поверхности твердых частиц Ван-дер-Ваальсовыми силами притяже- притяжения, достигающими десятков тысяч атмосфер, что обусловливает 16
ее передвижение только при переходе в парообразное со- состояние. Рыхло связанная вода образуется на поверхности частиц породы поверх проч- прочно связанной воды. Молекулы этой воды удерживаются у твердой фазы пород менее прочно силами Ван-дер- Ваальса, силами поля ориен- ориентированных молекул прочно связанной воды и образуют шдратные оболочки обмен- обменных катионов. Содержание связанной во- воды в породе определяется минеральным составом и сте- степенью дисперсности (удель- (удельной поверхностью) ее твердой фазы. Величина последней в значительной мере зависит от глинистости. В связи с этим для более глинистых пород характерно и большее содержание связанной воды. К свободной (подвижной или извлекаемой) воде относят воду, которая не подвержена действию адсорбционных сил и способна сравнительно лег- легко передвигаться в порах и пустотах пород под дей- действием гравитационных сил или напора. Содержание свя- связанной воды в поровом про- пространстве пород количест- _ венно оценивается коэффи- коэффициентом связанной водонасы- щенности Ъ __ в св /т on «в. св — — у A 21) "пор где VB св — объем связанной воды, соответствующий доле объема пор, занятого связан- связанной водой; Fnop — общий объем пор породы. 000000000000 300000000000 о о о CS1 17
Связанная и условно подвижная вода обусловливает остаточ- остаточную водонасыщенность, характеризуемую коэффициентом kBm 0Ст остаточного водонасыщения. Количество подвижной (извлекаемой) воды в породах оцени- оценивается^ коэффициентом kB подв подвижной водонасыщенности К подв = Ув'П0ДВ = 1 - К ост, A.22) ^пор а общее содержание воды в поровом пространстве — коэффициентом kB водонасыщенности ? = kB, подв 4~ #в. ост == ^в. подв ~Ь «в. усл. подв ~Г «в. св> (^«23) * пор На практике часто содержание подвижной и условно подвижной воды рассматривается совместно и характеризуется коэффициентом kB Эф эффективной водонасыщенности V а ^в. эф — —¦— =1 — kB, св = ^в. подв + &в. усл. подв- A.24) »^пор Нефтегазонасыщенность В поровом пространстве пород-коллекторов наряду с водой мо- может содержаться нефть или газ. При этом количественное соотно- соотношение и характер распределения воды, нефти и газа в поровом пространстве могут быть различными. Распределение нефти (газа) и воды в порах зависит от того, гидрофильна или гидрофобна по- порода. Твердая фаза осадочных пород сложена в основном гидрофиль- гидрофильными минералами, что обусловливает распространение в природе преимущественно гидрофильных коллекторов нефти и газа. В общем случае A.26) >'пор мпор где FH, Vr и VB — доли объемов порового пространства, занятых соответственно нефтью, газом и водой; kHj kT и kB — коэффициенты нефте-, газо- и водонасыщенности. Нефть или газ, находящиеся в поровом пространстве породы, извлекаются лишь частично. В связи с этим различают коэффи- коэффициенты извлекаемого и остаточного нефтенасыщения (газонасыще- (газонасыщения), сумма которых равна kH (kT): «н '==z %. извл 4~ Ян. ост- A*27) В частном случае для гидрофильной нефтенасыщенной породы &н + &в — ^н. извл + ^н. ост + &в- подв + ^в. ост ^ 1- A.28) 18
Если поры нефтеносной породы не содержат свободноподвижной ВОДЫ, {К подв = °)> что чаЩе Всег0 Имеет месТ0> Т0 &н + &в == kn. извл + Ьн. ост+ k*. ост = 1. A.29) Определение коэффициентов водо-, нефте- и газонасыщенности Определение коэффициента водонасыщенности kB и нефтенаеы- щенности kH заключается в анализе образцов естественных пород на содержание в них воды и нефти. Коэффициент газонасыщенно- газонасыщенности можно вычислить из соотношений kv = 1 — (kB + ku) или kr z=z \ — kB (если порода не содержит нефти). Существует несколько прямых и косвенных способов определе- определения коэффициентов водо- и нефтенасыщенности: отгонки, критиче- критической температуры растворения, титрования, экстракции пента- ном, хлоридный и др. [10]. Из перечисленных способов наиболее широкое применение в практике получил способ прямой отгонки. Сущность этого способа заключается в том, что образец породы, отобранный при бурении скважины на промывочной жидкости с нефтяной основой и предохраненный от потерь нефти, т. е. с ес- естественным насыщением, помещают в бюксу и определяют его массу M2 = M1~M0i где Мг — масса бюксы с образцом; Мо — масса бюксы (измеряется отдельно). Затем с помощью специального аппарата из образца отгоняют нефть и воду. При этом воду улавливают в градуированной ло- ловушке, в которой определяется ее объем FB, а следовательно, и масса М3, так как плотность бв = 1. После отгонки образец сушат до постоянной массы М4 и по полученным результатам рассчиты- рассчитывают массу содержащейся в нем нефти Мц и ее объем VH' где бн — плотность нефти. Определив отдельно для исследуемой породы коэффициент об- общей пористости kn и плотность ее в сухом состоянии бс, рассчиты- рассчитывают: 1) объем пор Vnop = Mtkn/6c, (I.30) 2) коэффициент нефтенасыщенности 2Ь <'3'> *— г2Ь <'-3'> Упор Msfcn 3) коэффициент водонасыщенности k - У* У*Ь* л 321 В Vnop M4kn ' У ' ' 19
«==44 Рис. 4. Аппарат Закса для определения водонефте- насыщенности образцов пород 1 — кипятильник; 2, 10 — трубки; 3 — ловушка; 4 — o6pai- ный холодильник; 5 — стаканчик, 6,9—- пробки пришлифо- пришлифованные; 7 — крючок, 8 — рубашка теплоизолирующая; // — сифон; 12 — кран; 13 — трубка с охлаждающей водой; 14 — трубка для отвода нагревшейся воды, 15 — сифон стаканчи- стаканчика; 16 — дужка Если УПор ф Vn + VB» делают вывод о со- содержании в порах породы газа, объем кото- которого и рассчитывают коэффициент газонасыщенно- газонасыщенности Vr VTbc / ^ АЛ U # (L05t3) kr =¦ Аппарат для определения водонефтенасы- щенности (аппарат Закса или ВН-3) показан на рис. 4. В качестве нагревателя используется элек- электрическая печь ПЭ-3 или электроплитка закры- закрытого типа. Работу по определению коэффициентов водо- и нефтенасыщенности выполняют в следующем порядке. 1. Подготавливают к работе все составные части аппарата ВН-3. 2. Из середины нефтеводонасыщенного керна откалывают ку- кусочки породы, помещают их в бюксу известной массы Мо и опреде- определяют взвешиванием массу бюксы с породой Мг. Взвешивание про- производят на технических весах 1-го класса или аналитических ве- весах типа АДВ-200. 3. Устанавливают в вытяжном шкафу на асбестовом листе элек- электроплитку. 4. Заполняют кипятильник на 2/3 емкости растворителем — толуолом или спиртобензольной смесью (90 % бензола и 10 % спирта) и устанавливают его на электроплитку. Растворитель вы- выбирают в зависимости от типа нефти в образце. 5. Переносят породу в гильзу из фильтровальной бумаги, ко- которая помещается в стаканчик, подвешиваемый на крючок 7. 6. Закрывают кипятильник пробкой 9 и закрепляют рубашку 8 на штативе. 7. Соединяют холодильник с ловушкой и, подключив с помощью резинового шланга трубку 13 к водопроводу, устанавливают цир- циркуляцию охлаждающей воды. 8. Включают в сеть через латр электроплитку. Переключате- Переключателем латра регулируют нагрев плитки так, чтобы кипение раство- растворителя в кипятильнике было небурным. 20
Пары растворителя и воды образца породы поступают через трубку в холодильник, где они конденсируются; конденсат, пред- представляющий смесь растворителя и воды, стекает в ловушку. Вода, вследствие большей плотности, собирается в нижней части ловушки, а "растворитель по мере заполнения ловушки до уровня сифона 11 сливается в стаканчик (на образец) и заполняет его до тех пор, пока не сработает сифон 15. При срабатывании сифона 15 растворитель с растворенной в нем нефтью из образца сливается в кипятильник. Экстрагирование (отгонка) считается законченным, если объем воды в ловушке не увеличивается, а растворитель над ней и сли- сливающийся из стаканчика через сифон 15 в кипятильник становится прозрачным (без видимой окраски). 9. По окончании экстрагирования выключают плитку, дают аппарату остыть и определяют объем VB = М3/бв воды, скопив- скопившейся в ловушке (бв = 1). 10. Разбирают аппарат, переносят образец породы в бюксу мас- массой Мо и сушат его в сушильном шкафу при температуре 105 °С до постоянной массы (М4 + Мо). 11. Определяют значения: М2 = Мг—Мо\ М4 = (Af4 + + Мо) — Мо\ Мп = М 2 — (М8 + М4) и при известных коэффи- коэффициенте пористости kny плотности бс сухой породы и плотности бн нефти по формулам A.30) — A.33) рассчитывают искомые вели- величины &в, kn и кг. § 3. ПЛОТНОСТЬ Плотность б — физическая величина, которая для одно- однородного вещества определяется его массой 1 в единице объема. Для практических целей часто используют относительную плотность. Для жидких и твердых веществ она устанавливается по отношению к плотности дистиллированной воды при 4 °С, для газов — по от- отношению к плотности сухого воздуха при нормальных условиях (р - 101325 Па, Т - 0). Средняя плотность тел (в кг/м3) 6 = M/V. A.34) Плотность достаточно тесно связана с рядом физико-химиче- физико-химических свойств горных пород и оказывает влияние на показания радиоактивных, акустических и других геофизических методов исследования скважин. Удельный вес равен отношению веса тела к его объему и может быть определен как произведение плотности б на ускорение свободного падения g. Следовательно, удельный вес является фи- физико-химической характеристикой вещества, так как зависит от значений g. 1 Масса — величина, характеризующая количество вещества в теле, и равная весу тела, деленному на ускорение свободного падения g. 21
Плотность твердой, жидкой и газообразной фаз пород Плотность твердой фазы 8М пород зависит от плотностей со- составляющих ее минералов и их соотношения в единице объема этой фазы. Минералы, наиболее часто встречающиеся в осадочных по- породах, по плотности могут быть условно подразделены на три группы [2J: 1) глинистые с 8М A,5—2,6)-103 кг/м3, некоторые сульфаты и хлориды [гипс B,32-103 кг/м3), галит B,20-103 кг/м3) и др. ]; 2) ос- основные породообразующие с 8М от B,6 + 2,65)-103 кг/м3 (полевые шпаты, кварц) до B,71— 2,85)-103 кг/м3 (каль- (кальцит, доломит), реже до 2,95-103 кг/м3 (ангидрид); 3) сопутствующие тяжелые с бм, C,5—5)-103 кг/м3 и более. Твердая фаза пород может быть мономинеральной и по- лиминеральной. В первом случае ее плотность совпа- совпадает с плотностью породооб- породообразующего минерала, во вто- втором определяется величиной 0,30 Рис. 5.- Зависимость плотности бр рас- растворов NaCl от их объемной концентра- концентрации Су и температуры Т [2] средней взвешенной плотно- плотности минералов и рассчиты- рассчитывается по формуле i=n A.35) где 8Mt-— плотность /-го минерала; VMn — объем, занимаемый /-м минералом в единице объема твердой фазы; п — число минералов, составляющих твердую фазу породы. Жидкая фаза пород обычно бывает представлена пластовой во- водой или пластовой водой и нефтью. Плотность пластовых вод в сво- свободном и рыхло связанном состояниях зависит в основном от хи- химического состава, минерализации и температуры (рис. 5) и изме- изменяется в пределах @,95—1,2)-103 кг/м3 [2]. Природные нефти характеризуются незначительным пределом изменения плотности. В нормальных условиях (р = 0,1 МПа, Т = = 20 °С) в зависимости от химического состава 8Н изменяется в пределах @,7—1,06) • 103 кг/м3. Повышенная плотность обычно свойственна окисленным нефтям с высоким содержанием асфаль- тенов и смол. В пластовых условиях 8Н зависит от температуры, давления и газового фактора. Различие плотностей нефти в пласто- пластовых и поверхностных условиях учитывается с помощью объемного коэффициента. 22
Если жидкая фаза состоит из воды и нефти, то ее плотность рас- рассчитывается как средняя взвешенная величина 8« = SBVB-f ShVh, О-36) Где 6В и бн — плотности воды и нефти; VB и 1/н — занимаемые во- водой и нефтью объемы в единице объема жидкой фазы. Плотность природных газов бг в нормальных условиях зависит от их химического состава и обычно определяется отношением к плотности воздуха при тех же условиях. В отличие от твердой и жидкой фаз бг более существенно зависит от температуры и давле- давления. Плотность пород Плотность пород бп зависит от содержания в единице объема породы твердой Мш жидкой Мж и газообразной Мг фаз и соответст- соответственно их плотностей бм, бж, бг. Плотности отдельных фаз опреде- определяются следующим образом: Ум мт Ут A.37) 2,75 2,50 2,25 2,00 /,75 10 Рис. 6. Зависимость плотности бп от коэффициента пористости kn известняков верхнего мела Восточного Предкавказья Плотность породы Уп A.38) 23
или бп = где 8Mi, бж* и 8rt—плотности минералов твердой фазы породы, жидкостей и газов, заполняющих ее поровое пространство; VKu Уж h Vh — их объемы в еди- единичном элементе породы. Если плотности твердой фазы и пластовой жидкости постоян- постоянны, то при полном насыщении Vs/ ////// бп является функцией пористо- пористости породы (рис. 6): A.40) Это позволяет использовать за- зависимость 8n — f (kn) для опре- Рис. 7. Номограмма для определе- определения пористости неглинистых пород деления плотности пород по ве- по данным б„, бм личине коэффициента пористости и наоборот. На рис. 7 приведена номограмма, построенная в соответ- соответствии с уравнением A.40) при условии, что бж = бв = 1. Эта номограмма позволяет определять пористость неглинистых водона- сыщенных пород по величине их плотности. В глинистых породах твердая фаза состоит из основных породо- породообразующих минералов, составля- составляющих скелет породы с плотно- плотностью бм, и глинистого цемента с плотностью бгл, заполняющего частично поровое пространство между зернами скелета. В этом случае A.41) где krjl — объемная глинистость. Откуда ¦ о г 4 е у _ ^м Оп у _ Ом Огл Рис. 8. Зависимость плотности пород от глубины их залегания П. 2] J2 & ? 51 Ом — Ож °м — Ож A.42) С увеличением глубины зале- залегания плотность пород, как пра- правило, возрастает, что связано с их уплотнением и, как следствие этого, уменьшением пористости под давлением вышележащих толщ. Зависимость плотности пород от глубины залегания, по данным / — известняки; 2 — песчаники; 3 — глины 24
работ [2, 3], в среднем описывается эмпирическим уравнением вида 6пЯ = бм - (*„ ~ бв) *пН=0 «*Р ( - *ТН) A.43) и графически представлена на рис. 8. Здесь knHQ — пористость породы на поверхности; Н — глубина, на которой определяется плотность; ат — температурный коэффициент, изменяющийся от 0,2- Ю~3 м (породы с жестким скелетом и низкой глинистостью) до 0,7-Ю~3 м (породы высокоглинистые). Наибольшее уплотне- уплотнение характерно для глин и сильноглинистых разностей пород. Хо- Хорошо отсортированные чистые песчаники уплотняются в меньшей степени. Определение плотности В лабораторных условиях обычно определяют не плотность пород в естественном залегании, а плотности сухой породы бс и ее твердой фазы бм, пластовой воды бв, нефти бн и газа 6Г. При из- известных плотностях составляющих фаз и объемном содержании их в породе, характеризуемом коэффициентами пористости, водо-, яефте- и газонасыщенности (kBy ka, 4), плотность пород с естествен- естественным насыщением рассчитывают по формуле A.39). Рассмотрим способы определения бс и бм. Для оценки бс доста- достаточно знать массу сухого образца породы Мг и его внешний объем Уобр, тогда бс = MJVoev Внешний объем образца можно опре- определить способом парафинирования с последующим гидростатиче- гидростатическим взвешиванием или вычислением, используя данные, получен- полученные при оценке открытой пористости способом насыщения. Способ парафинирования заключается в следующем. 1. Экстрагированный образец породы помещают в бюксу и су- сушат в сушильном шкафу при температуре 105 °С до постоянной массы. 2. Определяют массу сухого образцаjAfj, взвешиванием в воз- воздухе на технических (или аналитических) весах 1-го класса. 3. Образец обвязывают тонкой гибкой проволочкой длиной 20—30 см и покрывают его слоем парафина известной плотности толщиной 0,5—1 мм. Парафинирование производят несколькими повторными погружениями образца в расплавленный парафин, имеющий температуру не выше 60 °С. При большей температуре парафин может проникнуть глубоко в поры образца и исказить результаты определения объема. 4. Находят массу парафинированного образца взвешиванием в воздухе М2- Массу проволочки устанавливают отдельно и вычи- вычитают из результатов взвешивания. 5. Определяют массу парафинированного образца гидростати- гидростатическим взвешиванием Л43. Для этого образец с помощью проволочки подвешивают к серьге коромысла весов и погружают в стакан с ди- дистиллированной водой, устанавливаемый на мостике или кольце штатива. Массу проволочки определяют гидростатическим взвеши- взвешиванием отдельно и вычитают из результатов взвешивания. 25
6. Вынимают образец, обтирают его фильтровальной бумагой, повторно находят массу парафинированного образца взвешиванием в воздухе М'2 и убеждаются, что вода в него не проникла. При МГ2—М2 > 0,02 г результаты взвешивания считают непригод- непригодными. Величину плотности сухого образца вычисляют по формуле где бв — плотность дистиллированной воды при температуре из- измерений; бпар — плотность парафина, наиболее часто равная, @,87—0,9). 103 кг/м3. Если плотность парафина неизвестна, то ее устанавливают пикнометрическим способом или способом гидростатического взве- взвешивания на аналитических весах. Из данных, получаемых при определении пористости способом насыщения, очевидно, что (Ms—Мг) есть масса рабочей жидкости в порах, а (М3—М2) — масса жидкости в объеме образца (см. § 1). Следовательно, зная плотность рабочей жидкости легко вы- вычислить = (Afа - М2)/6Ж Vu - Voep - FnoP = (Mt - М2)/6Ж. A.45) Тогда плотность сухой породы 6с = —; —~у A.46) м3 — м2 а кажущаяся плотность твердой фазы (кажущаяся минералогиче- минералогическая плотность) *п.о«1 г^-. A.48) Ом. к Кроме того, можно записать Откуда бс = A~?л.оMм.к. A.49) Плотность твердой фазы (минералогическую плотность) опреде- определяют пикнометрическим способом. Для этого в пикнометр вводят исследуемое твердое вещество, затем доливают до риски жидкость известной плотности и устанавливают взвешиванием на техниче- технических весах 1-го класса его общую массу. Для измерения бм обычно используют часть образца сухой породы, оставшуюся после определения величины бс. Минималь- Минимальная масса ее должна составлять 10—15 г. В качестве рабочей жид- 26
кости применяют дистиллированную воду; если же порода содер- содержит растворимые или набухающие в воде компоненты, применяют очищенный керосин. _ Рассмотрим последовательность, в которой^выполняются из- измерения. 1. Экстрагированный и высушенный до постоянной массы об- образец измельчают в ступке до порошкообразного состояния. 2. Определяют массу MQ сухого пикнометра и засыпают в него 10—15 г породы. Рекомендуется использовать пикнометры емкостью 50 см3. 3. Определяют массу Мг сухого пикнометра с породой. 4. Заливают в пикнометр до половины объема рабочую жид- жидкость (при этом ее уровень должен быть выше породы примерно на 1—1,5 см), слегка взбалтывают и помещают в вакуумную уста- установку для удаления воздуха. Вакуумирование производят до пре- прекращения выделения воздуха. Затем вакуум снимают, доливают в пикнометр жидкость и вновь вакуумируют в течение 30 мин. По- После снятия вакуума пикнометр извлекают и термостатируют в во- водяной бане при 20 °С и доводят в нем уровень жидкости до риски, тщательно обтирают фильтровальной бумагой снаружи и опреде- определяют его массу М2. 5. Содержимое пикнометра сливают, заполняют его рабочей жидкостью, термостатируют в водяной бане при 20 °С и вновь до- доводят уровень жидкости до риски, обтирают снаружи фильтро- фильтровальной бумагой и определяют массу М3. 6. Вычисляют плотность твердой фазы по формуле Щ? бж 20 — плотность рабочей жидкости при температуре 20 °С. § 4. ПРОНИЦАЕМОСТЬ Свойство пород пропускать через себя жидкость, газы и их смеси при перепаде давлений называется проницаемостью. Проницаемость является мерой фильтрационной проводимости породы. Ее подразделяют на физическую (абсолютную), фазовую (эффективную) и относительную. Физическая проницаемость Физическая проницаемость соответствует про- проницаемости породы при фильтрации через нее однородной жидко- жидкости или газа, химически инертных по отношению к твердой фазе, н количественно оценивается коэффициентом физической проницае- проницаемости knp. В уравнении Дарси knp является коэффициентом про- пропорциональности между скоростью фильтрации v$ однородной 27
жидкости (газа) и градиентом давления Ap/L: Уф=йпр___=_, -(Ш) где Ар — перепад давления (в Па); L — длина пористой среды (в м); \i — динамическая вязкость жидкости (газа) [в Па-с]; Q — объемный расход жидкости (газа) в единицу времени (в м3/с) через сечение F (в м2) пористой среды. Отсюда коэффициент проницаемости (в м2) За единицу проницаемости значением в 1 м2 принимается про- проницаемость такой пористой среды, при фильтрации через образец которой площадью 1 м2, длиной 1 м и перепаде давления 1 Па рас- расход жидкости вязкости 1 Па-с составляет 1 м3/с. Физически эта единица измерения проницаемости характеризует величину пло- площади сечения каналов пористой среды, по которым происходит фильтрация.1 Горные породы условно подразделяются на проницаемые, (с knp > 10,2-10~15 м2), полупроницаемые [с knp = @,1 — 10,2)- •10~15 м2] и практически непроницаемые (с kup < 0,Ь105 м2). Физическая проницаемость коллекторов колеблется в очень боль- больших пределах: от 0,Ы0~15 до 3-10~12 м2 и более. Широкое рас- распространение имеют коллекторы с проницаемостью @,2—1,02) • • Ю-12 м2. Основным фактором, влияющим на коэффициент проницаемо- проницаемости пород, является структруа их порового пространства, характе- характеризуемая формой и размером пор, извилистостью и удельной по- поверхностью каналов фильтрации. Эти параметры определяют объем фильтрующего агента, траекторию его отдельных струй и силы поверхностного взаимодействия, препятствующие фильтрации. Теоретически согласно уравнению Козени-Кармана где &п. д — динамическая пористость образца породы, доли еди- единицы; S$ — удельная поверхность каналов фильтрации (в м2/м3); Тт — гидравлическая извилистость каналов фильтрации, равная отношению средней статистической длины поровых каналов LK к длине образца породы L; / — коэффициент, учитывающий форму сечения пор и изменяющийся для гранулярных коллекторов от 2 до 3 при наиболее частом значении 2,5. 1 Практической единицей измерения проницаемости является Дарси, равная 1,02-Ю-12 ма « 1 мкм2. Величина, равная 0,001 Д, называется мил- лидарси (мД). 28
Качественно зависимость проницаемости от пористости и струк- структуры порового пространства зернистых пород может быть показана на примерах статистических связей. По среднестатистическим дан- данным для многих нефтегазоносных коллекторов большей пористости обычно соответствует и большая проницаемость (рис. 9). Структура порового пространства обломочных пород до некоторой степени может быть охарактеризована удельной поверхностью зерен Syfl, их эффективным или медианным диаметром, коэффициентом от- отсортированное™, средним радиусом капиллярных пор и извили- 50 100 500 knp'10~™nz Рис. 9. Соотношение между проницаемостью &пр и открытой пористостью кп. о (песчано-алеврито-глинистые породы-коллекторы нижнего мела При- кумского нефтегазоносного района Ставрополья) стостью поровых каналов. Между коэффициентом проницаемости и указанными параметрами обычно наблюдаются статистические корреляционные связи, часто используемые совместно с зависи- зависимостью knp = f (kn) при интерпретации результатов геофизиче- геофизических исследований скважин (рис. 10). Разброс точек на графиках объясняется зависимостью knp от нескольких из перечисленных параметров, значение которых в каждом отдельном случае раз- различно. Согласно работам В. М. Добрынина и других исследовате- исследователей, с ростом эффективного давления (увеличением глубины за- залегания) проницаемость пород k^3^ уменьшается в соответствии с приближенным уравнением где Г 4Рэф) Т L kn J з+ас kl^' «k^ I — 1 , A.54) — структурный показатель, связывающий среднюю ста- 29
тистическую длину поровых каналов с их радиусом U = сгг ст. Здесь с — постоянная величина, изменяющаяся от (— 0,24) до <-^8> [3L На рис. 11 даны кривые относительного изменения коэффици- коэффициента проницаемости от глубины залегания (а) и эффективного дав- давления (б). Анализ многочисленных экспериментальных исследова- Рис. 10. Зависимость коэффициента проницаемости от структурных пара- параметров пород Зависимость &пр от а- удельной поверхности 5уд песчаников и алевритов Грозненского района (по А М Нечаю), б— медианного диаметра Md песчаников (по М. К. Калинко), в — остаточной водонасыщенности kB o карбонатных пород Куйбышевского Поволжья о (по данным К И. Багринцевой); г — связь отношения &п/?пр с УДельн°й поверхностью песчано-алевритовых пород Грозненского района (по А. М. Нечаю) ний показывает, что проницаемость пород, особенно глинистых, залегающих на больших глубинах, может оказаться почти в 2 раза меньше ее значения, определенного на поверхности. Проницаемость трещиноватой породы в общем случае зависит от геометрии систем трещин и направления фильтрации. Для приближенной оценки трещинной пористости и проницаемости используются шюскопа- ралдельные шлифы увеличенных размеров. Метод шлифов, хотя 30.
Рис. 11. Относительное изменение коэффи- коэффициента Л#пр/&пр проницаемости песчано-гли- нистых пород в зависимости от глубины зале- залегания Н (а) и эффективного давления рЭф (б) а — среднеотсортированные песчано-глинистые поро- породы (по В. М. Добрынину); б — песчано-алеврито- глинистые породы нижнего мела Ставрополья; А&Пр — уменьшение проницаемости под воздействием эффек- эффективного напряжения и дает приемлемые результаты, не находит широкого применения, так как требует исследования большого числа образцов. В связи с этим для оценки параметров трещиноватых пород все чаще ис- используются геофизические и гидродинамические методы. Фазовая и относительная проницаемости В нефтегазонасыщенных породах-коллекторах одновременно присутствуют две или три фазы (нефть—вода, газ—вода, газ—нефть—вода). Способность пород, насыщенных водонефтега- зовыми смесями, проводить отдельно нефть, газ, воду называют фазовой (эффективной) проницаемостью. Последняя характеризуется коэффициентом фазовой проницае- проницаемости, рассчитываемым для каждого компонента смеси по формулам ^пр. н = «пр г — FAp FAp A.55) A.56) A.57) где QB, QH, Qr, ^iB, }xHj jj,r — соответственно расходы и вязкости воды, нефти и газа. Отношение фазовой проницаемости к физической называют от- относительной проницаемостью k'up. „ = knp, B/knp; &пр. н = knp,Jkup; ^пр.г = kup.T/kupt которую выражают безраз- 31
мерной величиной в долях единицы или в процентах. Фазовая и относительная проницаемости зависят от нефте-, газо- и водона- сыщенности порового пространства породы, а также от физико-хи- Рис. 12. Зависимости коэффициентов от- относительной проницаемости от коэффици- коэффициентов нефтеводонасыщенности для песчано- глинистых пород 0,2 - Рис. 13. Области распространения одно-, двух- и трехфазного потока в пористой среде (по Леверетту) Вода W0°/o 90 80 70 60 60 W 30 20 /О О Водой асыщенногть) °/с мических свойств пористых сред и компонентов, насыщающих их смесей. Если часть пор занята какой-либо фазой, то совершенно очевидно, что проницаемость породы для другой фазы становится меньше. В реальных условиях, в зависимости от объемного соот- соотношения в поровом пространстве породы воды, нефти, газа, воз- возможны одно-, двух- и трехфазные потоки. При этом расход много- многофазных смесей всегда меньше расхода однородной фазы. На рис. 12 показаны зависимости коэффициента относительной проницаемости для двухфазной смеси (вода—нефть) от величины 32
коэффициентов нефте- и водонасыщенности. Участки диаграмм* в пределах которых относительная проницаемость практически становится равной нулю, определяют области фильтрации только одного из компонентов смеси. Положение этих областей и харак- характер кривых &пР.в = / (&в); ?пР н = / (?в); &пР г = / (kB) в опреде- определенной мере зависят от структуры порового пространства породы. При одновременном насыщении пород газом, нефтью и водой области одно-, двух- и трехфазных потоков изображают на треу- треугольной диаграмме, как это показано на рис. 13. На этой диаграмме нанесены кривые, соединяющие точки с одинаковым содержанием E °Л) соответствующего компонента смеси воды A), нефти B) и газа C) в потоке. Вершины треугольника соответствуют 100 %-ому насыщению породы одной из фаз, а его стороны, противолежащие соответствующим вершинам,— нулевому насыщению этой фазой. Из диаграммы (см. рис. 13), полученной для несцементированных лесков, видно, что при газонасыщенности среды меньше 10 % и нефтенасыщенности ниже 23 % в потоке практически будет со- содержаться только вода. Область трехфазного потока обычно мала и расположена в пределах насыщенности песка нефтью от 23 до 50 °'о, водой от 33 до 64 %, газом от 14 до 30 %. Для других типов пород-коллекторов эти пределы насыщенности, ограничивающие области одно- двух- и трехфазных потоков, могут быть несколько, отличными. Определение проницаемости пород Определение проницаемости пород в лабораторных условиях производят на специальных установках по газу (воздуху) или жидкости. Для одних и тех же пород проницаемость по жидкости несколько меньше, чем по газу, что связано с разбуханием глинистых частиц и образованием адсорбционных слоев воды (нефти) на поверхности каналов фильтрации. В связи с этим абсолютную (физическую) проницаемость пород чаще принято устанавливать при фильтра- фильтрации через них газа (обычно азота или воздуха). Для этого сущест- существует ряд установок, одна из которых (ГК-5) показана на рис. 14. Определение проницаемости выполняют в следующем порядке. 1. Изготавливают образцы пород цилиндрической формы диа- диаметров 20—30 мм, длиной 20—40 мм; экстрагируют их от нефти й солей в аппарате Сокслета и сушат до постоянной массы при температуре 105 °С. 2. Проверяют герметичность всей системы и работу отдельных узлов установки. 3. Измеряют штангенциркулем диаметр и длину образца, уста- устанавливают и надежно зажимают его в струбционном зажиме. 4. К нижнему штуцеру струбционного зажима через хлоркаль- цневую трубку и манометр на входе подают осушенный от влаги азот (воздух). 2 Заказ № 789 33
5. Фиксируют по соответствующим манометрам установив- установившиеся давления на входе рг и выходе р2 из образца и с помощью расходомера (например, реометра) измеряют расход газа в еди~ ницу времени. 6. Находят температуру фильтрующегося газа с помощью тер- термометра, установленного в хлоркальциевой трубке, и по табл. 1 Рис. 14. Схема расположения узлов в установке ГК-5 1 — манометр на входе; 2 — фильтр, 3 — редуктор; 4 — хлоркальциевая трубка; 5 — термометр, 6 — образцовый манометр; 7 — струбщюнный зажим; 8 — реометр определяют его вязкость при этой температуре; с барометра сни- снимают значения атмосферного давления рб. Таблица 1 Вязкость азота и воздуха в зависимости от температуры Температу- ^ [pa, °C 15 16 17 18 19 20 Вязкость азота 0,01741 0,01746 0,01751 0,01756 0,01761 0,01766 , МПа с воздуха 0,01783 0,01788 0,01793 0,01798 0,01803 0,01812 Температу- Температура, °С 21 22 23 24 25 0 0 0 0 0 Вязкость, азота ,01768 ,01771 ,01773 ,01776 ,01778 МПа с воздуха 0,01818 0,01822 0,01829 0,01834 0,01840 Коэффициент проницаемости вычисляют с учетом давления газа» его расхода и температуры. Для этого на основании уравнения состояния газов расечиты* вают: а) действительный расход газа Q по формуле где Q' — расход газа, определенный по показаниям реометра; 34
g — плотность азота по воздуху при температуре Т", при V = J= 20 "С и плотности азота по воздуху, равной 0,9836: Q « 0,058Q' У B73 -h Г)/Рб, A-59) б) средний расход газа по длине образца при атмосферном дав» лении рв " температуре измерений Т по формуле Qcp = 2QP6i(Pi + Рг + 2Рб) = Qpd{P6 + АР/2)- A-60) Тогда при фильтрации газа формула A.52) принимает следующий вид: Pi -г Р2 + 2Рб (Pi — Pa) ^ В случае определения расхода газа газометром р2 = 0 и При детальных исследованиях керна иногда необходимо опреде- лягь проницаемость образцов пород, моделируя радиальный по- гок жидкости или газа, т. е. как бы воспроизводя условия их при- притока в скважину. Для этого изготавливают образцы пород цилин- цилиндрической формы с осевым отверстием — «скважиной» и создают в нем с помощью специальных кернодержателей фильтрацию жид- когти или газа в радиальном направлении от наружной поверхно- поверхности к внутренней. Коэффициент проницаемости в этом случае рассчитывают: а) при фильтрации жидкости по формуле ^Пр Ж > /т /^QV 2jt/i (pH — рв) (l.Do; при фильтрации газа по формуле где [1Ж и fxr — вязкости жидкости и газа; Q>K — расход жидкости; Qr* Qr. cP — расходы газа при атмосферном и среднем давлениях в образце; гн и гв — наружный и внутренний радиусы кольца; Рп и ръ — давление у наружной и внутренней поверхностей кольце- кольцевого образца; h — высота образца. Кроме этого, специальные исследования проницаемости пород могут выполняться с помощью моделирования условий, близких к пластовым. Одной из установок для таких измерений является установка УИПК, последние варианты которой позволяют нахо- находить коэффициент проницаемости по газу, жидкости и их смесям при температуре до 90 °С, пластовом и горном давлениях соответст- соответственно до 60 и 100 МПа. * 35
Глава II ЭЛЕКТРИЧЕСКИЕ, РАДИОАКТИВНЫЕ, АКУСТИЧЕСКИЕ И ДРУГИЕ СВОЙСТВА ПОРОД § 5. УДЕЛЬНОЕ ЭЛЕКТРИЧЕСКОЕ СОПРОТИВЛЕНИЕ Свойство горных пород проводить электрический ток характе- характеризуется их удельной электропроводностью о или величиной, обратной ей,— удельным электриче- электрическим сопротивлением где R — полное электрическое сопротивление образца породы (в Ом); S и L — площадь поперечного сечения (в м2) и длина (в м) образца. Из формулы (II. 1) следует, что величина р измеряется в оммет- омметрах. Удельное электрическое г сопротивление в 1 Ом-м равно пол- полному сопротивлению в Ом 1 м3 породы с основанием 1 м2 и высотой 1 м, измеренному перпендикулярно к плоскости куба. Удельное сопротивление горной породы определяется удельным сопротивле- сопротивлением твердой фазы, жидкостей и газов, насыщающих поровое про- пространство, их объемным соотношением, характером распределения в породе и температурой. Удельное электрическое сопротивление твердой фазы, водных растворов солей, нефтей и газов Удельное сопротивление твердой фазы пород зависит от ее ми- минералогического состава и температуры. Минералы весьма разно- разнообразны по своему удельному сопротивлению, которое изменяется в широких пределах A0~6—1015 Ом-м). Однако основные минералы, образующие скелетную часть твердой фазы осадочных пород, {кварц, полевые шпаты, кальцит, слюды и др.), характеризуются удельным сопротивлением от 1010 до 1015 Ом-м и практически не проводят электрического тока. Присутствие в скелетной части твердой фазы полупроводящих минералов (графит, пирит, магнетит и др.) снижает ее удельное сопротивление в зависимости от их количественного содержания и характера распределения. В природных условиях содержание в осадочных породах минералов повышенной электропроводности невелико и, как правило, не приводит к существенному изменению удельного сопротивления скелетной части. В связи с этим ее р принято считать практически бесконечным. Пластовые поровые воды представляют собой растворы солей (электролиты) и относятся к классу ионных проводников. Их удель- 1 В дальнейшем для сокращения изложения слово «электрическое» опу- опускается. 36
яое сопротивление зависит от химического состава растворенных солей, концентрации и температуры. С увеличением общей концен- концентрация солей удельное сопротивление растворов уменьшается. При этом электропроводность раствора приблизительно равна сумме электропроводностей, обусловленных каждой из солей в отдельно- отдельности. Зависимость удельного сопротивления рв растворов от хими- химического состава растворенных солей и концентрации ионов выра- выражается формулой 10 Рв^———— — , A12) L,{CaiVtfai+CKiuifKi) где С& if CK i — концентрации (в молях) t-ro аниона и i-ro катиона в растворе; vu ui — электролитические подвижности анионов и 2 1 0 ¦f _N03 "Br I " Ctt "к7" "щ нсо ¦—-'"«. 3 ——— ¦— " — ' — — ^\ — — / / viaaCl ч Ca К\ V \ 0,01 0,05 0,1 0,2 0,5 1 2 5 W 20 50 Ссум,г/л Рис. 15. Номограмма для пересчета суммарной концентрации раствора Ссум в эквивалентную по NaCl с помощью множителя а катионов при бесконечном разбавлении (в Ом-см2/моль); fat, fKi — коэффициенты электропроводности анионов и катионов. В пластовых водах обычно преобладает содержание хлористого натрия, которое достигает 70—95 %. Если содержание других со- солей не превышает 10 %, то для практических целей удельное со- сопротивление пластовой воды оценивается по общей концентрации, приравненной к концентрации NaCl. Если же содержание других солей превышает 10 %, то определение удельного сопротивления воды необходимо проводить с учетом солевого состава. Для практических расчетов используется график (рис. 15), позволяющий по результатам химических анализов пластовых вод находить суммарную концентрацию раствора, эквивалецтную кон- 37
фоог-о а.
центрации NaCl. Рассчитав эквивалентное содержание NaCl, с по- помощью множителей а (см. рис. 15) по зависимости рв = / (CNaC1) при заданной температуре определяется искомая величина рв. Пример. Согласно химическому анализу пластовая вода при 29 °С содержит Са 0,54 г/л, SO^ 1,3 г/л и NaCl 15,5 г/л. Суммарная концентра- концентрация солей равна 17,34 г/л. Множителем для Са является 1,21; для SO4 — п 6э По полученным данным находим приближенное эквивалентное содер- содержание NaCl: 0,54-1,21 + 1,3-0,65+ 15,5= 0,66+ 0,84+ 15,5= 17 г/л. Согласно номограмме (см. рис. 16) при 29 °С рв = 0,35 Ом-м. С увеличением температуры удельное сопротивление водных растворов закономерно уменьшается: Р ^ " (IL3> где ат и |3Г—температурные коэффициенты, равные для раство- растворов NaCl соответственно 216-10~4 и 8-10~6; Т и То — температуры (в С) соответственно раствора и стандартная (равная 20 или 25 °С); 1> — удельное сопротивление раствора при Го; рвГ — удельное сопротивление раствора при Т. На рис. 16 представлена номограмма зависимости удельного сопротивления растворов NaCl от их концентрации и температуры, наиболее часто используемая для определения величины рв. За- Зависимость удельного сопротивления пород от температуры анало- аналогична зависимости для водных растворов: с увеличением темпера- температуры удельное сопротивление пород также уменьшается. Удельное сопротивление природных нефтей и газов во много раз превосходит удельное сопротивление пластовых вод и соизме- соизмеримо с р скелетной части твердой фазы пород. Практически элек- электропроводность нефтей и газов принимается равной нулю. Удельное сопротивление неглинистых пород Горные породы проводят электрический ток в основном за счет наличия в их поровом пространстве водных растворов солей. В связи с этом удельное сопротивление неглинистой породы рвп грануляр- гранулярного строения, поры которой полностью насыщены водой, зависит от рв этой воды, ее количества и характера распределения в породе, определяемых соответственно коэффициентом пористости kn и структурой порового пространства. Для исключения влияния удельного сопротивления пластовой воды вместо рвп для водонасыщенных пород принято рассматривать величину Р-Рвп/Рв, (Н.4) называемую относительным сопротивлением. Для чистых (неглинистых) пород Р не зависит от удельного сопро- сопротивления насыщающих вод, а связано с величиной пористости п структурой порового пространства. В связи с этим его называют также и параметром пористости. 39
Для пород с однородной структурой, сложенных частицами пра- правильной формы, связь между относительным сопротивлением и по- пористостью может быть рассчитана аналитически. На рис. 17 пока* заны примеры простейших форм порового пространства породы, для которых зависимости Р от kn имеют следующий вид: для а р — п для б для в для г Р = ¦ 1- l-k \2/3 A15) A16) Рис. 17. Примеры простейших форм порового пространства породы Для пород, содержащих частицы различных размеров, зависи- зависимость между относительным сопротивлением и пористостью выра- выражается в более сложной форме: " L 2 A — v) J lg A - v) 2(l-v)J ' <IL9> где v — коэффициент заполнения пор основного скелета породы более мелкими частицами. Анализ теоретических расчетов, выполненных для идеальных пород, подтверждает, что относительное сопротивление зависит не только от коэффициента пористости, но и от структуры порового пространства, обусловливающей характер распределения воды в породе. Расчеты и экспериментальные исследования показывают что зависимость Р = f (ku) имеет обратный степенной характер и может быть аппроксимирована в диапазоне средних значений kn (от 3—5 до 30—40 %) уравнением общего вида где а — некоторая постоянная; т — структурный показатель, за- зависящий от формы поровых каналов. 40
Численные значения а и т для различных пород могут соот- соответственно изменяться от 0,4 до 1,0 и от 1,3 до 2,4. Эти параметры определяются экспериментально на образцах породы. На практике широкое применение находят и частные виды этого уравнения (рис. 18): ^ = l/C (H.11) где структурный показатель т изменяется от 1,3 (для слабосцемен- хированных пород) до 2,3 (для крепких хорошо сцементированных пород); для среднесцементиро- ванных песчаников часто ис- р 'JO }С :и w 5 2 1 . \\\ - 2,2 1 \ \Л\\\ пользуют выражение A1.12) 0 го 50кг% Рис 18. Зависимости относительного сопротивления Р от коэффициента по- ристосли ku, рассчитанные по формуле A1.11) Зависимость, показанная пунктирной линией, соответствует уравнению Хамбла Р = Рис, 19. Схематическая модель анизотропной породы при распространении тока перпендикулярно к наслоению (а) и параллельно ему (б) /<!t и /?2~~ сопротивления частей породы с удельными сопротивлениями р^ и Рр при рас- распространении тока / перпендикулярно к наслоению; R3 н Ri — то же, при распростране- вии тока параллельно наслоению В слоистых породах удельное сопротивление в направлении, параллельном наслоению pBn[j, отличается от его значения, из- измеренного в направлении, перпендикулярном к наслоению pBnL. Такие породы называются анизотропными по удельному сопротив- сопротивлению. Явление анизотропности более заметно проявляется в гли- 41
нистых сланцах, мергелях, тонкочередующихся песчано-глини- стых и других породах. На примере единичного элемента породы (рис. 19), состоящего' из прослоев удельных сопротивлений рр и ps, легко установить, что (Pp-PsJ •=Л+- A1.13) Отсюда следует, что удельное сопротивление породы, измерен- измеренное перпендикулярно к наслоению, всегда больше удельного со- сопротивления этой породы в направле- направлении, параллельном наслоению, т. е. Рвп± ^ РВп о • Степень анизотропности породы оценивается коэффициентом ани- анизотропии Ч = УРвп±/РвпЦ>Ь (И.14) а величина ее удельного сопротивле- сопротивления характеризуется РсР==УРвп.1_Рвпц- (П.15) В нефтегазонасыщенной породе нефть или газ, частично замещая Рис. 20. Зависимости коэффициента увели- увеличения сопротивления Рн породы от содер- содержания в них воды kB (нефти kH или газа kT) §,% ф SO К к % (по 1 — 3 — песчано-глинистые соответственно гидро- фильные, слабогидрофобные и гидрофобные породы ( В ** Дахнову), 4 — карбонатные породы (по в поровом пространстве воду, повышают ее удельное сопро- сопротивление. В этих условиях рнг зависит от содержания в ее порах нефти, газа и воды, характера их распространения в поровом про- пространстве, минерализации пластовой воды, пористости и структуры порового пространства и т. п. Для полного или частичного исклю- исключения факторов (минерализации пластовых вод, пористости и струк- структуры порового пространства), влияющих на величину рнг, вместо него рассматривают отношение рнг Рвп A1.16) где рнг — удельное сопротивление породы, поры которой запол- заполнены нефтью (газом) и остаточной водой; рвп — удельное сопро- сопротивление той же породы при условии 100 %-ного заполнения ее пор водой. Величина Рп показывает, во сколько раз увеличивается удель- удельное сопротивление породы, насыщенной нефтью или газом при ча- частичном заполнении ее пор водой, и называется коэффициен- коэффициентом увеличения сопротивления. Для неглини- 42
¦тых пород Ри зависит не только от степени их насыщения водой, но if от характера распределения в поровом пространстве воды, нефти и газа. В связи с этим величина Рн часто называется п а - р а м е т р о м насыщения. Между Рп и kB существует обратная степенная зависимость Рн= —, A1.17) К rie n показатель, характеризующий структуру токопроводящих каналов нефтегазонасыщенной породы, зависящий от ее литолого- оегрографпческих особенностей и структуры порового пространства, физлко-хилшческих свойств нефти (газа) и воды и их распределе- распределения в порах. Численное значение показателя п может изменяться от 1,73 до 4,33. Если частицы воды, насыщающей породу, представляют собой связанную систему (что характерно для гидрофильных по- пород}, показатель п имеет меньшие значения, при разобщении ча- спш воды в порах (гидрофобные породы) п возрастает. На прак- практике для гидрофильных пород п часто принимается равным 2. Так как kB = 1 — knr (где kHr — коэффициент нефтегазонасы- шепности), то Рн = . A1.18) \1 кяг) На рис, 20 показаны усредненные кривые зависимостей Ри = ~~ f (kB) для различных нефтегазоносных пород. Анализ много- многочисленных исследований этих зависимостей показывает, что для их описания не всегда достаточно одного уравнения. В таких слу- случаях они аппроксимируются по крайней мере двумя уравнениями соответственно в области высоких и низких значений kB: Рн-—; (И.19) Рн = —. (П.20) Зависимости Ри == / (kB), полученные экспериментально для конкретного типа продуктивных отложений, используются для оценки их нефтегазонасыщенности методом сопротивлений и на- находят широкое применение в практической работе. Удельное сопротивление глинистых пород Для глинистой водонасыщенной породы пропорциональности между ее удельным сопротивлением рвп гл и удельным сопротив- сопротивлением насыщающей воды рв нарушается. Это связано с тем, что электропроводность такой породы определяется не только прово- проводимостью воды, но и поверхностной проводимостью глинистых ча- 43
стиц, точнее, гидратационной пленки, покрывающей их поверх- поверхность. Поверхностная проводимость проявляется тем значитель- значительнее, чем выше глинистость породы и меньше минерализация насы- насыщающей воды. Вследствие этого относительное сопротивление гли~ нистых пород в отличие от неглинистых зависит не только от их пористости и структуры пор, но и от их глинистости и минерализа- минерализации насыщающих вод. Относительное сопротивление глинистой породы, соответствующее насыщению высокоминерализованной во- водой, при которой поверхностная проводимость минимальна, на- называют предельным Рп. Учет влияния поверхностной про- проводимости глин на относительное сопротивление осуществляется при помощи коэффициента поверхностной проводимости П == *п QJ OS 0,5 0,2 О 10 15 20 Рис. 21. Зависимости коэффициента поверхностной проводимости Я от глини- глинистости kVJl и удельного сопротивления поровой воды рв для нижнемеловых продуктивных отложений Прикумского нефтегазоносного района Ставро- где Рк — кажущееся относительное сопротивление пород, насы- насыщенных менее минерализованной водой. Зависимость параметра П от глинистости пород и минерализа- минерализации насыщающих вод изображена на рис. 21. На основании экспе- экспериментальных работ [7 ] установлено, что электропроводность гли- глинистой породы с рассеянным глинистым материалом или где А 44 1 Ргл =л+-^, (Н.22> (П.23> Рвп. гл гъ отрезок, отсекаемый продолжением прямой на оси орди-
ат и характеризующий долю проводимости глин; В — угловой коэффициент; величина В/рв определяет долю проводимости пор шхс. 22). Как видно, зависимость 1/рвп = f (l/pB) для чистой породы шямолинейна во всем интервале рассматриваемой функции, а 17Я глинистой с рассеянным глинистым материалом,— начиная 1ишь с некоторого значения проводимости воды. Согласно формулам (П.4) и A1.23) имеем р Рвп гл 1 A1.24) Для высокоминерализованных растворов, когда рв -> О, = lim 1 В Из уравнений A1.21) и A1.24) рк в A1.25) В тех случаях, когда гли- глинистый материал присутствует в породе в виде прослоев Ае Оп 0,2 10 ' A1.26) Рис 22. Зависимости проводимости водонасыщенной породы 1/рПв от В породе с рассеянным глини- проводимости воды 1/рв, заполняю- стым материалом условия рас- щей поровое пространство ПрОСТраненИЯ ТОКа ПрИМерНО Та- Песчаник- 1 — слабоглинистый, 2 — силь- 1 1 ноглинистый кие же, как в водонасыщенной породе с kn = krjl и рв = ргл. Исходя из этого, согласно формуле (IIЛ 0) получим, что электропроводность такой породы равна k™JaprJl> где ант — коэффициенты, зависящие от свойств породы. В этом случае 1 Рвп.гл ЯРгл 1 Если принять а = 1, т — 2, тогда Рвп гл _гл Ргл РпРв A1.27) A1.28) 45
На основании формул A1.26) и A1.28) уравнения проводимости для глинистых нефтегазонасыщенных пород при слоистом и рассе- рассеянном содержаниях глинистого материала соответственно будут иметь вид Рнг гл Ргл РпРв 1 k2 1—k где kB — коэффициент водонасыщеиности чистого компонента. Согласно формулам A1.26), A1.29) и A1.30) коэффициенты уве- увеличения сопротивлений для этих случаев будут Р Как видно из формул A1.31) A1.32), при одном и том же значе- значении къ (knr) чистого компонента коэффициент увеличения сопро- сопротивления может колебаться в широких пределах в зависимости от изменения глинистости, удельного сопротивления глинистого ком- компонента и характера распределения глинистого материала в по- породе. В глинистых породах по данным работы [2] повышение темпе- температуры приводит к более заметному снижению удельного сопро- сопротивления по сравнению с чистыми породами благодаря повышению ионной концентрации свободных растворов в порах, за счет разру- разрушения слоя адсорбированных ионов с увеличением температуры, и относительно низкому A02—106 Ом-м) удельному сопротивлению большинства глинистых минералов, являющихся полупроводя- полупроводящими. Удельное сопротивление пород с трещинной и каверновой пористостью Породы с трещинной и каверновой пористостью весьма разно- разнообразны по составу и строению. Наряду с межзерновой (первичной) пористостью kn. M значительную роль играют поры вторичного про- происхождения kn вт — трещины, каверны и другие пустоты выщела- выщелачивания. Каверны (изолированные и полуизолированные пустоты) заметного влияния на удельное сопротивление пород не оказы- оказывают. Наличие трещин, заполненных электролитом, вызывает су- существенное снижение сопротивления по сравнению со снижением сопротивления, обусловленным межзерновой пористостью такого же объема. 46
В природе встречаются породы с различными системами трещин / 1С 23) [7]. Пользуясь правилом Кирхгофа для параллельных и последовательных проводников, можно получить выражения, связывающие удельное сопротивление трещиноватой породы ртр - удельным сопротивлением ее нетрещиноватой части рм, с величи- величиной трещинной пористости kn.TP и удельным сопротивлением рж жидкости, заполняющей ее. Эти выражения для различных систем трещин BВ, 2ГВ, ЗВП, X) будут иметь соответственно следующий вид: i п тр . п. тр , /т т оо\ __ 1 ? [И.Од) Ртп 2Рж Рм /// / Рис. 23. Различные системы трещин Системы: 1В — вертикальная, 1Г — горизонтальная; 2В —- две взаимно перпендикуляр- перпендикулярные чертикальные, 2ГВ — вертикальная и горизонтальная; ЗВП — взаимно перпенди- перпендикулярные, X — хаотическая трещиноватость 2k п. тр ^п. тр A1.34) Приведенные уравнения проверены А. М. Нечаем в лаборатор- лабораторных условиях на образцах с искусственными системами трещин; установлено удовлетворительное совпадение экспериментальных и расчетных данных. Как видно, порода с двумя взаимно перпен- перпендикулярными трещинами (см. рис. 23, система 2В, 2ГВ) является анизотропной, с тремя взаимно перпендикулярными трещинами и хаотическим их распределением (см. рис. 23, системы ЗВП, X)—• изотропной. Для хаотического распределения трещин формула A1.34) в уп- упрощенном виде без существенной погрешности может быть запи- записана 1 п. тр ЛР* Pi/ A1.35) 47
Множитель А в зависимости от геометрии систем трещин изме- изменяется от 1 до 2; для изотропной породы он равен 1,5. Формулу A1.35) можно представить в виде К трРм + ЛРж Так как удельное сопротивление водонасыщенной матрицы, со- согласно уравнению A1.10), рм = apjk% M, то выражение A1.36) будет Ртр - — в- • A1.37) Кп Mi ^п.тр а Рж Анализ уравнения A1.37) показывает, что: 1) с уменьшением трещинной пористости kn, тр удельное сопро- сопротивление чисто трещинной породы ртр приближается к величине удельного сопротивления ее матрицы рм; при ku тр = 0 ргр = рш и формула (II.37) преобразуется в уравнение A1.10). 2) с увеличением удельного сопротивления жидкости в трещи- трещинах рж (при насыщении трещин нефтью или пресным раствором) ртр и рм незначительно отличаются друг от друга; 3) влияние трещиноватости на удельное сопротивление пород возрастает с увеличением удельного сопротивления ее матрицы рм. При межзерновой пористости матрицы более 8—10 % влияние тре- трещинной пористости на удельное сопротивление заметно снижается и порода может рассматриваться как чисто поровая. В природных условиях чисто трещинные породы встречаются крайне редко. На- Наряду с трещинами в породе развиваются и другие пустоты вторич- вторичного происхождения: каверны, карсты, каналы растворения и т. п. Эти пустоты могут существовать в породе изолированно или сооб- сообщаясь между собой трещинами, которые сами по себе имеют тен- тенденцию к сужению, расширению, образуя разновидности вторич- вторичных пустот. Наличие указанных пустот в породе оказывает допол- дополнительное влияние на ее удельное сопротивление, которое учиты- учитывается с помощью формулы Максвелла [7 ] Hlia Vl C + 21/) Pl + 2l/p2 где pi, 2 — удельное сопротивление среды, в которой беспорядочно расположены сферические включения с удельным сопротивлением р2; V—доля объема, занятого средой с удельным сопротивлением pi. После введения в формулу A1.38) обозначений рх = ртр, р2 = = Рж, V = 1—kn. кав получим выражение для определения удель- удельного сопротивления породы с хаотически распределенными в ней трещинами и изолированными пустотами (кавернами): Ртр + 2Рж1 + ^п.кав(Рж1-Ртр) n m от Ртр. кав = — — — Ртр , {IL.W) РтР + 2Рж1~2*п.кав(Рж1-Ртр) где рж1 — удельное сопротивление жидкости, насыщающей изоли- 48
рованные пустоты (каверны); ртр — удельное сопротивление по- ооды, обусловленное влиянием только трещин, заполненных жид- жидкостью с удельным сопротивлением рж. Экспериментальные работы А. М. Нечая показали, что резуль- результаты расчетов по формуле A1.39) удовлетворительно согласуются с данными лабораторных измерений, размеры и формы каверн ока- оказывают малое влияние на величину удельного сопротивления, а воздействие отдельных каверн на распределение силовых линий тока незначительно, если расстояние между ними в 2—3 раза и бо- более превышает их размеры. В предельном случае, когда изолированные пустоты заполнены высокошшерализованной водой (рж1 ->- 0), выражение A1.39) при- приобретает вид 1 А п кав /тт .лч Ртр кав - Т~"IT Ртр ¦ A1.40) 1 + 2k Наибольший интерес представляют породы, в которых каверны сообщаются между собой по трещинам. Оценить удельное сопротив- сопротивление такой породы можно путем подстановки в уравнение A1.40) значения ртр, определяемого по формулам A1.35) — A1.37). При этом A1.42) A1.43) A1.44) тр {~тр. Рт) — . кав ^ ' В з кав А. п.тр 1 р. кав рм+. +" - тр . ^тр. кав ^ Где В — A — *п.кав)/A + 2 ^п.кав); ПрИ kn, Кав< C — 5 %) ЗНа- ченне В приравнивается к единице. Формулы A1.41) — A1.43) применимы для пород с водонасы- щенной практически непроницаемой гранулярной матрицей. Фор- Формула A1.44) позволяет определять удельное сопротивление пород со смешанной пористостью (трещинной, каверновой и грануляр- гранулярной). Это выражение может быть использовано для оценки удель- удельного сопротивления смешанного коллектора с нефтенасыщенной матрицей; в качестве гранулярной пористости матрицы (блока) принимается условное значение kn. M, соответствующее ее водона- сыщенному объему. 49
Влияние эффективного давления на удельное сопротивление пород В естественных условиях залегания горные породы находятся под действием эффективного давления, которое определяется no- формулам A.12) и A.15). В результате происходит обратимая де- деформация пород, сопровождающаяся уменьшением их пористости и изменением геометрии порового пространства, т. е. увеличением извилистости поровых каналов преимущественно из-за их суже- сужения. Изменение этих пара- о и >?9 jo /7зф, чип метров с ростом давления обусловливает и соответству- соответствующее обратимое изменение (в сторону увеличения) удель- удельного сопротивления пород, которое зависит от величины давления и их литолого- петрографических особенно- особенностей. Наиболее интенсивное изменение удельного сопро- сопротивления происходит под дей- действием /?Эф до 30 — 40 МПа. При дальнейшем увеличении рЭф интенсивность роста со- сопротивления снижается и стремится к некоторому пре- пределу до наступления необра- необратимых разрушений в ске- скелете породы (рис. 24), Рис. 24. Зависимости относительных изменений пористости (/) и электропро- электропроводности Аа/ог B) песчано-глинистых пород от эффективного давления (ниж- (нижнемеловые отложения Прикумского нефтегазоносного района Ставрополья) Экспериментально установлено, что максимальные изменения сопротивления при /7эф — const характерны для низкопористых глинистых разностей известняков, песчаников и алевролитов. Ва* риадии пористости, структуры поровых каналов и удельного со- сопротивления под действием давления приводят и к соответствую- щему изменению характера зависимости Рп = / (kn). Если под действием давления пористость породы уменьшается на Д&п, а т и а соответственно увеличивается и уменьшается на А/n и Да, то уравнение A1.10) принимает вид а — Да A1.45) Выражение A1.45) разделим на A1.10), тогда э(^эф) — (Аа/а) A1.46) 50
Исходя из допустимой погрешности, л АМДот I к) * Тогда ,_ ч При одном и том же давлении величина Ри изменяется значи- значительно больше, чем пористость, что объясняется преимуществен- преимущественным влиянием увеличения извилистости токопроводящих каналов. Зависимости Рп = f (ku) при различных давлениях показаны на рис. 25. Определение удельного сопротивления пород Сопротивление пород в лабораторных условиях может изме- измеряться различными способами [10]. Наиболее широкое примене- применение из них получили мостовые двухэлектродные способы регистра- регистрации на переменном токе, при которых осуществляется прямое оп- определение сопротивления образца породы в равновесном режиме. Одна из мостовых схем, позволяющая измерять сопротивление образцов пород с компенсацией емкостной составляющей, показана на рис. 26. Схема четырехплечного моста включает три активных сопротивления Rc, Rly R2 (магазины сопротивлений Р517 типа МСР) и Rx исследуемого образца породы. В плече сравнения па- параллельно с активным сопротивлением Rc включается и реактив- реактивное (емкость С — магазин переменной емкости Р544 типа МЕРП). Питание мостовой схемы осуществляется от звукового генератора ЗГ-10 (или ЗГ-12, ЗГ-34). Измерение величины Rx производится уравновешиванием моста по напряжению и фазе. В качестве ин- индикатора равновесия используется осциллограф или электронной индикатор нуля переменного тока Ф-510, включаемый в диагональ моста (см. рис. 26). Рассмотрим последовательность измерений. 1. Образцы пород правильной геометрической (цилиндрической или кубической) формы, предварительно экстрагированные от нефти и отмытые от остаточных пластовых солей, насыщают под вакуумом раствором заданной минерализации и химического со- состава. Образцы неустойчивых пород перед насыщением покрывают эпоксидной смолой или клеем БФ. 2, Непосредственно перед измерением образец породы извле- извлекают из раствора, обтирают с боковой поверхности фильтроваль- фильтровальной бумагой и устанавливают между электродами кернодержа- теля.При этом необходимо, чтобы с помощью пружин или винто- винтового зажима электроды прижимались к торцам образца с постоян- постоянным усилием. Для улучшения контакта между торцами образца и электродами устанавливают прокладки из двух-трех слоев филь- 51
тровальной бумаги, насыщенной тем же раствором, что и образец^ При необходимости сопротивление прокладок может быть изме- измерено и учтено. Электроды кернодержателя должны обладать минимальным переходным сопротивлением, которым можно пренебречь по срав- нению с сопротивлением образца. Согласно этому условию приме- применяют электроды из нержавеющей стали, периодически обрабаты- обрабатываемые в течение 5—10 мин переменным током частотой 50 Гц и на- напряжением до 60 В в растворе NaCl, которым насыщался образец. После работы электроды хранят в том же растворе. wo 5 10 го /гп,% Рис. 25. Зависимости относительного сопротивления Ри от пористости ku при различных эффективных давлениях рЭф для нижнемеловых песчано* глинистых отложений Прикумского нефтегазоносного района Ставрополья 1 - Рэф - 0, Рп = 1/&1/; 2 - рэф= 10 МПа; Ри = 0,8/а?; 3 - рэф= 20 МПа, Ри ~ = 0,65/4'15 ; 4-Яэф =30 Ь2.3 Рис. 26. Принципиальна!* схема моста переменного тока для измерения со- сопротивления образцов пород с компенсацией емкостной составляющей 3. Измеряют сопротивление образца, для чего: а) устанавли- устанавливают переключатели магазинов сопротивлений и емкости в нулевое положение; б) от звукового генератора на измерительный мост по- подают ток частотой 1000 Гц, плотность тока через образец не должна превышать 2 мА/см2, что обеспечивает минимальное переходное сопротивление электродов; в) при включенном осциллографе сна- сначала на грубом пределе, а затем на пределах более высокой чувст- чувствительности добиваются изображения минимальной амплитуды си* 52
нусонды переменного тока с помощью Rc. Изменяя поочередно Сопротивления RU R2 и емкость С, достигают полной компенса- компенсация 0 которой судят по преобразованию на экране осциллографа нуС0Пды в прямую линию на самом чувствительном пределе. Rx определяется из соотношения: п Rx=Rc-^-. (П.48> Удельное сопротивление образца породы рассчитывают по фор- муЛе (II. 1). Значение емкости С в плече сравнения, соответствую- соответствующее моменту компенсации, при необходимости можно использо- использовать для расчета диэлектрической постоянной образца (см. § 9). § 6. ЭЛЕКТРОХИМИЧЕСКАЯ АКТИВНОСТЬ Электрохимические процессы, протекающие в горных породах, вызывают их поляризацию. К ним относятся диффузионно-адсорб- диффузионно-адсорбционные, фильтрационные, окислительно-восстановительные про- процессы и процессы, связанные с действием внешнего электрического поля. В зависимости от фактора, вызывающего поляризацию, раз- различают диффузионно-адсорбционную, филь- трационную, окислительно-восстанови- трационную, окислительно-восстановительную и вызванную электрохимические активности пород. Окислительно-восстановительная ак- активность по сравнению с другими перечисленными активностями для осадочных пород имеет подчиненное значение и здесь не рас- сматривается. Диффузионно-адсорбционная активность Водные растворы солей представляют собой электролиты, со- содержащие положительно и отрицательно заряженные ионы. На контакте двух электролитов разной концентрации, вследствие раз- различной электролитической подвижности их ионов, возникают диффузионные потенциалы (рис. 27). Электродвижущие силы (ЭДС) диффузии этих потенциалов в случае контакта простых и не отличающихся по химическому составу электролитов оцениваются уравнением Нернста рт пи — п/о С, Ел=-Щ !^_^_1п-1, (II 49) F nKzKti + nazav C2 где R — универсальная газовая постоянная [8,314 Дж/(моль- Кл) ]; F ~~~ число Фарадея (96 500 Кл/моль); Т — абсолютная темпера- температура (в К); пк ияа — числа катионов и анионов, на которые диссо- диссоциирует одна молекула электролита; zK и га — валентности ка- катиона и аниона; и и v — электролитические подвижности катиона и аниона (в Ом-см2/моль); Сх и С2 — концентрации контактирую- контактирующих электролитов (в моль/л); Сг>С^ тогда ?д определяет потен- потенциал второго электролита по отношению к потенциалу первого. 53-
Из формулы A1.49) следует, что диффузионные ЭДС возникают лишь при условии, если пки Ф nav и Сг ф С2. Для одно-однова- одно-одновалентного электролита (например, раствор NaCl ->- Na+ + Cl~)* когда лк = да = 1 и zK = za = 1, формула A1.49) принимает вид RT u—v - v (П.50) С учетом числовых значений констант R и F, замены натураль- натурального логарифма десятичным при Т = 291 К и ЭДС в м'В получим NciGl п е ф ? NaCL Рис 27. Схема возникновения диф- диффузионных потенциалов на контак- контакте двух электролитов С± и С2 — концентрации двух растворов хлористого натрия (Сг > С2), П — изме- измерительный прибор; 5i и Эо, — электроды или .^58(Мк-Мг)Ы^Л1Ш) где Л^к = и/(и + t/), Na ^ = v/(u -\- v) — числа переноса катионов и анионов.** Выражение перед логариф- логарифмом в формулах A1.50) и A1.51) при заданной температуре не изменяется вследствие постоян- постоянства подвижностей ионов. После обозначения этого выражения через величину /Сд формула A1.52) принимает вид A1,53) Величина Кд, называется коэффициентом диффу- диффузионной ЭДС и зависит от температуры и химического состава контактирующих электролитов. Для растворов NaCl при Т = 18 °С подвижность катиона Na+ и = 43,3, аниона Cl~~ v = 65,2 (табл. 2) и /Сд = — 11,6 мВ. Формула A1.53) справедлива для растворов низкой концентра- концентрации, в которых практически отсутствует взаимодействие между ионами. Для реальных растворов в формуле A1.53) необходимо использовать вместо концентрации Сг и С2 значения активностей * В растворе NaCl ион С1~ примерно в 1,5 раза подвижнее иона вследствие чего раствор с концентрацией С2 < Сг заряжается отрицательно (см рис. 27). * * Число переноса характеризует относительную подвижность или скорость движения аниона либо катиона (при сумме их скоростей, равной единице) и ту долю общего количества электричества, прошедшего через электролит, которая будет перенесена соответственно анионом или катионом* 54
Таблица 2 Л 1ШН0СТЬ ионов в воде при бесконечном разбавлении (Омсм2/моль) Ион К* €!~~ 1.2SOf~ 0 240,0 40,4 27,0 41,1 41,0 105,0 Температура Т 18 315,0 64,5 43,3 65,2 68,5 174,0 25 350,0 73,5 50,11 76,3 80,0 198,0 50 465,0 115,0 "82,5 116,0 125,0 284,0 100 644 206 154 207 234 439 Температурный коэффициент подвижности ионов, °о на 1 СС 1,54 2,17 2,44 2,16 2,27 1,80 d ~, f-fii и &2 — !<2ръ гДе /i и /2 — коэффициенты активности, учитывающие влияние сил взаимодействия между ионами в реаль- реальных растворах. При этом ?Д = /Сд^ —. A1.54} В значительных пределах изменения концентраций отношение активностей а1/а2 растворов практически равно обратному отноше- отношению их удельных электрических сопротивлений p2/pi, следова- следовательно Е д =/сд lg J?*-. (П.55) Pi Значения диффузионного потенциала для растворов NaCl раз- различной концентрации, вычисленные по формуле (П.55) при 18 °С приведены в табл. 3. Величина диффузионного потенциала для разных по химическому составу растворов различна. В част- частности, для растворов КС1, вследствие почти одинаковой подвиж- подвижности ионов хлора и калия, диффузионный потенциал близок к нулю. Если два раствора различной концентрации разделить тонко- дисперсной породой (например, глиной), то твердые (особенно глинистые) частицы последней из-за ионно-сорбционной способно- способности адсорбируют ионы определенного знака и изменяют их концен- концентрацию в свободных растворах. Под ионно-сорбционной способностью понимается свойстве влажных пород избирательно адсорбировать из насыщающих элек- электролитов наряду с молекулами воды ионы определенного знака и, как следствие этого, образовывать двойные электрические слои на границах раздела фаз. Двойной электрический слой состоит из потенциалоп- ределяющего слоя ионов и слоя противо- ненов. Потенциалопределяющий (адсорбционный) слой ионов 55'
Таблица Величина Яд для электролитов NaCl различной концентрации С3, нормаль 1 0,1 0,01 0,001 С1? нормаль 1 — 10,7 —21,8 —33,1 0,1 -11,1 —22,5 0,01 — 11,3 Коэффициент активности / 0,65 0,80 0,92 1,00 Примечания. 1 Нормальная концентрация (нормальность) раствора выража- выражается количеством вещества, растворенного в 1 л раствора, в молях. 2. Знак минус означает, что электрод, находящийся в менее концентрированном растворе, имеет отрицательный знак. прилегает наряду с молекулами прочно связанной воды непосредст- непосредственно к поверхности минеральных частиц твердой фазы породы. Далее располагается прочно связанный с ним неподвижный слой противоионов и затем более рыхлый, менее прочно связанный и более подвижный диффузный слой. Установлено, что между ионами адсорбционного и диффузного слоев происходит непрерывный равновесный обмен. Это приводит к изменению диффузионных ЭДС как по величине так и знаку. Электродвижущая сила, возникающая в цепи электро- электролит с концентрацией С1 — тонкодисперсная порода — тот же элек- электролит с концентрацией С2, носит название диффузионно- адсорбционной (рис. 28). Ее величина при наличии двух растворов одинакового химического состава, но различной кон- концентрации изменяется от значения диффузионной ЭДС до предель- предельного, определяемого уравнением рт р ?Да = —In-^, (П.56) zF C2 где г — валентность электролита. Уравнение A1.56) получено из уравнения A1.49) при условии, что анионы неподвижны (v — 0) и валентности катионов и анионов равны (гк = га). Из выражения A1.52) при NK = 1, Na = 0 по- получаем предельное значение Г±. A1.57) По аналогии с диффузионной диффузионно-адсорбционная ЭДС рассчитывается с помощью уравнения Яда = /Сда lg -^- = Кда lg -^ , а2 pi где /Сда — коэффициент диффузионно-адсорбционной ЭДС. 56 A1.58)
Величина /Сда в отличие от Кд зависит не только от температуры: химического состава растворов, но и от свойств породы, через которую происходит контакт растворов. Численно /Сда ~ Кд + Лда, и, следовательно, Pi где АДа — диффузионно-адсорбционная активность породы. Из уравнения A1.59) вытекает, что диффузионно-адсорбцион- диффузионно-адсорбционость лая активность А да "¦ Pl -/Сд == Кда — A1.60). Pi NaCL • •• • С I 1 I Рис. 28. Схема возникновения диф- диффузионно-адсорбционных (мембран- (мембранных) ЭДС, когда два раствора раз- различной концентрации разделены тонкопористой (глинистой) перего- перегородкой (Сг > С2) Рис. 29. Зависимость диффузионно- адсорбционной активности Лда от приведенной емкости обмена qn для песчано-глинистых пород [1] Таким образом, диффузионно-адсорбционные потенциалы, воз- возникающие на границе пород с электролитами, могут существенно отличаться от диффузионных. Свойство же горных пород создавать диффузионно-адсорбционные потенциалы, вследствие их способ- способности адсорбировать из растворов наряду с молекулами воды ионы определенного знака, называется диффузионно-адсорб- диффузионно-адсорбционной акт и-в н о с т ь ю. По результатам многочисленных лабораторных исследований установлена достаточно тесная связь между диффузионно-адсорб- диффузионно-адсорбционной активностью Лда и ионно-сорбционной способностью по- род> характеризуемой приведенной емкостью обмена qu (рис. 29). Приведенная емкость обмена (поглощения) характеризует число молей поглощенных катионов, приходящихся на единицу объема 57
(м3, см3) порового пространства породы ^п="г"=—'~^~i бм> (H.etj; ^п лп где 6М — плотность твердой фазы породы; QM. e — массовая ем* кость обмена, характеризующая число молей поглощенных катио- катионов на единицу массы (кг, г) породы и зависящая от химического состава жидкой и твердой фаз, а также степени дисперсности по- последней; Qn — объемная емкость поглощения. v ' 30 20 10 .<i >; * О V' 30 - 2U - Ю ~ 3 о у р ./ .7 1 у /• ^ 1 I I 10 г{7 JG (9 U,Z Oh 0,6 0,8 уп Рис. 30. Зависимость диффузионно-адсорбционной активности Лда пород от их объемной ?гл (а) и относительной %л (б) глинистости [1] / — породы-коллекторы, 2 ~ породы-неколлекторы, 3 — линия регрессии Степень дисперсности пород в значительной мере зависит от иж глинистости. Исследования образцов различных терригенных по- пород показали, что отличие диффузионно-адсорбционного потен- потенциала от диффузионного тем больше, чем выше содержание в нш глинистого материала. Глинистость пород, способствующая фор* мированию субкапиллярных пор, является основным фактором, обусловливающим различия величин ?да и ?д. Зависимости Ат от krJl и т]гл, полученные экспериментально на образцах, представ- представлены на рис. 30 [1 ]. Величина Лда пород зависит не только от содержания в ни глинистого, но и другого по вещественному составу цементирую- цементирующего материала. Наиболее высокой активностью характеризуется железисто-сидеритовый и глинисто-серицитовый цементы, мевь* шей — опалово-халцедоновый и карбонатный. Влияние глинистости на удельную поверхность 5УД, эффектив- эффективный диаметр 4$, коэффициенты пористости kn, проницаемости &пр, поверхностной проводимости П и другие, с одной стороны, и на величину Лда — с другой, обусловливает наличие зависимо- зависимостей между диффузионно-адсорбционной активностью, и перечне* 58
параметрами пород. На рис. 31 показана связь между %1 1Ф Численные значения диффузионно-адсорбционной активности. 21Я различных горных пород изменяются от — 25 до 70 мВ. Наи- Наименьшими величинами Лда обладают чистые разности песков, пес- песчаников, алевролитов, пористых известняков и доломитов. Фильтрационная активность iipn фильтрации жидкости через пористую среду (породу) воз- возникают потенциалы течения — фильтрационные по- х е н и и а л ы. Их образование обусловлено наличием электро- "дш" 3 О 20 10 X 10 50 100 Рис* 31. Связь между диффузионно- адсорбционной активностью Лда и коэффициентом проницаемости ?пр для песчано-алевритовых пород Ильско-Холмского района Красно- Краснодарского края [2] Рпс. 32. Падение потенциала в двой- двойном электрическом слое при переме- перемещении жидкой фазы относительно твердой ц — потенциал частиц твердой^фазы, гх — юлщпна адсорбционного слоя (слой Гель- мгольда), / 2 — толщина диффузионного слон (слои Гуи); стрелкой показано на- иравление фильтрации кнштического потенциала ? на границе между подвижной и не- неподвижной частями двойного слоя (рис. 32, линия АВ) и свобод- свободного раствора. Механизм образования фильтрационных потенциа- потенциалов схематично можно представить следующим образом. Под воз- воздействием перепада давлений Ар в поровом канале происходит дви- движение свободного раствора, заполняющего среднюю часть поро- вого канала. Одновременно с ним в направлении течения переме- перемещается и подвижная часть внешней диффузной обкладки двойного электрического слоя. Это приводит к смещению электрокинетиче-
•ского потенциала ? в направлении течения, вследствие чего на концах капиллярного порового канала возникает разность потен- потенциалов. Если толщина двойного электрического слоя мала по сравнению •с радиусом капиллярного порового канала, величина фильтра- фильтрационного потенциала определяется формулой Гельмгольца где еж—диэлектрическая постоянная жидкости; рж — удельное сопротивление фильтрующейся жидкости; \х — вязкость; Ар —. избыточное давление, под действием которого происходит филь- фильтрация жидкости. Согласно формуле A1.62), величина фильтрационного потен- потенциала пропорциональна перепаду давлений Ар и удельному со- сопротивлению жидкости рж. С увеличением минерализации жидко- жидкости значение ?ф уменьшается вследствие снижения величин рж и ?. Параметры 8, рж и \i, входящие в формулу, являются положи- положительными при всех условиях и не изменяют знака Еф. Последний зависит от знака дзета-потенциала ? и направления движения жидкости, т. е. направления избыточного давления Ар. В общем случае величина фильтрационного потенциала для породы где Лф.?п — фильтрационная электрохимическая 'активность по- породы, зависящая от структуры порового пространства и свойств •фильтрующейся жидкости. Установлено, что потенциалы фильтрации возникают в различ- различных по своей литологии породах, в том числе и глинах. Однако наибольший интерес они представляют в породах-коллекторах. Определение диффузионно-адсорбционной активности пород Диффузионно-адсорбционную активность Лда пород опреде- определяют по результатам измерений диффузионных ЕА и диффузионно- адсорбционных ?да потенциалов с помощью схем, изображенных на рис. 33. Схема для регистрации ?да (см. рис. 33, а) включает электро- электрохимическую цепь, состоящую из растворов NaCl с концентрациями Сх и С2 (удельные сопротивления рх и р2), которые разделены об- образцом породы 2, насыщенным раствором концентрации Съ элек>. тролитическую ячейку 1 (парафиновую или плексигласовую ван- ванночку), стаканчики 3 с насыщенным раствором КС1, сифоны 5, а также приборы 7 и в, служащие для измерения разности потенциа- потенциалов. Контакт измерительных приборов с электрохимической цепью осуществляется с помощью каломелевых (или хлорсеребряных типа ЭВЛ-1М1, хлорталлиевых типа ЭВЛ-5М1 и 5М2) электродов 4Р 60
аемых в СТаканчики 3. Эти электроды обладают небольшим •¦тойчивым во времени потенциалом. Чтобы исключить перете- П * растворов между ячейкой и стаканчиками, сифоны 5 запол- Кяют смесью насыщенного раствора КС1 с агар-агаром. пчя поддержания стабильности концентраций растворов на оцах образца (особенно с хорошей проницаемостью) целесооб- чзво создавать проточные границы. Для этого над ячейкой слева ^* справа от образца устанавливаются капельницы (на рисунке не ^оказаны), наполненные соответствующими растворами, а в стен- стенах отделений ячейки делаются отводы 6. рис, 33. Схемы для измерений диффузионно-адсорбционных (а) и диффузионных (б) ЭДС а. / __ электролитическая ячейка {парафиновая ванночка), 2 — об- образец породы, 3 - стаканчики с насыщенным раствором КС1, 4 — кахмелевые электроды, 5 — агар- «rjpOBbie сифоны, 6 — отверстия a,Wстока раствора; б. 1 — элек- тоолишческая ячейка 2 — солевой ЯОС1ИК (сифон), 3 — стаканчики <г насыщенным раствором КС1, 4—- шшомелевые электроды, 5 — агар- «гвро»ые сифоны, 6 — резиновая ?шгшгка, 7 — потенциометр 5 — гальванометр, 9 — переключатель Отличительной особенностью схемы измерения ?д (см. рис. 33, б) является то, что растворы с концентрациями Сх и С2 (удельные со- сопротивления рх и р2) в ячейке 1 полностью разобщены, а непосред- непосредственный контакт между ними возникает с помощью мостика (си- (сифона) 2, заполненного менее концентрированным из контактирую- контактирующих растворов. Определение Лда производят в следующем порядке. L Образец породы, экстрагированный от нефти и остаточных солей, насыщают под вакуумом пластовой водой или раствором NaCl с концентрацией Сх (удельное сопротивление рг) и заделы- заделывают в парафиновой ячейке так, чтобы сообщение между ее отделе- отделениями происходило только через образец. 2. Наливают в левое отделение ячейки раствор NaCl с концен- концентрацией Сх и, убедившись, что образец заделан герметично, в пра- правое отделение заливают раствор NaCl с концентрацией С2 = ~~~ о j/ 1U. 3. С помощью капельниц в отделениях ячейки устанавливают одинаковые уровни медленно протекающих растворов. 61
4. Собирают схему полностью, как это показано на рис. 33, а и производят измерение величины Еда1 (в мВ) до получения ее стабильных во времени значений. Затем электроды в стаканчиках меняют местами и, повторив измерения, получают ?да 2 (в мВ), Значение Ет (в мВ) получают как алгебраическую полуразность Таким образом исключают влияние электродной разности потен- потенциалов на результаты измерений. 5. Заменяют электролитическую ячейку для измерения Ет на ячейку для измерения Ед. Наливают в левое ее отделение раствор NaCl с концентрацией Сг (удельное сопротивление рг) , а в правое— с концентрацией С2 = Сх/10 (р2 = 10 рх) и после создания при помощи перекидного мостика (сифона) их непосредственного кон- контакта измеряют величину ?д (в мВ). 6. По формулам A1.55), A1.58) и A1.60) рассчитывают соответст- соответственно величины /Сд, /Сда и Л да (в мВ). Вызванная электрохимическая активность Горные породы обладают свойством поляризоваться под воз- воздействием внешнего электрического поля. Это свойство связано с электрохимическими явлениями и называется вызванной электрохимической активностью Е (H.64) где Е — напряженность внешнего (поляризующего) поля; Еви — напряженность поля вызванной поляризации. Так как то Здесь AU и At/вп — падения соответствующих напряжений на об- образце; / — длина образца. Величина вызванной электрохимической активности Ав пород зависит от ряда факторов, основными из которых являются химико- минералогический состав, удельная поверхность составляющих зе- зерен, минерализация и химический состав насыщающих поровызс вод. В породах с электронной проводимостью величина Лв возрас- возрастает с увеличением объемного содержания в них электропроводя- электропроводящих минералов. Для осадочных пород, обладающих преимуществ 62
венной ионной проводимостью, при постоянстве минерализации химического состава поровых вод, а также условий измерения /с!ИЫ поляризующего тока, времени воздействия поляризующего тока и времени измерения после его отключения и др.) величина Ав основном зависит от удельной поверхности частиц породы тпс. 34, а). Последнее обусловливает наличие связи между Ав и ^оэффйциентом проницаемости пород (рис. 34, б), что представляет значительный интерес. Определение вызванной электрохимической активности пород Вызванную электрохимическую активность Ав пород опреде- 1яют по результатам измерения приложенной A?/np и вызванной Д(/ва разности потенциалов: Ав = A1.66) / Л8»7о / 0,5 / 0,05 I 0,01 1 5 10 10 50 fOO 600 Рис, 34. Зависимость вызванной электрохимической активности Лв от удельной поверхности зерен 5уд (и)" п коэффициента проницаемости knp (б) пород [2] Pvc. 35. Схема измерения вызванных потенциалов } — образец породы; 2 — электролитическая ячейка (парафиновая ванночка); 3 — элек- электродные стаканчики с насыщенным раствором KCI, 4 — агар-агаровые сифоны; 5 — ба- хареи, 6 — пульсатор, 7 — эталонное сопротивление; 8 —реостат; 9 — потенциометр типа ЭП-1 или Р-307, 10 — гальванометр типа Ml95/3; 11 — переключатель; А и В, М и Л — соответственно токовые (свинцовые) и измерительные (каломелерые) электроды Для измерения величин А(/Пр и А(/вП применяют несколько схем, одна из которых (схема с пульсатором) показана на рис. 35. Сущность методики заключается в том, что через породу при по- помощи пульсатора пропускают короткие импульсы тока, а в проме- ж\ гках между ними измеряют потенциалы вызванной поляризации. Определение Ав производят в следующем порядке. 1. Образец породы 1 цилиндрической (кубической) формы, эк- экстрагированный от нефти и остаточных солей, насыщают пласто- пластовой водой или раствором NaCl эквивалентной концентрации и гер- мешчно заделывают в электролитической ячейке 2 (парафиновой ванночке), заполненной раствором, насыщающим образец. 63
2. Включают пульсатор 6 и при отключенной батарее 5 изме- измеряют разность потенциалов в цепи, соответствующую собственной и электродной поляризациям. 3. От батареи через токовый коллектор пульсатора запитываю? токовые (свинцовые) электроды А и В и при помощи эталонного сопротивления 7 определяют силу тока 4. Находят разность потенциалов, приложенную к образцу Л[/Пр = /#обр> (Н.68). где i?o6p — замеренное сопротивление образца породы. 5. Подключают к измерительным приборам измерительные (ка- ломелевые) электроды М я N и регистрируют вызванную разность потенциалов AUBU. Измерения А(/пр и Д?/вп обычно выполняют при нескольких значениях силы питающего тока, которые регу- регулируют с помощью реостата 8. 6. Рассчитывают Ав по формуле A1.66). § 7. РАДИОАКТИВНОСТЬ Под радиоактивностью понимают самопроизвольное превраще- превращение неустойчивых изотопов химических элементов в другие более устойчивые, которое сопровождается выделением энергии с испу- испусканием а-, |3- и у-лучей. Различают естественную и ис- искусственно вызванную радиоактивности горных по- пород, широко используемые для изучения геологических разрезов скважин методами радиоактивного каротажа. Самопроизвольное превращение одного изотопа в другой называется радиоактивным распадом. Естественная радиоактивность В естественных радиоактивных превращениях основными ви- видами распада являются: а- и р-распады, захват ядром электрона одной из оболочек атома, самопроизвольное деление некоторых тяжелых ядер и др. При распадах радиоактивных ядер и их пере- переходах из более возбужденного энергетического состояния в менее возбужденное или основное возникает ^-излучение. Все виды радиоактивных (а-, Р-, у-) излучений, попадая в ма- материальную среду, испытывают в той или иной мере поглощение* Наибольшему ослаблению подвержены а-лучи, обладающие боль- большой ионизирующей способностью. Поток а-лучей почти полностью поглощается даже листом бумаги или слоем пород толщиной в не- несколько микронов. Поток Р-лучей обладает большей проникающей способностью и полностью поглощается слоем алюминия толщи- толщиной до 8 мм или слоем породы в несколько миллиметров. Гамма- излучение представляет собой высокочастотное коротковолновое электромагнитное излучение, граничащее с жестким рентгенов^ 64
- -им ривае нтс возникает в результате ядерных процессов, и рассмат- рассматтся как поток дискретных частиц у-квантов; у-лучи не откло- отклоива няются ни в электрическом, ни в магнитном полях. % Благодаря своей высокой проникающей способности у-излуче- *ше имеет практическое значение при исследовании разрезов сква- скважин (?~ЛУЧ и полностью поглощаются лишь слоем пород толщиной коло*1 м);; большой интерес представляет изучение характера взаимодействия у-частиц с веществом. Основными процессами взаи- взаимодействия являются: 1) эффект а комигоповского рассеяния; 2) |Ьотоэлектрнческое поглощение; |) эффект образования элек- электрон-позитронных пар. Компто- новский эффект заключается в рассеянии 7~кванта на элек- тронах атома. При соударении с электроном у-квант передает ему часть своей энергии и от- отклоняется от первоначального направления на угол 0 (рис. 36, д). В свою очередь элек- электрон выбрасывается из атома под углом ф к направлению па- падающего кванта. Коэффициент ослабления ^излучения за счет комптонов- ского рассеяния A1.69) rjt$ fiu — плотность поглощаю- поглощающего вещества; Z—атомный (по- (порядковый) номер; А — атомная масса; NQ—число Авогадро, рав- равное 6,023-1023 молекул/(г моль); \it рассчитанный на один электрон. Число электронов в единице объема вещества Рис. 36. Схемы процессов взаимодей- взаимодействия ^-квантов с веществом а — фотоэффект, б — комптоновское рас- рассеяние, в — эффект образования пар коэффициент ослабления, A1.70) Для элементов, составляющих горные породы, отношение ZIA B<;30) является достаточно постоянным и практически равно 0,5. Соогветственно число электронов в единице объема, как и коэффи- коэффициент ослабления, пропорционально плотности вещества. При фотоэффекте у-квант взаимодействует с атомом в целом (см, рлс. 36, б) и передает всю свою энергию одному из электро- электронов его оболочки. Возникающий при этом фотоэлектрон уносит меть энергии у-кванта. Наиболее интенсивное проявление фото- 3 Заказ № 789 65
эффекта наблюдается в веществах, содержащих тяжелые элементу~- и при облучении их ^-квантами малой энергии (близкой к энергий связи электрона). !}\ Электрон-позитронные пары образуются при взаимодействий ^-кванта с ядром атома (см. рис. 36,в). Образовавшийся позитрон через короткий промежуток времени соединяется со свободным электроном, испуская два кванта аннигиляционного излучения энергия каждого из которых равна 0,82-10~13 Дж. Эффект образен вания пар возможен только для жестких у-лучей с энергией, боль-' шей 1,63-10~13 Дж. Под воздействием фотоэффекта, эффекта Комптона и образова- образования пар энергия ^-излучения при прохождении через слой вещества ослабляется. В горных породах поглощение у-излучения природ, ных радиоактивных элементов происходит в основном за счет кошь тоновского рассеяния [см. формулу A1.69)]. Для тонкого моно- монохроматического параллельного пучка ^-лучей его поглощение про- исходит по экспоненциальному закону. После прохождения Y-лучей через слой вещества интенсивность излучения снижается до величины /7, которая может быть рассчи- рассчитана по формуле /7 = /7оехр(-^7/), (ILTOa) где /Vo — первоначальная интенсивность у-излучения; I — тол* щина слоя; \лу — коэффициент ослабления, показывающий, какая часть у-лучей поглощается слоем вещества толщиной, равной еди- единице длины. Коэффициент ослабления пропорционален плотности и зависит от энергии испускаемых у-квантов. На практике процесс поглоще- поглощения оценивается толщиной слоя вещества, в котором интенсив- интенсивность потока у-квантов уменьшается в 2 раза (табл. 4). Таблица 4 Толщина слоя некоторых веществ, поглощающих половину энергии у-квантов (в см) Энергия 7-квантов, МэВ 0,2 1,0 5,0 Вода 5 10 23 Плотный известняк 2,1 4,6 10,0 Железо 0,65 1,58 2,80 Свинец 0,14 0,86 1,47 Согласно экспериментальным данным, число ядер dn радиоак- радиоактивного элемента, распавшихся за бесконечно малый промежуток времени dty пропорционально числу ядер п, не распавшихся к мо- моменту времени t, следовательно где Я — коэффициент пропорциональности, характеризующий ве- 66
ятяость расПада в единицу времени и называемый постоян- постоянной распада. После интегрирования уравнения A1.71) с уче- учетом того, что при t = 0 п = п0, т. е. числу атомов радиоактивного вещества в начальный момент времени, получим п= лоехр( — %t). (II.71 а) Постоянная распада Я, период полураспада Tip и средняя про- продолжительность жизни радиоактивного атома т являются характер- характерными взаимосвязанными величинами и имеют определенные чис- значения для каждого радиоактивного элемента: 0,693,. A1.72) Вероятность радиоактивного распада является свойством са- самого ядра и практически не зависит от физико-химических условий шедшей среды. Распад протекает по экспоненциальному закону Ш-71)) который прослеживается только для очень большого числа распадающихся атомов. Следовательно, интенсивность радиоактив- радиоактивного излучения определяется как статистическая закономерность, обусловленная наличием многочисленных однородных явлений. При неизменных условиях она претерпевает непрерывное измене- пне, колеблясь около некоторой средней величины. Это явление носит название статистической флуктуации (или просто флуктуации). Для", количественной оценки радиоактивности пользуются со- соответствующими единицами измерений. В практике петрофизиче- ских исследований горных пород наиболее часто применяется ве- весовая или объемная единица эквивалентной концентрации радия по у-излучению — г-эквивалент радия на 1 г породы (г-экв Ra/r). Такая единица соответствует концентрации радиоактивных эле- элементов в горной породе, при которой возникает у-излучение такой же интенсивности, как при распаде 1 г Ra. На практике пользуются единицей в 1012 раз меньшей грамм-эквивалента — пико-грамм- зквивалентом радия на 1 г породы^(пг-экв Ra/r = г-экв Ra/г-10~~12). За рубежом радиоактивность горных пород принято выражать в единицах API A пг-экв Ra/r = 16,5 API). Радиоактивность пород в основном связана с со- содержанием в них таких радиоактивных элементов, как уран, то- W 40 ТУ рии, актино-уран, продуктов их распада и изотопа калия шК (в природной смеси изотопов калия его содержится 0,012 %). Кроме этих элементов, источниками радиоактивных излучений, составляющих меньшую долю, чем вышеназванные, являются изо- изотопы рубидия, циркония, индия, лантана, самария, лютеция, ре- шия, висмута и др. Па рис, 37 приведены спектры у-лучей радиоактивных минера- минералов: для калия iqK они имеют одну и ту же энергию, равную 1,46 МзВ; а для рядов тория и урана — радия—множество ее значений.
В литосфере известно более 200 минералов, в состав которых входят уран, торий, радий и калий. Среди изверженных кисльщ породы (граниты, кварцевые диориты, кварцевые порфириты, л&/ париты и др.) обладают наибольшей радиоактивностью, а ультра- ультраосновные (группа перидотитов и пироксенитов) — наименьшей. Радиоактивность осадочных горных пород зависит от радиоак- радиоактивностей их твердой, жидкой и газообразной фаз, а также от со* держания этих фаз и в общем случае Q fa A1.7$ Калии Ряд торил где ць — удельная радиоактивность i-го компонента твердой фазы; ki — объемное содержание /-го компонента в твердой фазе; qjK — удельная радиоактивность флюи- флюидов, насыщающих поровое про- пространство; kn — коэффициент по- пористости. В большинстве случаев плае- товые воды и нефти характери- характеризуются незначительной радиоак- радиоактивностью, а у природных углево- углеводородных газов она практически равна нулю. Исключение состав- составляют лишь пластовые воды, обога- обогащенные солями радия и калия. Радиоактивность твердой фазы пород обусловлена наличием в ее составе собственно урановых^ ториевых и калиевых минералов (уранинит, торит, бреггерит, калий- калийные соли) и минералов (циркон, сфен, ксенотим, монацит, полевые шпаты, слюды, глауконит, гли- глинистые минералы и др.)» содержащих адсорбированные радиоак- радиоактивные элементы. Обычно глинистая составляющая твердой фазы, особенно в кварцевых песчаниках и карбонатных породах, обладает значительно большей удельной радиоактивностью qrjli чем ее соб- собственно твердая фаза (скелетная часть). Если радиоактивность твердой фазы^м постоянна, ее плотность бм = бгл, а ^ж = 0, то I Ряд урана-радия JjlL 0 5 Г i г г 2,5 j Знергия у-лучеИ 10 , Дж Рис. 37. Спектры 7"лУчей радио- радиоактивных минералов Зависимость вида A1.74), получаемая для конкретных типов пород, устанавливает практически однозначную связь между qn и &Гл и широко используется при интерпретации результатов гамма- каротажа скважин (рис. 38). В ряде случаев в формировании есте- естественной радиоактивности терригенных пород наряду с глинистой существенную роль может играть и мелкоалевритовый компонент породы (рис. 39). В связи с этим радиоактивность песчано-алеврнто* 68
ШСТЬ1Х пород не всегда является однозначной функцией глини- глинистости и определяется выражением а 7п — ^пес^пес I #ал«ал г <7гл 5 ?п'-Кг (И 75) — г- • 1 • * • г* > 1 7,0 оО/° S/ j L рис. 38 Зависимости радиоактивности пород qn от глинистости /ггл (а) и со- содержания нерастворимого остатка Сно (б) ^ — посчаио-глиннстые породы [2], б —^карбонатные отложения верхнего мела Восточ- Восточного Предкавказья: / — мергели и глины, отобранные из пропластков и трещин в извест- известняках, 2 — известняки Рис. 39. Гистограмма изменения удельной радиоактивности фракций альб-аптских пес- чано-алеврито-глинистых пород Прикумского нефтегазоносного района Ставропольского края 4^ — диаметр фракции <0,01 0,05 0,1 0,2 0,063 0,18 0,3 If, м _| 1 Z N \ 6- / If / 8 Рис. 40. Блок-схема лабораторной радиомет- радиометрической установки со сцинтилляционным детектором. / _ детектор излучения, 2 — усилитель; 3 — дис- криминатор, 4 — счетчик импульсов; 5 — источник высоковольтного стабилизированного питания де- текюра; 6* — источник низковольтного питания; 7 — защитный свинцовый домик; 5— внешний блок цифропечати гДе $ше, ^ал, ^7гл — удельные радиоактивности соответственно песчаной, алевритовой и глинистой составляющей твердой фазы породы; kUQC, ?ал, &гл — их объемные содержания. 69
Определение естественной радиоактивности пород Естественную радиоактивность пород qyn определяют путем измерения интенсивностей у-излучения образца 1уп и эталона 1уш известной радиоактивности qy3T совместно с фоновой /Тф, обус- обусловленной космической радиацией, собственным излучением ма- материалов детектора и окружающих предметов. Величину /7ф ущщ тывают путем самостоятельного ее измерения. Удельную массовую радиоактивность породы рассчитывают по формуле п"qyЭТ где Мп и Мэт — массы породы и эталона. Формула A1.76) справедлива при идентичности плотности спек- спектра излучения и геометрических размеров образца и эталона [10], В качестве эталонов обычно применяют порошковые препараты равновесного радия активностью 60—100 пг-экв Ra/r с наполните- лем из неактивного кварцевого песка (при изучении терригеннщ пород) или кальцита (при изучении карбонатных пород). Наряду с радиевым могут применяться и эталоны равновесного урана, то- тория или калия. Блок-схема одной из установок со сцинтилляционным детекто- детектором (счетчиком), служащих для измерения интенсивности -у-излу- чения образцов пород, показана на рис. 40. Исследуемый препарат (образец породы, эталон) помещают в плексигласовый стакан, устанавливаемый непосредственно на сцинтиллятор. Для снижения влияния внешней радиации детектор помещают в защитный свин- свинцовый домик 7. Регистрацию импульсов, получаемых на выходе с детектора, производят с помощью счетных устройств, включаю- включающих усилитель 2, дискриминатор 3 и счетчик импульсов 4. Совре- Современные счетные устройства (например, ПСО2-2ем и др.) позволяют измерять среднюю частоту следования импульсов, производить их счет в течение заданного времени в автоматическом режиме с за- записью числа зарегистрированных импульсов во внешнем цифропе* чатающем блоке типа Б3-15м. Питание измерительной схемы уста- установки осуществляется от блоков высокого 5 и низкого 6 напряжений, В качестве источника высокого напряжения могут использоваться высоковольтные стабилизаторы типа ПВ-2-2, ВС-2-2 или блоке БНВ-06 и другие, а низкого — блоки типа БНН-30, БНН-33, БНН-42, B5-47ji др. Ь'- Определение5^ естественной]! радиоактивности пород проводят в следующем^порядке. 1. Образец сухой породы измельчают до размеров частиц, мень- меньших" 0,25 мм, и находят его массу Ми. На практике наиболее часто применяют детекторы, позволяющие измерять у-активность образ- образцов пород массой 30, 50 или 150 г. 2. Подготавливают установку к измерениям, для чего: 70
\ подключают детектор к источникам питания и счетному устройству; б) включают в сеть источники питания и счетное устройство; в) подают на схему детектора необходимое низкое и высокое заражения, величины которых зависят от типа и характеристики детектора; г) включают счетное устройство и дают установке прогреться тече1Ше 30—40 мин, после чего проверяют правильность ее ра- 3, Измеряют поочередно в течение фиксированного времени t число импульсов, поступающих от: а) эталона и фона #эт+ф; б) образца породы и фона Л^п+ф; в) фона Л?ф. Время измерения t выбирают таким, чтобы погрешности из-за статистических флуктуации были минимально допустимыми. Кроме того, для снижения погрешности число импульсов Л^эт+ф, Мп+ф II j\L измеряют по несколько раз и устанавливают их средние зна- значения. 4, По результатам измерений рассчитывают: а) регистрируемые интенсивности у-излучения h п*Ф = ^п+ф/^ h эт+Ф = Л/эт+ф/fc h Ф = ЩН> (П. 77) б) интенсивности, соответствующие 7-активности образца по- породы и эталона, свободные от фонового излучения 1у и = fy п+ф — ^7Ф» 1у эт = 1у эт+ф — ^7Ф' A1.78) в) удельную массовую Y-активность породы по формуле A1.76). При детальных исследованиях пород изучают не только общую у-активность, но и содержание в них радиоактивных элементов, таких, как уран, торий, калий fgK и др. Концентрация радиоактив- радиоактивных элементов в породах определяется способом гамма-спектро- гамма-спектрометрии» сущность которого заключается в сравнении при одинако- одинаковых условиях измерений спектра ^-излучения (распределение у-из- лучения но энергиям) исследуемой породы со спектрами эталонов равновесного урана, тория, калия и др. Спектр у-излучения определяется с помощью гамма-спектро- гамма-спектрометра по распределению амплитуд сигналов, поступающих от спек- спектрометрических детекторов. Последние (сцинтилляционные или полупроводниковые) преобразуют энергию регистрируемых у-кван- тов в электрические импульсы, амплитуда которых пропорцио- пропорциональна анергии, а количество — потоку 7~квантов. Регистрация шшштудного распределения импульсов производится многока- многоканальными амплитудными анализаторами типа АИ-100, АИ-128, АИ-256 и др. 71
Искусственная радиоактивность (нейтронные свойства ropHt*g пород) ; Искусственная радиоактивность связана с радиоактивным падрм искусственных радиоактивных изотопов химических зе% ментов, образующихся при облучении их элементарными части- частицами (электронами, протонами, нейтронами, 7"частиДами и дрЛ в результате изменений в ядре, происходящих вследствие пронад новения в него заряженной частицы или нейтрона. ; Нейтрон — электрически нейтральная ядерная частица с сой (Мп = 1,0086654-10~24 г), примерно в 1836 раз большей esjj электрона (позитрона) и приблизительно равной массе протощ (ядра водорода). Так как он представляет собой нестабильную ядерную частицу, то распадается с TV2~ 16,83 мин на протон, электрон и антинейтрино с выделением энергии 0,78-10~~13 Д^ Нейтроны не взаимодействуют с электронными оболочками атомов и не отталкиваются кулоновским полем ядра, что обусловливав их высокую проникающую способность. По величине энергии раз* личают нейтроны: холодные — 10~21 Дж; тепловые — 25-10~21 Цщ медленные — 5-10~~20 Дж; надтепловые—0,3—5-Ю"8 Дж, ре* зонансные — 10~17 Дж, промежуточные @,5—2)-10~14 Дж и бы- быстрые — 2 • Ю-14 — 2 • 10-12 Дж. Нейтроны, получаемые при помощи нейтронных источников^ распространяются в окружающей среде и взаимодействуют с яд- ядрами ее химических элементов. При этом наиболее существенный процессами являются рассеяние и поглощение (захват). : Рассеяние нейтронов может быть упругим и н е у п р у *. г и м. Сущность этого процесса состоит в изменении направления движения и уменьшении кинетической энергии нейтронов при щ столкновении с ядрами элементов окружающей среды. При уяру»; гом рассеянии происходит перераспределение энергии между щ*: летевшим нейтроном и неподвижным ядром в соответствий с hi массами и углом рассеяния по принципу соударения упругих ша- шаров. При этом внутреннее состояние ядра и кинетическая энергия системы нейтрон—ядро остаются неизменными. Энергия Е2 и ош>, рость v2 нейтрона после столкновения с ядром определяются во формулам v: (Л . где Ех и v± — соответственно энергия и скорость нейтрона до столк- столкновения; А — массовое число ядра элемента, на котором пронс: ходит рассеяние нейтронов; 6 — угол между начальным и послед дующим направлениями движения нейтрона. > А 72 Л
Величина потери энергии нейтроном зависит от характера его кновения с ядром и массы последнего. Максимальная потеря Нисходит при центральном столкновении F = л), когда масса ipa равна или соизмерима с массой нейтрона. Так, при централь- ом соударении последнего с ядром водорода (А = 1) нейтрон мо- Н потерять всю свою энергию в одном акте. Следовательно, во- д является аномальным замедлителем нейтронов. При неупругом рассеянии нейтрон сначала захватывается, а затем выбрасывается ядром, но уже с меньшей энергией и под не- некоторым углом к направлению начального движения. Ядро же, захватившее и потерявшее нейтрон, остается на некоторое время з возбужденном состоянии и затем возвращается в основное, испу- испуская т-квант: r\e z% n z^* — ядРа ИСХ°ДНОГО элемента и того же элемента в воз- возбужденном состоянии; In и In — нейтроны поглощенный и выбро- выброшенный ядром. Быстрые нейтроны, распространяясь в окружающей среде, 1% процессе неупругого и упругого рассеяний сравнительно быстро (за 10~*—Ю~5 с) теряют свою энергию (до 25-10~21 Дж) и скорость (до 2200 м/с) и превращаются в тепловые. Последние поглощаются ядрами вследствие реакции радиационного захвата с образованием на первой стадии составных ядер, которые затем (черед 10~~12— 10~~*4 с) переходят в основное состояние с испусканием у-квантов: Наиболее типичной реакцией (я, у) является реакция захвата теплового нейтрона водородом Распределение нейтронов в среде (породах), т. е. плотность нейтронов на различном расстоянии от источника, зависит от ней- нейтронных свойств этих пород, в основном связанных с химическим составом последних. Для большинства горных пород поглощающие и замедляющие свойства определяются водородосодержанием: чем оно выше, тем быстрее убывает плотность нейтронов с удалением от источника. Вероятность той или иной реакции взаимодействия нейтронов с веществом количественно характеризуется нейтронным эффективным сечением, численное значение кото- которого выражается в м2. Эффективные сечения процессов рассеяния и захвата, относя- относящиеся к одному ядру, называют соответственно микроскопи- микроскопическими ядерными сечениями рассеяния ар и з а х в а т а а3' элемента. Эффективные же сечения рас- рассеяния и захвата, относящиеся к единице объема вещества (породы), 73
называют макроскопическими сечениями рас' сеяния У-р и захвата 2з- Для вещества, состоящего из разных химических элементов (П.81 п 2Я=- 100 Ac где Лг и Р( — атомная масса и содержание (в %) t-ro элемента; Gpi и сг31- — соответственно эффективные ядерные сечения paccet* 0,8 0,6 0,2- t\ :\ 70 20 JO 50 А Рис. 41. Зависимость средней логарифмической потери энергии | нейтроном от атомной массы элемента А *> ния и захвата i-то элемента; б — плотность вещества; Л/^ — число Авогадро. Макроскопические сечения рассеяния недостаточно полно хв/; рактеризуют свойства пород замедлять нейтроны. Более полно оеш1 оцениваются комплексными параметрами — макроскопи-' ческой замедляющей способностью пород н коэффициентом замедления. Первый параметр ра-} вен произведению макроскопического сечения рассеяния 2Р на среднюю логарифмическую потерю энергии ? при одном соударе^ нии. Последняя зависит от атомной массы элемента (рис. 41), с кхИ торым сталкивается нейтрон \ 6 = In- ¦ = 1 + (А + IJ 2А In А — ] А + 1 AШГ где Е± и Е2 — энергии нейтрона до и после соударения. Коэффициент замедления представляет собой отношение ляющей способности породы к ее макроскопическому сечению заж-| 74 ^
v В процессе замедления нейтрон совершает сложный путь, - оящий из отдельных прямолинейных участков, представляю- х собой отрезки его пути между последовательными соударе- ялш с ядрами элементов среды. Эти отрезки имеют различную Ш1*шУ и направление и определяют траекторию движения нейтрона источника до точки замедления, т. е. места, где он становится епловым. Расстояние, на которое смещается нейтрон в процессе замедления, оценивается средним квадратом r% 3 смещения по "и длиной замедле- а $ прям пня U, связанных между собой "соотношением: 'ср. з (II 84) Длина замедления L3 характеризует интенсив- интенсивность замедления нейтро- нов и зависит в основном от содержания в породах водорода (величины пори- пористости) [рис. 42, а]. Вре- зия с момента излучения нейтрона (Еп 0 = 2 • 10~ — -—2- Ю2 Дж) до момента его превращения в тепло- тепловой (?,1т = 25.10-21 Дж) называется временем за- замедления, которое рассчи- рассчитывается по формуле 6 2 0 \ _\ ~2 - к 1 1 г 20 кп°/о 0 20 кПУ% Рис. 42. Зависимости длины Гзамедления L3 (а) и среднего времени жизни т^ теп- тепловых нейтронов (б) от коэффициента по- пористости kn для песчаников (/) и извест- известняков B) [по В. В. Ларионову]. Минерализация пластовых вод 7 — 0, II — 250 г/л NaCI (II 85) где Е„о и нейтронов; Еп, М- - соответственно начальная и тепловая энергии абсолютная масса нейтрона A,667-10"4 г); 4р р-^ i/jVp — средняя длина свободного пробега нейтрона до замедления. Нейтроны, достигшие теплового состояния, в процессе их диф- диффузии в конечном счете поглощаются (захватываются) элементами с высоким эффективным сечением захвата. Поглощающие свойства пород характеризуются временем (временем жизни теплового ней- нейтрона) T<f, длиной Ld и коэффициентом D диффузии. Под временем жизни теплового нейтрона понимается время су- существования нейтрона в породе с момента его превращения в теп- тепловой до момента его захвата ядром и определяется по формуле ^ = -^-2-=-^, A1.86) 75
где vr — скорость тепловых нейтронов B200 м/с); /ср 3 = 1/'2з <~~~ средняя длина свободного пробега нейтрона до захвата. Длина диффузии характеризует среднюю длину пути нейтрона от точки, где он стал тепловым, до точки захвата: (Пщ где ГсР а — средний квадрат расстояния по прямой от точки за- замедления нейтрона до точки его захвата ядром. Отсюда коэффициент диффузии L2 Интенсивность поглощения тепловых нейтронов зависит от со- держания в породах элементов с высоким эффективным сеченйед захвата, основным из которых в осадочных породах является хлор (рис. 42, б). Замедляющая и поглощающая способности горных пород определяют пространственное распределение нейтронов на различных стадиях их взаимодействия с породами, на изучение которого основаны нейтронные методы исследования скважин, Определение нейтронных параметров пород Нейтронные параметры горных пород можно устанавливать как экспериментальным, так и расчетным путем. В петрофизике макро- макроскопические сечения рассеяния 2р и захвата (поглощения) 2&* длину La, коэффициент D и время т^ диффузии тепловых нейтро- нейтронов получают с помощью расчета. Макроскопические сечения рассеяния 2р (см™1) и захвата ^ (см) нейтронов в породе вычисляют по формулам A1.81), A1.82), с учетом объемного соотношения твердой фазы и заполнителя пор, определяемого коэффициентом пористости &п, их плотностей, хи- химического состава, а также микроскопических сечений ядерных реакций с химическими элементами породы. Химический состав твердой фазы и заполнителя пор породы устанавливают по результатам их химических анализов, значения микроскопических сечений рассеяния ар и захвата а3 для отдель- отдельных элементов (окислов), входящих в состав породы,— по табл, 5. Затем, используя известные 23, 2р и vt ~ 2,2-105 см/с рассчиты- рассчитывают: 1) время диффузии (жизни) тепловых нейтронов xd (в с) по формуле A1.86): 2) коэффициент диффузии (в см2/с) по формуле D = Щ , (II Щ 2Sp(l-cose) где cos 8 — средний косинус угла рассеяния нейтронов, принимае- принимаемый для воды и водосодержащих пород равным 0,2; 76
Таблица б Значения микроскопических сечений рассеяния ор и захвата о3 нейтронов для ядер некоторых элементов и окислов [10] Элемент (окисел) Н С О К а Mg S к Сг Fe С1 МаО8 CO., Саб Fe2O3 Н2О к2о MgO МпО2 Р9О5 SiO2 Na2O Атомная (молекуляр- (молекулярная) масса 1 12 16 23 24 32 39 40 56 35 101,9 44,01 56,08 159,7 18,02 78,2 40,32 70,93 142,0 60,06 61,99 Сечение, . Ю8 м2 38±4 4,8+0,2 4,2+0,3 4,0+0,5 3,6+0,4 1,1+0,2 1,5+0,3 3,2+0,3 11 + 1 16+3 15,4 13,2 13,2 34,6 44,92 7,2 7,8 6,5 42,0 10,1 12,2 0,328+0,002 C,78+0,07). Ю-3 < 2-10-4 0,515+0,008 0,063+0,003 0,52+0,02 2,07+0,07 0,44+0,02 2,62+0,06 33,8+1,1 0,43 0,0049 0,43 4,86 0,66 3,94 0,059 12,6 0,38 0,13 1,0 3} длину диффузии Ld (в см) по формуле A1.87). Нейтронные параметры пород обычно рассчитываются для стан- стандартных температурных условий (То = 293 К) [10]. § 8. УПРУГИЕ СВОЙСТВА ГОРНЫХ ПОРОД Компоненты горных пород — твердая фаза, жидкость и газы — обладают резко отличными упругими свойствами. В связи с этим последние определяются упругими свойствами составляющих ком- компонентов и фазовой однородностью. Условно различают иде- идеально и дифференциально упругие горные по- породы. Упругие параметры В породе, на которую действуют внешние силы, возникают внут- внутренние силы, стремящиеся к восстановлению ее начальной формы п размеров. Величина последних сил, приходящаяся на единицу площади сечения тела, измеряемая в паскалях, называется н а - пряжение м,1 которое является векторной величиной, зави- * Н пористых породах напряжения концентрируются в основном в об" даст и контакта минеральных зерен. 77
сящей от действия внешних сил, внутренних свойств и формы об- образца породы. В зависимости от характера приложенных внешня* сил, образец породы может испытывать линейное, плоскостное и объемное напряженные состояния. Под воздействием внешннд сил изменяются линейные размеры, объем или форма горной по- породы. Эти изменения называются деформацией. Деформации, соответствующие нормальным напряжениям, вы- выражаются через относительное изменение е линейных размеров образца (рис. 43) к к в,--^=-^-. (Н.9!) где 1г и /2 — начальные продольный и поперечный размеры об- образца; 1[ и й — то же, при одностороннем сжатии. Рис. 43. Пример простого одноосного сжатия бруска породы а — брусок породы до сжатия; б — тот же брусок после де- деформации под воздействием напряжения р Деформации 8С, соответствующие касательным напряжениями выражаются через тангенс угла сдвига а-граней образца: ес = tgou При увеличении напряжений можно наблюдать три вида де- деформации породы — упругую, пластическую и разрушающую* Для каждого из приложенных напряжений существует свой ко- коэффициент пропорциональности между напряжениями и упругими деформациями, являющийся упругим параметром породы. Коэф- Коэффициент пропорциональности между продольным (сжимающим или стягивающим) напряжением р и соответствующей ему относитель- относительной деформацией е называется модулем упругости или модулем Юнга Е: р^Ее. {IIЩ. Коэффициентом пропорциональности между касательным на- напряжением тк и соответствующей деформацией сдвига гс является 78
сдвига G: При объемном напряженном состоянии породы, что соответст- соответстт действию всестороннего гидростатического давления, связь ЖДУ величиной напряжения и относительным изменением объема \yIV выражается через модуль всестороннего сжа- сжатия Кс или сжимаемость рс: (II 94) Связь между относительными продольными и поперечными де- сЬормашшш! сжатия (растяжения) устанавливается коэффициентом Пуассона Скорость распространения упругих волн Если к ограниченному участку упругой породы приложить внешние силы, которые вызовут напряжения, не превышающие пре- предел упругости, то на этом участке возникнут деформации и прои- произойдет смещение частиц вдоль направления действующей силы» Смещение одной частицы вызывает сдвиг других более удаленных. Распространение упругой деформации проис- происходит с определенной ско- скоростью. Если на породу дей- действуют знакопеременные или кратковременные силы, то впей возникают упругие ко- _ лл п «г гт J Рис. 44. Схема смещения частиц среды лебания. Процесс последова- прп распространении продольных (а) и тельного распространения поперечных (б) волн В ПОрОДе Деформаций (упру- Направление движения. / — волн, 2 — час- хик колебаний) называется тиц среды упругой волной. В за- зависимости от вида деформации в породе возникают различные типы волн, основными из которых являются продольные и поперечные. Первые связаны с объемной деформацией среды. Распростране- Распространение продольной волны представляет собой перемещение зон растя- растяжения и сжатия, при котором частицы среды совершают колебания около своего первоначального положения в направлении, совпа- совпадающем с направлением распространения волны (рис. 44). Про- Продольные волны распространяются в любой среде — твердых те- 79
лах, жидкостях и газах, так как все вещества обладают сопротив- сопротивлением объемному сжатию. Поперечные волны обусловлены деформациями сдвига в среде и присущи только твердым телам. Последнее связано с тем, цщ в жидкостях и газах отсутствуют сопротивления сдвигу. Их рас» пространение представляет собой перемещение зоны скольжения слоев среды относительно друг друга; частицы среды совершавд колебания около своего первоначального положения в плоскости перпендикулярной к направлению распространения волны (сц, рис. 44, б). Одним из важных кинематических параметров упругих волн является скорость их распространения. Для идеально упругих т т 7 i Рис. 45. Схематическая модель неглинистой (а) и глинистой (б) пород изотропных горных пород скорости продольных vP и поперечных волн определяются по следующим формулам: v) A — 2v) ' Sn(l +v) (IL96) A1,97) где 6П — плотность породы; Е и v — соответственно модуль Юнп и коэффициент Пуассона. Для горных пород величина Е изменяется от 0,15-10~4 до 0,6- 10~5МПа, а коэффициент поперечного сжатия v близок к 0,25* После подстановки средних значений упругих констант получш Vplvs ~ 1,73. Следовательно, продольная волна, распространяю- распространяющаяся приблизительно в 1,75 раза быстрее поперечной, достигает удаленной точки раньше. i-j Осадочные горные породы в большинстве своем являются даф* ференциально упругими и не обладают достаточно совершенной связью между фазами. Вследствие этого скорости распространения 80
иХ упругих волн отличаются от скоростей, вычисляемых по фор- в ам л.96 и 11.97. Известен ряд уравнений, с помощью которых ^танавливается зависимость скорости продольных волн в породах * т скоростей их распространения в отдельных фазах и коэффици- коэффициента пористости [2, 5]. Наиболее простым из них и получившим широкое применение на практике является уравнение Вилли (урав- (уравнение среднего времени), согласно которому vn рассчитывается лорнстой породе по времени ее прохождения через минеральный скелет tu и жидкость *ж, заполняющую поры (рис. 45): *п = *м + *ж (Н.98) пли 1 — Лгп A1.99) О *tO 80 /20 f?OT,°CO 20,0 Рис* 46. Оценка интерваль- интервального времени А^ск в мине- минеральном скелете породы по данным А^п и kn Рис. 47. Номограмма для определения скоро- скорости распространения продольных волн в жид- жидкости vB при заданных минерализации Св, давлении рил и температуре Т [5] Пунктиром показан пример определения v при где ku — коэффициент пористости; vM и иж — скорости продоль- продольных волн соответственно в минеральном скелете и насыщающей жидкости. В уравнении A1.99) вместо скорости практичнее исполь- использовать соответствующее ей интервальное время (величину, обрат- обратную скорости) А^п = A -- ku) Мм + къЫж, (П ЮО) где Д/м, A tn и А/ж — интервальные времена соответственно в ми- минимальном скелете, породе и насыщающей ее жидкости. Величина A tM зависит от минералогического состава скелета и для конкретных типов отложений является постоянной. В поро- породах с мономинеральной твердой фазой она соответствует интерваль- интервальному времени распространения упругих волн в породообразующем минерале (кварце, кальците, доломите и т. д.). При содержании в скелете породы нескольких минералов, отличающихся по своим 81
упругим свойствам, Л/м определяется как средняя взвешенная' величина. Среднее значение скорости распространения волн в оса- осадочных породах составляет 2500—4000 м/с (см. табл. 7). Основными факторами, влияющими на скорость распростране- распространения упругих колебаний в горных породах, являются: литолого* минералогический состав, поровое пространство, заполненное жщ* костью (форма и размер пор играют меньшую роль), степень насы- насыщения пор жидкостью или газом, степень цементации, текстурные и структурные особенности, разность горного и пластового давле- давлений (эффективное давление) и др. Часто возникает необходимость в определении Д/м для конкрет- конкретного интервала геологического разреза. В этом случае сопостав- сопоставляются времена, отсчитанные по диаграмме акустического каро- каротажа А/п, со значениями пористости ku, установленными по керщу или одному из геофизических методов. Полученные данные исполь- используются для построения графика зависимости Д/п от kn (рис. 48), Осредненная прямая, проведенная через нанесенные точки, отсе- отсекает на оси времени значение AtM при ku = 0. Если пористость па разрезу изменяется слабо, то значение А4, для каждого однород* ного пласта рассчитывают по формуле 1 — kn Скорость (интервальное время) в воде зависит от ее минерали- минерализации, температуры и давления и по данным из работы [5 ] опреде- определяется с помощью номограммы (рис. 47) или эмпирического выра- выражения Д*ж = [1557 — 0,0245 G4 - ТJ + 0,8С + 1,9/?пл]-1.106, A1.10$ где Т — температура воды (в °С); С — минерализация (в г/л); /?пл — пластовое давление (в МПа). Скорость распространения упругих волн в нефти и газе меньше,- чем в воде. Это объясняется большей сжимаемостью углеводоро- углеводородов, чем воды. Так, скорость распространения волн в песке, пол- полностью насыщенном нефтью, на 15—20 % меньше, чем в песке* заполненном водой. Выражения A1.99) и A1.100) справедливы для однородных неглинистых пород, обладающих межзерновой пори- пористостью. На рис. 48 показана зависимость интервальных времен и скоростей от коэффициента пористости для пород различного литологического состава. На величину скорости (как было сказано выше) оказывает влияние также тип цемента, который принято делить на вязкий (глинистый) и жесткий (карбонатный, кварцевый и пр.). Увеличе- Увеличение количества жесткого цемента соответствует росту доли твердой фазы в единице объема среды, уменьшению пористости, росту мо- модуля упругости н повышению скорости распространения волн* При возрастании объемной доли глинистого цемента, обладающего высокой сжимаемостью и пластичностью, наблюдается снижение 82
Съемного модуля упругости среды и скорости. При значительной ечнчине krn для "определения kn используют обобщенное уравне- уравнение среднего времени для пород с любой глинистостью: А^п = (I —kn — ктл) А?м + бглА^гл + ЬиЫж, A1.103) -у интервальное время в пластовых глинах. В случае со- содержания глин в рассеянном виде величина Д?гл берется меньшей to сравнению с А^гл в пластовых глинах; для приближенных рас- расчетов применяется А^гл = 500 мкс/м. Наряду с уравнением сред- среднего времени применяется и более общее — уравнение степенной связи [21, имеющее следующий вид: А/ = Л* + (Лгж - Д/м) С + Кл ~ А/м) СЛ- (".104) W- Atn, мкс/м Рис. 48. Зависимости интервального времени про- пробега упругой волны А^п от коэффициента пористости kn пород vM (в м/с). /, 2~доломиты — 7000 — 8000, известняки — 6400 — 7000, 5 — 5 ¦ Здесь показатели степени тп и тгл зависят от структуры и сте- степени консолидации породы (обычно 0,7—1,5) и возрастают с ее увеличением. При шп = 1 и krjl = 0 уравнение степенной связи преобра- преобразуется в уравнение A1.100), а при тп = 1 и тгл = 1 — в уравне- уравнение A1.103). Распространение упругих волн в горных породах сопровож- сопровождается постепенным уменьшением их интенсивности по мере удале- удаления от источника возбуждения. Уменьшение интенсивности упругих колебаний в основном свя- связано с поглощением части энергии упругих колебаний породой и превращением ее в тепловую вследствие взаимного трения частиц породы, совершающих колебательные движения; с рассеиванием якуаической энергии и неодиородностями породы. 83
Амплитуда А упругих колебаний связана с пройденным волной расстоянием / экспоненциально: Л = Лоехр( — аАК/), (ПЛ«) где Л о — амплитуда упругих колебаний вблизи источника дения; аАК — коэффициент поглощения упругих волн. Коэффициент поглощения упругих во#н аАК характеризует интенсивность поглощения энергии волн в среде и может быть определен по формуле аАк^ ln ln ~J^> ШЩ где Аг и Л 2 — амплитуды волн, регистрируемые приемниками расположенными на расстоянии / друг от друга. Размерность осАК выражают в децибелах на 1 м или м", от- относя величину ослабления к единице длины. Перевод величины затухания, измеренной в дБ/м, в м, производится по формуле аАК [дБ/м] - 8,68адк [м-*] Величина аАК в горных породах зависит от монолитности их скелета, пористости, трещиноватости, вещественного состава за- заполнителя пор и других параметров и вместе со скоростью широко используется при изучении пород. Поглощающие свойства пород связаны с литологией еще более тесно, чем скорость распространения упругих волн. Интенсивность поглощения породой упругих колебаний зависит также от харак- характера жидкости, заполняющей поровое пространство. В слабосце- ментированных нефтеносных и газоносных породах с хорошей по- пористостью затухание колебаний происходит более интенсивно, чет в таких же породах, но водоносных. Это особенно заметно в газо- газоносных породах из-за большой разницы скоростей распространения упругих волн в минеральном скелете породы и заполняющем но- норовое пространство газе. Наибольшее затухание претерпевают уп- упругие волны в трещиноватых и кавернозных породах. В связи е этим акустический каротаж по затуханию весьма эффективен при изучении разреза скважин, вскрывающих карбонатные породы. На скорость и поглощение упругих волн в породах существенно влияет эффективное давление. Его увеличение приводит к сжатию скелета породы, уменьшению ее пористости и увеличению кон- контактной упругости между зернами, а следовательно, к возраста- возрастанию скоростей упругих волн и уменьшению коэффициентов их по- поглощения. С увеличением эффективного давления скорость в породах воз- возрастает сначала быстро. По мере сближения и переукладки зерен породы под действием давления и улучшения акустического кон- контакта между ними рост скорости замедляется. Дальнейшее ее уве- увеличение определяется контактной упругостью зерен и при давле- давлении, превышающем некоторое предельное, растет очень медленно* 84
Поя низком давлении даже незначительное содержание газа в жид- vOM заполнителе пор приводит к резкому уменьшению скорости vripvrux волн в пласте. С ростом давления наблюдается ее плавное увеличение. В воде, находящейся под давлением около 60 МПа, ?К0роеть увеличивается примерно на 5—7 % по сравнению со скоростью в воде, находящейся при атмосферном давлении. Определение скорости распространения упругих волн Скорость распространения упругих продольных ультразвуко- ультразвуковых волн vP в образцах горных пород обычно определяют одним из импульсных ультразвуковых способов: прямого прозвучивания, продольного профилирования, кратных отражений или критиче- критических углов [10]. На практике наиболее широкое применение из них получил способ прямого прозвучивания, при котором с по- помощью специальной аппаратуры измеряется время t пробега уп- упругой ультразвуковой волны через образец породы длиной L. По данным этих измерений рассчитывают скорость vp^LH (П. 107) Скорость vP, измеренная на образце породы, соответствует ско- скорости г>рм в массиве при условии, если г/Я>1 (где г — радиус образца; Я — длина ультразвуко- ультразвуковой волны); при г/Я<0,3 устанав- устанавливается скорость vP cT в стер- Рис. 49. Блок-схема ИСУ-1 жяе породы [8]. Величины vP M 1\~\п и 1>рст получают на образцах раз- j Г лого диаметра с использованием иьезодатчиков с различной часто- частотой, которые обеспечивают соот- соответствующие отношения г/К. Ско- Скорость v^ определяют с помощью специальных датчиков, предназна- предназначенных для этой цели. Для прямого прозвучивания образцов пород применяют аппара- аппаратуру типа ИПА, ДУК-20, УЗИС, УК-ЮП, ИСУ-1 и др. Из пере- перечисленной аппаратуры измеритель скорости ультразвука ИСУ-1, разработанный во ВНИИГеофизике, позволяет определять ско- скорость ультразвука в образцах пород с наименьшей погрешностью A %). Аппаратура ИСУ-1 (рис. 49) состоит из устройства 1 для прозвучивания образцов, генератора 2 высоковольтных импульсов и комплекта стандартных приборов, включающего усилитель 3 (ламповый вольтметр ВЗ-4), импульсный осциллограф 4 типа С1-20, двухканальный генератор импульсов 5 типа Г5-4Б, источ- источники питания 6 и 7 типов УИП-1 и УИП-2. Основными частями устройства для прозвучивания 1 являются пьезоэлементы (излуча- (излучатель И и приемник Я), в качестве которых используются столбики 85
керамики ЦТС-23 высотой 3 и 15 мм, служащие соответственно для получения высоко- и низкочастотных колебаний. При измерениях с помощью аппаратуры ИСУ-1 образец породы О цилиндрической формы размером L — 20—40 мм и d = 30— 40 мм, насыщенный пластовой жидкостью, устанавливается в ус? ройстве для прозвучивания между излучателем И и приемником Я При этом на экране осциллографа наряду с горизонтальным еле- дом электронного луча появляются и первые периоды электриче- электрического сигнала от приемника ультразвуковых колебаний. Для улу*ь шения акустического контакта между пьезоэлементами и торцама образца последние смачиваются насыщающей жидкостью. Замер времени t пробега упругой волны через образец породи производят дважды. Расхождение в отсчетах допускается не более 0,02 мкс. При регистрации первого вступления его амплитуда на экране с помощью усилителя осциллографа и лампового вольтметра ВЗ-4 устанавливается равной 1 или 2 см. Измеряют штангетщр- кулем или микрометром расстояние между пластинами устройства для прозвучивания LnjI, соответствующее длине образца Lo6p, В случае прозвучивания жидкости между излучателем И и прием- приемником П помещают специальный стакан, заполненный исследуе- исследуемой жидкостью. Измерение величин t и LB производится так же, как и при прозвучивании образца породы. Время регистрации первого вступления включает время по* стоянной задержки в аппаратуре и время прохождения упругого импульса через образец породы (жидкости) [табл. 6]. Для исклш* чения постоянной временной задержки определение величины % проводится по методике относительных измерений. При этом от- относительные измерения могут выполняться двумя способами. Пер- Первый способ заключается в последовательном измерении времени задержки tx и t2 на образцах различной длины L01 и L02, приготов- приготовленных из одной породы: A1ЛЩ U-U Таблица 6 Скорость ультразвука vpB в дистиллированной воде при различных температурах т, °с 15 16 17 18 19 20 Орв, М/С 1466,75 1469,70 1473,07 1476,35 1479,55 1482,66 т, °с 21 22 23 24 25 vPb, м/с 1485,69 1488,63 1491,50 1494,29 1497,00 т, °С 26 27 28 29 30 vPb, м/с 1499,64 1502,20 1504,68 1507,10 1509,44
при втором способе последовательно измеряют времена за- еожкн, помещая между излучателем и приемником сначала об- образец породы — tl9 а затем дистиллированную воду — t2. vp = — ^ , A1.109) гяе Lu и Ьъ — длина образца породы и толщина слоя воды; vP в — fKOPOCTb продольных ультразвуковых волн в дистиллированной воде при температуре измерений (см. табл. 6). § 9. ДРУГИЕ СВОЙСТВА ГОРНЫХ ПОРОД При геофизических исследованиях скважин изучаются диэлек- хрические, магнитные и термические (тепловые) свойства. Диэлектрическая проницаемость Вещества, которые поляризуются в электрическом поле и об- *шдают очень малой электропроводностью (практически не прово- проводят электрический ток), называются диэлектриками. Поляризуе- Поляризуемость вещества характеризуется диэлектрической про- проницаемостью e=l+4jta, (II.110) где а — поляризуемость среды. Абсолютная диэлектрическая проницаемость определяется со- соотношением 8 — DIE, показывающим, во сколько раз напряжен- напряженность электрического поля Е в данном диэлектрике меньше напря- напряженности поля индукции D в вакууме. Диэлектрическая проницае- проницаемость является одним из физических свойств горной породы и по- показывает, во сколько раз возрастает емкость конденсатора, если вместо вакуума между обкладками в качестве диэлектрика поме- поместить данную породу. Диэлектрическая проницаемость измеряется в фарадах на метр и определяется в виде произведения (ПЛИ) где вот — относительная диэлектрическая проницаемость, пока- показывающая, во сколько раз абсолютная диэлектрическая проницае- проницаемость данной среды превышает абсолютную диэлектрическую про- проницаемость вакуума е0. Вакуум обладает наименьшей диэлектри- диэлектрической проницаемостью, равной 10~9/Зб я ^ 8,85-10~12 Ф/м. Различают следующие основные группы поляризации: смеще- смещения, ориентационная (релаксационная, дипольная) и структурная (объемная). Поляризация смещения вызвана элек- электронной, ионной и атомной поляризацией, когда под действием внешнего поля происходит смещение упруго-связанных зарядов. Ориентационная поляризация связана с на- наличием в диэлектрике полярных молекул. Она характерна для воды и других полярных жидкостей, в которых сопровождается ориен- 87
тацией диполей в направлении поля. Ориентационная поляриза- поляризация является наиболее значительной и снижается с повышением температуры, так как тепловое движение дезориентирует моле- молекулы в пространстве. Поляризации ориентационная и смещения протекают в течение 10~12—10~15 с и являются быстрыми. Структурная поляризация наблюдается в не- неоднородных, многофазных породах и обусловлена движением ионов накапливающихся на межфазовых границах раздела. Для прояв* ления структурной поляризации требуется больше времени, для других видов поляризации. 20 " П 07Н 100 50 0 hO 80к$у°/о Рис. 50. Экспериментальные зависи- зависимости диэлектрической проницаемости образцов пород еп. Отн от коэффициен- коэффициента водонасыщенности kB 1 — кварцевый песок; 2 — гранулярный из- известняк; 3 — глина; рв = 18 Ом«м; ku = = 37,5 % Рис. 51. Зависимость 8П отн образца -породы от коэффициента водо- kB и нефтенасыщенности ka (kn = 25 %) Рис. 52. Принципиальная схема измерения диэлектрической проницаемости резонансным способом с применением куметра. Г — генератор; С — настроечный конденсатор переменной емкости; Сх — измеритель- ный конденсатор с образцом породы, L и R — соответственно индуктивность и сопротив- сопротивление контура, V — электронный вольтметр, тА — миллиамперметр Диэлектрическая проницаемость горных пород зависит от со- става, содержания в них твердой, жидкой и газообразной фаз, я также от частоты поля и температуры. Для главных породообразую- породообразующих минералов значения 8ОтН невелики D— 10); для воды прш 20 °С еотн = 80, а при 100 °С — 55, поэтому диэлектрическая про- проницаемость пород в большей степени зависит от их водонасыщен* ности. Для нефти 8ОтН = 2—6, а для нефтенасыщенной породы соответственно 6—10. Для горных пород значение soth опреде- определяется разными видами поляризации. На рис. 50 показано увеличение диэлектрической проницаемо- проницаемости с ростом коэффициента водонасыщенности &в для различные 88
нов яород. Высокая относительная диэлектрическая проницае- т ть глин (до 50—60), возможно, объясняется наличием в них Рачительного количества слабосвязаннои^воды. Экспериментально ^становлено, что диэлектрическая проницаемость воды мало за- йнсит от минерализации и снижается с увеличением температуры, так как в этом случае растет тепловое движение молекул воды. Наибольший интерес представляет зависимость величины еп Отн от коэффициента водо- и нефтенасыщенности для пород-коллекто- пород-коллекторов (рис, 51). Как видно, величина 8П. Отн почти линейно связана с водонасыщенностью. Определение диэлектрической проницаемости пород Диэлектрическая проницаемость пород определяется мостовыми (при измерениях на частотах до 106 Гц) и резонансными (при изме- измерениях на частотах от 105 до 108 Гц) способами [10 ]. На практике ши- широкое распространение получили резонансные способы, одним из которых является способ с использованием куметра и измеритель- измерительного конденсатора постоянной емкости (рис. 52). При заданной частоте емкость конденсатора С изменяют до значения С2, настраи- настраивают колебательный контур куметра в резонанс с колебаниями ге- генератора и измеряют его добротность Qx*. Затем параллельно настроечному конденсатору переменной емкости включают измерительный Сх с образцом породы емкостью Cv При этом контур куметра расстраивается. Затем уменьшением емкости настроечного конденсатора до значения С3 повторно до- добиваются резонанса и измеряют новое значение добротности кон- контура Q2- После этого, проводя аналогичные измерительные опера- операции с конденсатором емкостью Сэт, в который вместо образца по- породы устанавливается эталон, определяют резонансные значения емкости настроечного конденсатора Сг и Сз. Эталон изготавливается из вещества с хорошо известной диэлектрической проницаемостью ?зт (обычно из полиэтилена с еОтН — 2,3 или полистирола с 80тн = =~ 2,2), форма и размеры его должны соответствовать форме и раз- размерам исследуемых образцов пород. По результатам этих измерений рассчитывают: 1) емкость (в Ф) 2) емкость (в Ф) ю-'2; (п.из) л Добротность Q контура характеризует отношение полной энергии, 'шпасеиной контуром, к энергии его джоулевых потерь и равняется 1/tg б (где б ~~~ угол потерь, показывающий потери электрической энергии, назы- называемые диэлектрическими). Добротность определяет, во сколько раз напря- напряжение на конденсаторе последовательного контура при резонансе больше Действующей на контур ЭДС. 89
3) геометрическую емкость измерительного конденсатора (в Щ Со = Сэте0/8Эт> (ПЛИ) где 8 0— диэлектрическая проницаемость вакуума, равная 8,885-10~12 Ф/м; 4) диэлектрическую проницаемость образца породы (в (П ЛЩ где h — толщина образца породы (в м); S — площадь поперечного сечения образца (в м2); ~ 5) тангенс угла диэлектрических потерь где величина (Сг—Со) близка к емкости СОбр образца породы, Магнитные свойства горных пород Основными магнитными параметрами горных пород, исполь- используемыми в геофизике являются: намагниченность /, магнитная восприимчивость х и магнитная проницаемость jut. Намагниченность характеризует магнитное состоя- состояние намагниченного тела и оценивается для однородного намагни- намагниченного тела, как отношение магнитного момента Мм к единице его объема V: I = MJV (ILII7) В случае неоднородного намагниченного тела / определяется для каждой точки (физически малого объема) и представляет сред- средний магнитный момент единицы объема (равный геометрической сумме магнитных моментов отдельных атомов и молекул, заклю* ченных в этом объеме). Единица намагниченности — 1 А/м, при которой 1 м2 вещества обладает магнитным моментом 1 А-м2. Маг- Магнитный момент, основная величина, характеризующая магнитные свойства вещества. Магнитным моментом определяются силы, дейст- действующие на тело во внешнем магнитном поле. М а г н и т н ая восприимчивость к определяет связь между магнитным моментом (намагниченностью) породы / и ее магнитным полем Я: к = 1Ш, A1ЛЩ где % — величина безразмерная; размерность Я—А/м. Различают объемную х магнитную восприимчивость, отнесенную к единице объема, и удельную худ, рассчитанную на 1 кг вещества. Магнитная проницаемость \i характеризует связь между магнитной индукцией В в породе и магнитным полем Я где \х0 — коэффициент пропорциональности, принятый в качестве магнитной постоянной. 90
Индукция В магнитного поля и его напряженность Н связаны вакууме соотношением В = [г0Я. Величины \х и х связаны соот- В шеннем \i = 1 + х> ^ является безразмерной величиной. В ва- JJvvMe х = 0 и ц = 1. Все породы и минералы по величинам к и \i и характеру взаи- одеЙствия с внешним магнитным полем разделяются на диамаг- * ||КП? парамагнетики и ферромагнетики. Для диамагнети- о в (^<0; ^<1) характерно намагничивание, которое прояв- проявляется в противоположном направлении относительно действую- действующего на него магнитного поля. К диамагнетикам относятся многие ггородообразующие минералы (кварц, кальцит, гипс, ангидрит, ортоклаз, циркон/ галит). Для галита х = —0,8-10~3 А/м. Из металлов висмут имеет наибольшее значение" к = 1,4-10~3 А/м. Для ""парамагнетиков A >%>0 и |i> 1) подмагничива- §ше во внешнем магнитном поле совпадает с его направлением. У па- парамагнетиков атомы, молекулы, ионы, ядра атомов обладают соб- собственными магнитными моментами, которые во внешнем магнит- магнитном поле ориентируются в зависимости от его направления дан- данного поля; если последнее отсутствует, то магнитные моменты ори- ориентируются хаотически и намагниченность / = 0. Следует также отметить, что парамагнетики притягиваются к полюсам магнита, а диамагнетики отталкиваются. К парамагнетикам относятся мно- многие металлы (Fe, Ni, К, Na, Co, Ca и др.), соли элементов группы железа и породообразующие минералы (слюды, ортоклаз, роговая обманка, хлорит, эпидот, гранат, турмалин, сульфиды — пирит, халькопирит и др.). У ферромагнетиков (х> 1 и \i ^> 1) магнитные мо- моменты атомов параллельны друг другу, формируя единый магнит. Намагничивание ферромагнетиков внешним полем Н происходит по замкнутой кривой, именуемой петлей гистерезиса. При снятии "поля намагничивания у ферромагнетиков сохраняется остаточная (спонтанная) намагниченность /г, имеющая важное значение при геологических исследованиях. Остаточная намагниченность ферро- ферромагнетика может быть снята воздействием намагничивающего поля противоположного знака, • напряженность которого Яс получила название коэрцитивной "силы. Для различных "ферромагнетиков характерны предельные температуры (точка Кюри)^выше^"которых они превращаются в парамагнетики; дальнейшее снижение темпе- температуры ведет к полной потере остаточной намагниченности. Отли- Отличительной особенностью ферромагнетиков является магнитострик- ция — способность изменять объем тела при намагничивании. К основным химическим элементам ферромагнетиков относятся Fe, Co и Ni. Наиболее распространенные ферромагнитные мине- минералы, содержащие окиси железа,— это магнетит Fe3O4, титаномаг- нетит Fe3O4-TiFe2O4, пирротин FeflS^i (п = 6—11) "и некоторые разновидности гематита Fe2O3. По величине магнитной восприимчивости горные породы под- подразделяются на четыре группы: очень сильномагнитные, сильно- 91
магнитные, среднемагнитные и слабомагнитные. Осадочные породу обладают слабой магнитностью. Их магнитные свойства опреде- определяются содержанием частиц ферромагнитных минералов, сильны^ темноцветных парамагнетиков и слабых парамагнитных и диамаг* нитных минералов. Для глинистых пород % = A0-f-14)-10~3, %д% песчаников — A4 —15)-10~3, а для гидрохимических и карбо- карбонатных пород — менее 6-10~3 А/м. Под действием магнитного попп Земли горные породы в период своего формирования способны на- намагничиваться и сохранять приобретенную (остаточную) намагни- намагниченность в последующие геологические эпохи. По величине и на- направлению остаточной намагниченности пород определяют магнит- ное поле, существовавшее в данной точке земной поверхности при образовании породы. На этом основывается палеомагнетизм — об- область знаний, занимающаяся изучением эволюции во времени гео- геомагнитного поля. Определение магнитных свойств пород В лабораторных условиях обычно определяют магнитную вос- восприимчивость % и остаточную намагниченность /. Для этого раз- разработан ряд способов, среди которых наиболее широкое примене- применение получили магнитометрические и индукционные. Первые осу- осуществляются с помощью астатических магнитометров и основаны на взаимодействии магнитов, одним из которых является намагни- намагниченный образец породы, а другим — астатическая магнитная си* стема. Определение к индукционным способом производят измери- измерителями магнитной восприимчивости (каппаметрами) ИМВ-2 шли ИМВ-3, а остаточной намагниченности / — измерителем остаточ- остаточной намагниченности ИОН-1. Тепловые свойства Основными тепловыми свойствами горных пород являются теп- теплопроводность Я или тепловое сопротивление породы I = 1Д» теплоемкость или удельная теплоемкость с и температуропровод- температуропроводность породы а. Коэффициент теплопроводности % характери- характеризует способность горных пород к передаче тепла и численно пока* зывает поток тепла в ваттах в единицу времени через породу се- сечением 1 м2, высотой 1 м при разности температур 1 К и выражается в Вт/(м-К). В промысловой геофизике обычно пользуются величи- величиной обратной теплопроводности — тепловым сопротивлением по* роды I (в м-К/Вт). Теплопроводность способствует переносу теплоты (энергии теп* лового движения микрочастиц) от более нагретых частей тела к ме- менее нагретым, что приводит к выравниванию температуры* При температуропроводности перенос энергии осуществляется нело» средственной передачей энергии от частиц (молекул, атомов, элек- электронов), обладающих большей энергией, частицам с Ш%& энергией. 92
В результате многочисленных исследований тепловых свойств 1Х пород установлено, что тепловое сопротивление понижается увеличением их плотности. В связи с этим изверженные и мета- ° оЛнческие породы имеют меньшее тепловое сопротивление, чем ^садочные песчано-глинистые и главным образом глинистые от- °ожения. С глубиной плотность горных пород закономерно возрас- возрастает что приводит к уменьшению в этом направлении теплового сопротивления и геотермического градиента (см. § 27). Тепловое сопротивление зависит от слоистости горных пород: в направлении, перпендикулярном к напластованию, тепловое со- сопротивление выше, чем в направлении напластования. Это явление известно под названием тепловой анизотропии. По- Понижение теплового сопротивления по напластованию связано с цир- циркуляцией вод в этом направлении и возникающим в результате дополнительным переносом тепла — конвекцией. Этим же объясняется повышение геотермического градиента под сводами складок, содержащих мощные пачки проницаемых пород, в кото- которых перемешивается вода. Свойство среды поглощать тепловую энергию при теплообмене оценивается удельной теплоемкостью (массовой см и объемной су)- Под удельной теплоемкостью понимают количество тепла в Дж, необходимое для нагрева 1 кг данного вещества на 1 К, и выражают в Дж/(кг-К). Изменение температуры различных по- пород при поглощении или отдаче ими тепла может происходить с различной скоростью. Эта скорость изменения температуры пород характеризуется комплексным параметром, называемым темпе- температуропроводностью (II 120) где 8„ — плотность пород. Величина а показывает изменение температуры единицы объема породы за единицу времени и выражается в м2/с. У многофазных пористых пород (водо-, нефте- и газонасыщен- газонасыщенных) передача тепловой энергии происходит не только путем тепло- теплопроводности, но и путем конвекции заполнителя порового про- пространства. Величина теплопроводности (теплового сопротивления) в многофазных осадочных породах определяется объемным соотно- соотношением фаз, их химико-минералогическим составом, агрегатным состоянием каждой фазы, удельной поверхностью и др. Тепловые свойства основных породообразующих минералов изменяются не- незначительно. Несколько повышенным тепловым'"сопротивлением и* пониженной теплопроводностью обладают глинистые минералы твердой фазы [7]. Данные о тепловых свойствах горных пород широко используются при термических исследованиях бурящихся и эксплуатационных скважин и решении задач, связанных с раз- разведкой и разработкой месторождений нефти и газа, 93
Определение тепловых свойств пород Определение основных тепловых свойств пород — ффищ^ ентов теплопроводности X, удельной (массовой сш объемной # i теплоемкости и температуропроводности а — может производиться способами стационарного и нестационарного тепловых режимов, : Нестационарные способы (способы нерегулярного тепловоз поля) характеризуются простотой измерительных установок, боль* шей производительностью, и поэтому при массовых исследования^; пород чаще применяются на практике. Среди большого числа фикаций способов нерегулярного теплового поля наиболее страненными являются зондовые, которые могут различаться щ форме зондов (плоская, сферическая) и продолжительности дейог*. l^JrS^ I Рис. 53 Схема установки для определения тепловых свойств пород / __ образец породы; 2 — нагреватель; 3 — горячий спай термопары; 4 — кронштеШ* 5 — холодный спай термопары, 6 — автоматический электронный потенциометр tit!» ЭПП-09; 7 — электромеханический секундомер типа ЭСПВ-53 ГЦ, 8 ~ электронный се- секундомер типа ЭС-2; 9 — магнитный пускатель; 10 — амперметр, 11 — вольтметр; ?$ «— автотрансформатор типа РНО-250-2; 13 — струбционный зажим, 14 — ключ магнитного пускателя вия (мгновенно, импульсно, постоянно) источника тепла [10]. одной из установок с плоским и мгновенно действующим источни- источником тепла постоянной мощности показана на рис. 53. При работе с этой установкой определяют количество тепла Q^ выделенное нагревателем за короткий промежуток времени t ш максимальную температуру Ттах, которую приобретает за врем imax часть составного образца породы толщиной х и сечением 5» расположенная между нагревателем и горячим спаем дифферен- дифференциальной термопары (см. рис. 53). Количество выделенного нагревателем тепла: AШЦ где / — время действия нагревателя (в с); (/ — напряжение на за- зажимах нагревателя (в^В); / — сила тока ^проходящего через на- нагреватель (в А). 94
95
96
По результатам этих определений при условии, что Т < Ттах, рассчитывают коэффициенты: F 1) теплопроводности (в Вт/м-К) 0,5 ^-+0,37-А—?]; A1.122) r J max L rmax *max 2) температуропроводности (в м2/с) + 0,5 —*— + 0,042^—— Y]l0-*; A1.123) 3) удельной массовой (в Дж/кг-К) и объемной (в Дж/м3-К) теплоемкостей В табл. 7 приведены основные физические свойства наиболее распространенных осадочных пород и породообразующих минера- минералов* Заказ № 789
ЧАСТЬ В ТОР А Я ГЕОФИЗИЧЕСКИЕ МЕТОДЫ ИССЛЕДОВАНИЯ РАЗРЕЗОВ СКВАЖИН Изучение геологического разреза и выявление полезных иско- ископаемых в скважинах проводятся электрическими, радиоактив- радиоактивными, акустическими, ядерно-магнитными, термометрическими, гео- геохимическими и другими видами геофизических исследований. Объек- Объектом геофизических исследований являются скважины, не обса- обсаженные колонной, бурящиеся, действующие, обсаженные колонной эксплуатационные, нагнетательные и остановленные. Основным объектом, подлежащим изучению каротажем в сква- скважинах, не обсаженных колонной, является геологический разрез, вскрываемый скважиной. По данным каротажа определяют после- последовательность и глубину залегания пород в разрезе, ихлитолого- петрографическую характеристику; наличие полезных ископаемых, качественное и количественное содержание последних. На результаты измерений геофизическими методами влияют породы, вскрытые скважиной, и промывочная жидкость, запол- заполняющая ее ствол, физические свойства которых различны. Из-за воздействия на стенки скважины промывочной жидкости в ряде случаев наблюдается отклонение фактического диаметра скважины от номинального (см. § 32). При бурении плотность глинистого раствора обычно поддер- поддерживается такой, чтобы гидростатическое давление его столба в сква- скважине было выше пластового. Под действием разности давлений фильтрат глинистого раствора проникает в проницаемые пласты, а его твердые частицы откладываются на стенке скважины, обра- образуя глинистую корку. Фильтрат раствора вытесняет нефть, гаа и воду, насыщающих поры породы. В результате физические свой- свойства пласта, и в первую очередь его удельное сопротивление вблизи стенок скважины, существенно отличаются от свойств неизменной части пласта. Для правильной интерпретации результатов геофизических ис- исследований скважин необходимо учитывать изменение диаметра скважины, а также физические свойства промывочной жидкости и зоны проникновения его фильтрата в пласт (рис. 54). Глубина проникновения его зависит от пористости и проницаемости пласта» проницаемости глинистой корки, разности гидростатического и пластового давлений, а также от продолжительности фильтрация» Установлено, что определяющими факторами являются пористость пласта и проницаемость глинистой корки. Так как проницаемость корки на 3—4 порядка ниже проницаемости коллекторов, коли- количество проникшего фильтрата раствора будет в основном зависеть 98
от проницаемости корки. Вместе с тем проникший в пласт фильтрат, заполняя объемы порового пространства, распространяется на большую глубину в низкопористом пласге. В связи с этим наиболь- наибольшая глубина проникновения наблюдается в проницаемых, но низ- низкопористых пластах: известняках, доломитах, сцементированных алевритах и песчаниках. Высокопористые породы обладают мень- меньшей глубиной проникновения фильтрата промывочной жидкости в пласт. В водоносном пласте зсна проникновения имеет следующее строение. Непосредственно у стенки скважины формируется про- промытая зона, в которой фильтрат промывочной жидкости почти пол- полностью вытесняет пластовую во- воду. Считается, что размеры этой зоны изменяются от 5—6 до 10— 13 см, а степень вытеснения пла- пластовой воды зависит от количества профильтровавшейся жидкости, параметров пористой среды _ и промывочной жидкости. Обычно удельное электрическое сопротив- сопротивление фильтрата рф больше удель- удельного сопротивления пластовой во- воды рв» в результате чего удельное сопротивление промытой зоны рпп превышает удельное сопротивле- сопротивление неизменной части пласта рвп. За промытой зоной следует переходная, сопротивление насы- насыщающей жидкости в которой изме- изменяется от Рф до рв за счет посте- постепенного смещения фильтрата гли- глинистого раствора с пластовой во- водой. Очевидно, что удельное со- сопротивление в переходной зоне из- изменяется от рпп до рвп. Промытая и переходная зона образуют зону проникновения, которую условно считают концентрическим слоем диаметром D, сопротивлением рзп и сопротивлением насыщающей жидкости pBD (см. рис. 54). В качестве величин рзп и D принимаются такие значения, влияние которых на результаты измерений эквивалентно влиянию фактической неоднородной зоны проникновения. В нефтеносном пласте процесс проникновения более сложен. В промытой зоне нефтяного пласта происходит замещение пласто- пластовой воды и нефти фильтратом глинистого раствора, но в тонких и тупиковых порах нефть частично сохраняется. Принято считать, что в промытой зоне содержится 15—25 % остаточной нефти. В гли- глинистых коллекторах, а также при большой вязкости нефти остаточ- 4* 99 Рис. 54 Схема строения прони- проницаемого пласта, вскрытого сква- скважиной / — глинистая корка; 2 — промытая зона; 3 — переходная и окаймляющая зоны; 4 — неизмененная часть пласта; D —¦ диаметр зоны проникновения; dnn —диаметр промытой части плас- пласта; dQ — диаметр скважины, hrK — толщина глинистой корки, h — мощ- мощность пласта
ная нефтенасыщенность достигает 30 % и более.В газоносных пла- пластах остаточная газонасыщенность всегда больше остаточной нефте- насыщенности даже для очень вязких нефтей. Остаточная газонасы- газонасыщенность обычно принимается равной 30 %. При удалении от сте- стенок скважины фильтрат бурового раствора в зоне проникновения смешивается все с большими порциями пластовой воды и нефтщ. На процесс проникновения фильтрата глинистого раствора в нефте- нефтегазоносный пласт, представленный гидрофильными породами, су- существенно влияют относительная проницаемость пород и началь- начальное распределение флюидов. В связи с более высокой фазовой проницаемостью нефти в нефте- нефтеносном пласте сравнительно с фазовой проницаемостью воды у внешней границы зоны проникновения может образоваться зона фронтального вытеснения жидкости, так называемая окаймляю- окаймляющая зона с повышенным содержанием воды, по минерализации близкой к пластовой. Сопротивление этой зоны рок значительно ниже рзп. Образование окаймляющей зоны возможно, если насы- насыщенность пластовой водой несколько превышает остаточную и в пласте имеется подвижная соленая вода. При значительной во* донасыщенности пласта окаймляющая зона не образуется. Наличие такой зоны в гидрофильных породах не может служить признаком нефтеносности пласта и лишь свидетельствует о том, что в пласте имеется некоторое количество подвижной воды и скважина при эксплуатации может давать нефть с водой. Окаймляющая зона со временем обычно исчезает. В соответствии с конструкцией скважины ее ствол на опреде- определенных интервалах по мере углубления и при достижении проект- проектной глубины закрепляется стальной колонной. В результате от- отдельные интервалы скважины могут быть перекрыты двумя и бо- более колоннами. В целях закрепления колонны и разобщения ме- между собой продуктивных и водоносных пластов пространство ме- между стенкой скважины и колонной заполняется цементным камнем. Исследование разрезов таких скважин осуществляется радиоактив- радиоактивным каротажем, термометрическими и акустическими измерениями» На показания диаграмм, полученных этими методами, помимо пород, перекрытых колонной, дополнительное влияние оказывают обсадные колонны (толщина их стенок), толщина цементного камня и степень его схватывания с колонной и породой, наличие затруб* ной циркуляции жидкости, зона проникновения в пласте, наличие в колонне воды, нефти, газа. На завершающем этапе бурения перед вводом скважины в экс- эксплуатацию или под нагнетание жидкости проводится контроль тех- технического состояния колонны для выявления дефектов (негерме* тичность колонны, некачественное цементирование). Контроль осу- осуществляется с помощью акустического каротажа (АКЦ), радиоак- радиоактивных, термических и других видов геофизических исследований. Путем прострела отверстий в колонне и цементе (перфорацией) осуществляется сообщение между пластом, намеченным к опробо- 100
ию ЙЛи разработке, и скважиной. Выбор интервала перфора- перфорации проводится по результатам геофизических исследований. Действующие скважины (эксплуатационные и нагнетательные) гйстематически контролируются в целях изучения эксплуатацион- эксплуатационных характеристик пласта с помощью определения отдающих и по- пшдающих интервалов (профилей притока в эксплуатационных приемистости в нагнетательных скважинах); выявления обвод- неиных слоев; установления затрубной циркуляции и других за- задач, возникающих в процессе разработки залежи. Исследование действующих нефтяных и газовых скважин осу- осуществляется малогабаритными геофизическими приборами, чаще радиоактивного каротажа (диаметром 25—50 мм). Эти приборы опускаются в скважину через колонну насосно-компрессорных труб в фонтанирующих скважинах или через серповидный зазор между штангами глубинного насоса и обсадной колонной в скважинах, оборудованных штанговыми глубинными насосами. Одновременно в действующих скважинах проводится систематическое исследова- исследование состава флюида в стволе скважины геофизическими методами. В остановленных или неработающих скважинах, из которых из- извлечено технологическое оборудование, создаются условия для гео- геофизических исследований, такие же, как и в скважинах на завер- завершающем этапе бурения (перед вводом скважины в действие). Глава III ЭЛЕКТРИЧЕСКИЙ КАРОТАЖ Электрический каротаж состоит в основном из двух модифика- щй: метода сопротивлений и метода самопроизвольно возникаю- возникающего электрического поля (естественных, собственных потенциа- потенциалов), Основными видами каротажа по методу сопротивления яв- являются каротаж нефокусированными (обычными) зондами, в том числе боковое каротажное зондирование (БКЗ), боковой и индук- индукционный каротаж, микрокаротаж. Сущность электрического каро- каротажа заключается в проведении измерений, показывающих изме- изменения вдоль скважины кажущегося удельного сопротивления (КС) пород и естественных потенциалов (ПС) для изучения геологиче- геологического разреза скважины. Результаты измерений изображаются в веде кривых изменения параметров КС и ПС вдоль ствола сква- скважины, § 10. КАРОТАЖ СОПРОТИВЛЕНИЯ НЕФОКУСИРОВАННЫМИ ЗОНДАМИ Под каротажем сопротивления нефокусированными зондами понимают электрический каротаж, основанный на измерении ка- кажущегося удельного сопротивления горных пород трехэлектрод- ными нефокусированными зондами. 101
Основы теории потенциала электрического поля Для определения удельного сопротивления горных пород в сква- скважине используется источник тока, создающий в окружающей среде электрическое поле. Допустим, что в неограниченную проводящую среду при помощи электродов А и В вводится ток, создающий в ней электрическое поле (рис. 55). Такое поле тождественно электриче* скому полю зарядов электродов Л и В, помещенных в непроводя- непроводящую среду. Разница заключается лишь в том, что в электрическом поле заряды неподвижны, а в проводящей среде они находятся в движении, непрерывно возобновляясь источником тока. -F { '4: i У А — Рис 55. Схема ввода тока в неограничен- неограниченную проводящую среду с помощью зазем- заземлений (электродов) А и В Рис. 56. Направление напряженности ля Е Электрическое поле характеризуется напряженностью Е> которая является вектором, имеющим величину и направление. За единицу напряженности электрического поля принимается вольт на метр (В/м), т. е. напряженность электрического поля, при котором между точками, находящимися на расстоянии 1 м, вдоль линии напряженности поля (отражающей ее направление) соз- создается разность потенциалов 1 В. Под линией напряженности элек- электрического поля, называемой чаще силовой линией, под- подразумевают такую линию, в каждой точке которой вектор напря- напряженности направлен по касательной к ней (рис. 56). Силовые ли- линии соответствуют путям, вдоль которых должен был бы двигаться положительный заряд. При помощи этих линий можно наглядно изобразить силовое поле; при этом густота его линий выбирается пропорционально напряженности. Работа, совершаемая силами электрического поля при переме- перемещении единичного положительного заряда из некоторой точки в бесконечно удаленную, численно равна электрическому потен- потенциалу данной точки (с обратным знаком). Потенциал есть величина скалярная и в каждой точке поля имеет вполне определенное зна- значение, поэтому может служить характеристикой поля наравне с на- напряженностью Е. За единицу электрического потенциала пршш* 102
мается вольт (В)—разность потенциалов между двумя точками при постоянном токе силой 1 А, в котором затрачивается мощность 1 Вт. Потенциал электрического поля представляет собой функцию, которая изменяется от точки к точке и убывает в направлении хода силовой линии. В каждом реальном случае можно выделить сово- совокупность точек, потенциалы которых одинаковы. Геометрическое б Рис. 57. Эквипотенциальные поверхности (пунктирные линии) и силовые линии (сплошные) а — точечный заряд; б — два разноименных численно равных заря- заряда; в — два точечных одноименных заряда место точек постоянного потенциала называют эквипотенциальной поверхностью. Если путь перемещения заряда замкнут по эквипо- эквипотенциальной поверхности, то работа электрических сил равна нулю. Известно, что потенциал точечного заряда е в точке, отстоя- отстоящей от него на расстоянии /*, U = -j- ^Er. (III.l) Следовательно, эквипотенциальная поверхность с постоянным значением г есть сфера с центром в точечном заряде (рис. 57, а). 103
Между напряженностью поля Е и потенциалом U имеется деленная связь. Существование такой связи следует из того, работу электрических сил можно выразить через напряженность и разность потенциалов точек поля. Рассмотрим однородное электрическое поле напряженностью Ё (рис. 58). Расстояние Аг между эквипотенциальными поверхностямЕ Ux и U 2 бесконечно мало, следовательно, на всем расстоянии ме- между ними можно считать напряженность поля постоянной. Работа перемещения единичного заряда е на пути дг равна Едг. Эта же ра- работа может быть выражена через разность потенциалов начала U и конца 1/г с обратным знаком: t/a—t/i = -Едг. пит ^-Е Рис. 58. Однородное электростатическое поле напряженностью Е Аг — расстояние между эквипотенциальными поверхностями Ut и U2 Рис. 59. Элементарный объем среды с удельным сопротивлением р, длиной dr и сечением dS В неоднородном поле силовые линии не будут прямолинейными* а эквипотенциальные поверхности будут иметь сложную форму. Однако для бесконечно малых участков пространства можно дре* небречь кривизной силовых линий и эквипотенциальных поверх- поверхностей и на основании предыдущих рассуждений записать: -Едп дг (Ш.З} Знак минус указывает, что напряженность Е направлена в ту сторону, в которую действует сила на положительный заряд, т, е> в сторону убывания потенциала. Величина dUldr, характеризующая быстроту изменения потей* циала при перемещении в направлении, перпендикулярном к экви- эквипотенциальным поверхностям в сторону его увеличения, называется градиентом потенциала и обозначается grad U, Как видно из формулы (Ш.З), напряженность поля — это градиент потенциала с обратным знаком, т. е. ?= -gradt/. "(IIU) 104
Сила тока / представляет собой физическую величину, изме- изменяемую количеством электричества, перенесенного через данную площадку в единицу времени, независимо от того, в каком направ- направлении и под каким углом к площадке движутся частицы, несущие заряды. Для учета направления переноса зарядов вводится в рассмот- рассмотрение вектор плотности тока /, который направлен в сторону движения положительных зарядов, т. е. в направлении вектора напряженности Е. Под плотностью тока пони- понимается количество электричества, протекающее в единицу времени Ч0рез единичную площадку, перпендикулярную к направлению тока. Если ток / равномерно распределен по площади S, перпен- перпендикулярной к его направлению, то величина плотности тока / = //S. (III.5) В общем случае плотность тока определяется отношением силы тока dl, протекающей через перпендикулярный к направлению элемент сечения среды, к площади dS этого элемента (III.6) Распределение электрического поля в пространстве удовлетво- удовлетворяет двум основным законам: Ома и Кирхгофа, выраженным в диф- дифференциальной форме. Для пояснения закона Ома выделим эле- элементарный объем (рис. 59) среды с удельным сопротивлением р, длиной dr и сечением dS; через сечение dS и перпендикулярно к нему проходит ток dl, образуя на концах падение потенциала dU. По формуле (II. 1) сопротивление элементарного объема будет CLi\ ¦— Р , dS а падение потенциала на его концах — du = dlp — dS или dl I dU dS p dr Пользуясь уравнениями (III.2) и (III.6), получаем / = E/p. A11.7) Закон Ома в дифференциальной форме выражается так: плот- плотность тока в каждой точке проводника равняется напряженности электрического поля в этой точке, деленной на удельное сопротив- сопротивление вещества. Физическая сущность первого закона Кирхгофа в дифферен- дифференциальной форме заключается в том, что если какой-либо элемент объема не содержит источников, то сила тока, втекающего в этот 105
объем, равна силе тока, вытекающего из него. Этим выражается непрерывность потока токовых линий через любую замкнутую по- поверхность, не содержащую дополнительных источников тока. Если считать, что входящие и выходящие из данного объема токи имею? разные знаки, то алгебраическая сумма их равна нулю, т. е. эти токи по величине равны. Невыполнение этого условия привело бы к накоплению электрических зарядов в некоторых точках, что ис- исключается. Закон Кирхгофа в дифференциальной форме записывается ш виде следующего уравнения: div / = 0. (IILB) Рассмотрим электрическое поле в однородной изотропной среде. Предположим, что имеется такая среда с удельным сопротивле- сопротивлением р. Введем в нее электрод Л, из ко* торого вытекает ток силой /. Допустим» что размеры электрода Л малы и его можно рассматривать как точечный, а второй электрод удален (теоретически в бесконечность) и не оказывает влия- влияния на электрическое поле вблизи электрода Л. При этих условиях линии тока, исходящие из точки Л, будут прямыми, а эквипотенциальные поверх- поверхности — концентрическими сферами с центром в точке Л. В пересечении с плоскостью чертежа эти сферы дают окружности с центром в точке А (рис. 60). Определим потенциал в точке М, расположенный на расстоянии г от ис- источника тока Л. Плотность тока / в точке М согласно формуле (III.5) / = //4jtra, (II 1.9) так как, исходя из формул (III.5) и (III.9), ди Рис. 60. Электрическое по- поле в однородной изотроп- изотропной среде ТО Е = — dU дг dr 4лг2 Зная это, потенциал в точке М будет p/dr р/ 4зх/2 4лг Так как потенциал в бесконечности равен нулю, т. е. при г = = оо, U = 0, то постоянная интегрирования С = 0. Тогда выра- 106
жение потенциала в точке М, созданного точечным электродом Л, через который протекает ток /, в однородной изотропной среде с удельным сопротивлением р примет вид р/ (ШЛО) Легко видеть, что если поменять местами точки А и М, т. е. источник тока поместить в точку М и определять потенциал в точке ЛУ то его величина выражается также с помощью уравнения (ШЛО). Это положение справедливо и для неоднородной среды, оно нахо- находит важное практическое применение в электрическом каротаже и известно под названием принципа взаимности. Измерение кажущегося удельного сопротивления обычными зондами Для замера сопротивления пород, пересеченных скважиной, используют четырехэлектродную установку AMNB. Три электрода этой установки (А, М, N или М, Л, В), присоединенные к концам кабеля и опускаемые в скважину, пред- представляют каротажный зонд. Четвертый электрод В или N (заземление) уста- устанавливают на поверхности вблизи устья скважины (рис. 61). Через электроды Л и В, называемые токсвыми, пропускают ток /, создаю- создающий электрическое поле в породе; при помощи измерительных электродов М и N регистрируют разность потенциа- потенциалов AU между двумя точками этого электрического поля. Согласно прин- Рнс. 61. Схема измерения кажущегося удель- удельного сопротивления А, В и М, N — токовые н измерительные электроды; П — измерительный прибор, К — трехжильнын ка- кабель; Е — источник тока, R — сопротивление для установки силы тока в цепи питания, тЛ — милли- миллиамперметр ципу взаимности при каротаже сопротивления допускается взаимная замена токовых и измерительных электродов. Регистри- Регистрируемая величина при этом является одной и той же. Потенциал некоторой точки М, созданный точечным электродом Л, через который протекает ток /, в однородной изотропной среде с удельным сопротивлением р согласно формуле (ШЛО) будет 4яг м 107
Подставляя вместо гм величины AM или AN, получаем потен- потенциал в точке М: АлАМ и потенциал в точке N: - р/ N An AN Если считать электроды'зонда точечными, то разность потенциа- потенциалов между его измерительными электродами М и N будет р/ / 1 1 \ р/ (AN —AM) Лп .. -г, р/ / 1 1 \ 4л \ AM ANj An AM-AN где AN—AM - Тогда An AM AN ; По формуле (III. 12) можно вычислить удельное сопротивление однородной среды AU AnAM AN MN (Ш.13) Все величины, входящие в правую часть формулы (III. 13), можно измерить и таким образом определить величину удельного сопротивления среды, что и является целью электрического каро- каротажа по методу сопротивлений. При каротаже разность потенциа- потенциалов выражается в тысячных долях вольта — милливольтах (мВ), сила тока в тысячных долях ампера — миллиамперах (мА), а рас- расстояния MN, AM и AN в метрах (м), при этом удельное сопротив- сопротивление будет выражено в ом-метрах (Ом-м). Приведем формулу (III. 13) к виду, в котором она обычно при- применяется в практике электрического каротажа. Для этого, полагая, что 4nAMAN =K, (IIU4) MN получим р = /С-^-, A11.15) где К — коэффициент зонда — постоянный множитель, зависящий от расстояний AM, AN и взаимного расположения электродов; AU/I представляет собой сопротивление части среды, заключен- заключенной между двумя эквипотенциальными поверхностями, проходя- проходящими через точки М и N. 108
Выражение (III. 15) справедливо для вычисления истинного удельного сопротивления р изотропной и однородной среды. При этом условии значение удельного сопротивления должно оставаться постоянным при любых расстояниях AM и AN. При каротаже всегда имеют дело с неоднородной средой, состоя- состоящей из пластов различного удельного сопротивления и промывоч- промывочной жидкости, заполняющей скважину. В этих условиях получен- полученный по формуле (II 1.15) результат является условным (фиктивным) и назван кажущимся удельным сопротивле- й и е м (КС или рк). Кажущееся удельное сопротивление среды рк численно равно истинному удельному сопротивлению такой одно- однородной среды, показания в которой при заданных зонде и силе пи- питающего тока равны показаниям в данной неоднородной среде. Кажущееся удельное сопротивление зависит от многих факторов: удельного сопротивления и мощности пластов, против которых находится каротажный зонд, диаметра скважины и удельного со- сопротивления заполняющей ее промывочной жидкости, глубины проникновения фильтрата промывочной жидкости в пласт, харак- характера взаимного залегания и сопротивлений смежных пластов, типа и размера зонда, которым проводят измерения. Кажущееся удельное сопротивление, замеренное в однородной анизотропной среде (см. § 5), не зависит от типа и размера зонда; оно пропорционально среднему удельному сопротивлению анизот- анизотропной среды ртСр и зависит от коэффициента анизотропии %а и угла (~ ом, составленного скважиной и направлением падения пластов.1 Значение кажущегося сопротивления при этом определяется формулой рк = pmcp . (III.16) Из этой формулы следует, что кажущееся сопротивление ани- анизотропной среды, замеренное по перпендикуляру к напластованию р , равно истинному удельному сопротивлению этой среды, за- замеренному в плоскости напластования р((: (Ш.17) Если же замер произведен в плоскости напластования, то ка- кажущееся сопротивление рк [( равно среднему удельному сопротивле- сопротивлению этой среды ртср или ркц =Яарц. (III. 18) 1 Угол а—видимый угол падения пород (а—угол, дополняющий до пря- прямого между осью скважины и плоскостью напластования). При вертикаль- ной скважине а равен истинному углу падения пород. 109
Так как коэффициент анизотропии пород больше единицы (см, § 5), то кажущееся сопротивление в направлении, перпендикуляр- ном к напластованию, меньше кажущегося сопротивления по на- напластованию (Рк L < Рк л) ¦ Между тем известно, что истинное удельное сопротивление ани- анизотропных пород в направлении, перпендикулярном к напластова- напластованию, больше истинного удельного сопротивления по напластова- напластованию (р±>Р„). Такое несоответствие между значениями кажущихся сопротив- сопротивлений и истинных называется парадоксом анизотро- a I S Обращен- Лоследова - ный тельный. л: Обращен- Последова- ный тельный о >-г-ХВ МАА- д %-у-Ам N СИ Рис. 62. Обозначение зондов / — градиент-зонд а — кро- кровельный, б — подошвенный; // — потенциал-зонд; циф- цифры в кружках. 1 — однопо* люсный зонд (зонд прямого питания), 2—двухполюсный (взаимного питания) зонд; 1 — токовые электроды (Л, В); 2—измерительные элек- электроды (МАО; 3 — точка за- писи кажущегося сопротив- сопротивления, 4 — электроды для замера ПС п и и и вызвано увеличением плотности тока вдоль напластования в направлении повышения электропроводности анизотропных (слои- (слоистых) пород. Результаты измерения кажущегося сопротивления пород в скважине изображаются в виде кривой изменения КС вдоль ствола скважины (см. рис. 66, 70). Для измерения КС пород при каротаже применяют зонды (измерительные установки, содержащие три электрода: Л, М, N или М, Л, В) различных типов и размеров, из которых выделяют два основных типа: градиент-зонды и потенциал-зонды (рис. 62). Градиент-зондами называют зонды, у которых рас- расстояние между парными электродами М и N или Л и В мало по сравнению с расстоянием непарных электродов Л и М или М и Л. Замер кажущихся сопротивлений этим зондом сводится к измере- измерению градиента-потенциала электрического поля электрода Л. Дейст- Действительно, при бесконечно малом MN формулу (III. 13) можно за- 110
писать так: (AU/MN) E , (Ш. 19) гд6 ? — составляющая напряженности электрического поля по оси z в точке О или градиент-потенциала с обратным знаком. Размером градиент-зонда является величина АО; где О — точка записи кривой КС, находящаяся в середине между парными элек- электродами М, N (А, В). От величины АО зависит глубина исследова- исследования, которая тем больше, чем больше размер зонда. Градиент- зонд, У которого сближенные парные электроды расположены под непарным, называют подошвенным. При расположении сближен- сближенных парных электродов над непарным градиент-зонд называют кро- кровельным. Градиент-зонд, у которого расстояние между парными электродами бесконечно мало, называется идеальным. Потенциал-зонды — это зонды, у которых расстояние AM мало по сравнению с расстоянием между парными электродами М, N (А, В). Расстояние AM является размером потенциал-зонда. Замер кажущегося сопротивления относят к середине AM, Кажу- Кажущееся сопротивление в этом случае АЕ/ _ АМ-AN АЕ/ ^^ / AM-AN\ Um — Un Рк~ / ~ ^ MN I ~ V MN ) I (III.20) Потенциал-зонд с электродом N, удаленным в бесконечность, называется идеальным. Для такого зонда AN ~ оо, МА^^оо, Vh = 0, отсюда ^ (III.21) Кажущееся сопротивление при использовании потенциал-зонда определяется потенциалом электрического поля в точке М. Поэтому зонды такого типа и называют потенциал-зондами. Для условной оценки глубины исследования зондом применяют термин радиус исследования зонда — радиус сферы в однородной среде неограниченной мощности, оказывающей на показания зонда такое же влияние, как и та часть сферы, ко- которая расположена за ее пределами. Исходя из этого считают, что радиус исследования градиент-зондом приблизительно совпадает с его размером АО, а потенциал-зондом соответствует его удвоен- удвоенному размеру, т. е. 2АМ. Следовательно, при одинаковом размере зондов радиус исследования потенциал-зонда примерно в 2 раза превышает радиус исследования градиент-зонда. Зонд с одним питающим электродом и двумя измерительными называется однополюсным (или зондом прямого питания), а с двумя питающими электродами и одним измерительным — двухполюсным (или взаимного питания). Коэффициент К при двухполюсном зонде вычисляют по формуле (III.22) 111
Рис. 63. Кривые сопротивления для однородного пласта с большим (а, б) и малым (в, г) сопротивлениями а, в — подошвенный градиент-зонд, б, г — потеидиал-зонд
Зонды записывают по обозначениям электродов в порядке их расположения в скважине сверху вниз, проставляя между ними расстояния в метрах. Так, например, градиент-зонд двухполюсный, подошвенный, у которого верхний электрод является измеритель- измерительным, а на расстояниях 2,5 и 2,75 м ниже его расположены соответст- соответственно первый токовый электрод А и второй В, будет обозначаться М2, 5ЛО, 25В. Кривые кажущегося удельного сопротивления против пластов ограниченной мощности Величина кажущегося удельного сопротивления, определяющая форму кривой КС, зависит от мощности пласта, типа и размера зонда? его положения относительно границ пласта. На рис. 63 при- приведены кривые, полученные в результате экспериментальных и тео- теоретических исследований для обычных зондов против однородных пластов ограниченной мощности и различного удельного сопротив- сопротивления. Условно принято считать пласт мощным, если его размер превышает размеры зонда, тонким, если его мощность меньше или равна его размерам. Если удельное сопротивление пласта соот- соответственно больше или меньше удельного сопротивления вмещаю- щей~среды, то пласт квалифицируют как пласт высокого или низ- низкого сопротивления. Градиент-зонт. Пласт высокого сопротивле- п и я* На кривой КС такой пласт отмечается асимметричным макси- максимумом. При замерах подошвенным градиент-зондом кровля пласта соответствует минимальному сопротивлению, а подошва — макси- максимальному. В действительности для реального зонда граница по- подошвы пласта фиксируется ниже максимума на половину расстоя- расстояния между сближенными электродами. Тонкому пласту соответст- соответствует максимум со слабо выраженной асимметрией. Кровля его на- находится против точки наиболее крутого подъема кривой, а по- подошва — несколько ниже максимума. Ниже подошвы пласта на длину зонда наблюдается повышение сопротивления, вызванное экранным максимумом (см. рис. 63). Пласт низкого сопротивления. Мощный пласт фикси- фиксируется"^ кривой сопротивления асимметричным минимумом. При за- замерах подошвенным градиент-зондом кровля пласта приблизи- приблизительно отмечается максимумом, а точнее — ниже него на половину расстояния между сближенными электродами, подошва — мини- минимумом. Для тонких пластов подошва на кривой КС фиксируется но переходу кривой сопротивления от пониженных значений к мак- максимальным (см. рис. 63, в). При измерениях кровельным градиент-зондом кривые сопротив- сопротивления являются зеркальным отражением кривых, полученных по- подошвенным градиент-зондом. Определение границ пласта кровель- кровельным градиент-зондом производится по тем же правилам, что и в случае" подошвенного, но с учетом обратного хода кривой. 113
Потенциал-зонд. Пласт высокого сопротив^и н и я. Пласт мощный отмечается на кривой КС максимумом, с\ш метричным относительно середины пласта. Его границы проводите» симметрично относительно максимума, кровля — на половвд длины зонда выше точки перехода от плавного к более т-— подъему кривой, а подошва — на ту же величину ниже этой Тонкий пласт высокого сопротивления фиксируется снижевдш сопротивления: некоторое повышение последнего наблюдается вщщ кровли и ниже подошвы пласта на расстояниях, равных половнщ длины зонда из-за экранных явлений (см. рис. 63, б). Пласт низкого сопротивления. Такой пласт на кривой кажущегося сопротивления отмечается минимумом симметричным относительно середины пласта. Его границы про- проводятся по точкам перехода от крутого спада к плавному понижен- пониженному участку кривой с учетом того, что эти точки смещены относи* тельно кровли и подошвы на половину длины зонда. Таким образом ширина минимума превышает мощность пласта на длину зонда! Выделение границ тонкого пласта малого сопротивления в этом случае затруднительно (см. рис. 63, г). При чередовании пластов, имеющих различные сопротивления обычное распределение плотности тока в скважине нарушается, происходит перераспределение силовых линий тока и возникшее явления экранирования, которые оказывают влияние на величины кажущихся сопротивлений и должны учитываться при интерпрета- интерпретации кривых кажущихся сопротивлений. На измерения градиент-зондом значительное влияние оказы- оказывает соседний пласт высокого сопротивления, расположенный ш стороны удаленного электрода. Если расстояние между серединами соседних пластов больше длины зонда, то происходит повышенщ кажущихся сопротивлений, а если меньше — понижение по срав- сравнению с теми, которые наблюдались бы в случае одиночного пласта (рис. 64). Против пачки чередующихся пластов большого и малого сопро- сопротивлений форма кривой зависит от числа составляющих пачку пла- пластов, их мощности и удельного сопротивления, а также от типа н длины зонда. Форма кривых кажущихся сопротивлений при чере- чередовании пластов достаточно полно изучена на модельных установ- установках [8]. Кажущееся удельное сопротивление различно против разшж точек пласта. Для определения истинного удельного сопротивле- сопротивления необходимо выбрать наиболее характерные (существенные) значения КС. Существенными значениями кажущегося сопротивле- сопротивления принято считать среднее рКСр> максимальное рк тах пли ми- минимальное рк min и оптимальное рк опт (рис. 65). Среднее значение КС соответствует отношению площади, огр&* ничейной нулевой линией диаграммы и кривой КС против пласта* к его мощности. На практике визуально проводится линия, парал- параллельная нулевой и отсекающая прямоугольник с основанием у ву* И4
oft линии, равный мощности пласта. Если площадь полученного 1$«\юугольника равновелика искомой, то высота его соответствует Среднему значению КС. Максимальное и минимальное значения КС (экстремальные со- тпвлеяия) отсчитывают для пластов, удельные сопротивления ушрых соответственно больше или меньше, чем у вмещающих пород. Рис. 64. Кривые сопротивления для двух пла- пластов, мощность которых h меньше длины зонда, записанные подошвенным градиент-зондом (рп— 10 рс и рвм = Рс) Мощность прослоя малого сопротивления соответственно равна для а—•< — h, i,5 h\ Ah (занижающее экраниро- экранирование), г — 8 h (завышающее экранирование), AM— 7,5 h, MN = h Рис. 65. Пример отсчета среднего, рк ср макси- максимального рк тах и оптимального рк опт сопротив- сопротивлений Показан пласт высокого сопротивления, запись велась кровельным градиент-зондом ВО, 25А2, 5М О 20 40 60 ПЛпСГПП По кривой сопротивления, полученной потенциал-зондом, мак- максимальное и минимальное значения КС отсчитывают против сред- средней части пласта. По кривым КС, полученным кровельным и по- подошвенным градиент-зондами, максимальное значение сопротив- сопротивления отсчитывают соответственно в кровле пласта и его подошве, а минимальное — у границы пласта, расположенной со стороны удаленного электрода. Оптимальное значение КС наиболее близко к истинному сопро- сопротивлению пласта. Оно соответствует величине рк в точке, располо- расположенной выше или ниже середины пласта приблизительно на полови- половину длины зонда при использовании соответственно кровельного или подошвенного градиент-зонда. Величина кажущегося удельного сопротивления пласта конечной мощности зависит от его сопротив- сопротивления, типа зонда, соотношения длины зонда и мощности пласта, а также расположения зонда относительно границ пласта. Для И^аста высокого сопротивления наибольший интерес представляют 115
средние и максимальные величины рк, измеренные градиентной! дом, и максимальные — потенциал-зондом. «* Фактические кривые сопротивления, щ^ писанные в скважине, имеют более сложную форму, чем расчетный или полученные на моделях. Сложность формы фактической щш вой обусловлена неоднородностью пласта и вмещающих порою изменением диаметров скважины и зоны проникновения фильтрат^ промывочной жидкости в пласт, углом между осью скважины ш плоскостью напластования. Для получения представления о разрезе, вскрываемом скважи- скважиной (выделение пластов различного удельного сопротивления и оп- определение их границ), необходимо выделить на фактической ¦г г 130 2140 PI50 2160 _г~ X X Х- ___•*_• КС В 2,5 А 0,5 М 50 ЮО 150 Ом м — 10 20 30 Ом м i 1 i A4M0,5N --250 500 750 0мм 50 W0 150 Ot/ м 10 20 30 Ом м i i i N 0,5 М4А 2$Q 500 750 Омм 50 /00 /50 Ом м 10 20 30 Ом м i i i \ ПС 25 иВ ^.i. Рис. 66. Фактические кривые сопротивления для двух пластов высокого противления и мощностью, большей и меньшей длины зонда 1 — глина; 2 — песчаник вой сопротивления основные элементы, известные по расчетным кривым. Это можно сделать в том случае, если известны основные причины, вызывающие осложнение кривой КС, и особенности ее изменения под влиянием различных факторов. На рис. 66 приведены кривые кажущихся сопротивлений КС, зарегистрированные в скважине при помощи потенциал-зонда В2, 5А0, 5М, подошвенного А4М0, 5N и кровельного N0,5 Шк градиент-зондов. Выделение границ пластов большого и малого сопротивлений выполнено в соответствии с изложенными выше по- положениями. Кажущееся удельное сопротивление пласта неограниченной мощности (боковое каротажное зондирование) Результаты расчета кажущегося удельного сопротивления пласта неограниченной мощности, выполненного Л. М. Альпиным, представлены в виде кривых, выражающих зависимость рк от раз- 116
ичных определяющих его параметров: а) для непроницаемого ласта — от удельных сопротивлений пласта рп и промывочной жидкости рс диаметра скважины dc и длины зонда L3; б) для проницаемого пласта при наличии зоны проникновения, кроме перечисленных параметров,— от удельного сопротивления зоны проникновения рзп и ее диаметра D. Эти кривые называются кри- кривыми бокового каротажного зондирования (БКЗ). Такие кривые, сгруппированные по определенному признаку (двухслойные, трех- трехслойные) и выражающие зависимость рк/рс от LJdc для пласта Рис. 67. Палетка БКЗ-1А для градиент-зондов при рп >• рс (по Л. М. Альпину) Рис. 68. Палетка БКЗ-1Б для градиент-зондов при рп > Рс (по Л. М. Альпину) Л'А. Ц1 0,2 0,5 5 10 Z0 50 100 200 Ж AQ/dz 0,005 0,005 0J 0,2 Iff ?UAOjdz неограниченной мощности, называются палетками БКЗ. Разли- Различают кривые БКЗ двух основных типов — двухслойные и трехслой- трехслойные. Двухслойные кривые БКЗ. рассчитаньГдля усло- условий, когда проникновение промывочной жидкости в пласт отсутст- отсутствует. При этом возможны следующие случаи: а) сопротивление промывочной жидкости, заполняющей скважину, меньше сопро- сопротивления пласта (рс<рп); б) сопротивление жидкости больше сопротивления пласта (рс>рп). Двухслойные расчетные кривые БКЗ сгруппированы в палетки, обозначаемые БКЗ-1А (при рп>рс) и БКЗ-1Б (при рл<рс) [рис. 67, 68]. Как видно, кривые палеток БКЗ-1 в своей правой части асимптотически приближаются к значениям удельного со- сопротивления пласта. Изображенная на палетках кривая А характе- характеризует геометрическое место точек пересечения кривых БКЗ с их правыми асимптотами, кривая В — геометрическое место точек (максимумов и минимумов) кривых. Двухслойные кривые БКЗ обоз- 117