Full Text

Введение Выделение проницаемых и непроницаемых пропластков в толще пород является одной из основных задач нефтегазопромысловой геологии. Разделение толщи пород служит основой для построения геологической модели, подсчета запасов и дальнейшего планирования разработки месторождения. Для расчленения продуктивного разреза скважины необходимо выделить слои различного литологического состава, установить последовательность их залегания и в конечном итоге выделить коллекторы и непроницаемые разделы между ними. Эти задачи решаются с помощью комплексного метода изучения разрезов. Основными в этом комплексе являются геофизические методы, применяемые в обязательном порядке в скважинах всех категорий (поисковые, разведочные, эксплуатационные и т.д.). Геофизические данные связываются с имеющимися геологическими данными описания пород (керна, шлама), с результатами испытания скважин на приток и с результатами гидродинамических исследований. Качество разделения зависит от степени изученности геологического разреза, уровня теоретической разработки геофизических методов исследования скважин и общей геофизической характеристики района. Интерпретация геофизических кривых наиболее достоверна в комплексе с геологическими исследованиями, при этом следует понимать, что керн в некоторых случаях не дает полного представления о положении границ раздела пластов, что связано с низким процентом его выноса и сложностями с привязкой кернового материала по глубине [1]. Объект исследования За объект исследования принята поисково-оценочная скважина, пробуренная в 2014 г. на Северо-Ефремовском куполе Софьинского месторождения. В административном отношении Софьинское месторождение расположено на территории Уинского, Чернушинского и Октябрьского районов Пермского края, в тектоническом приурочено к Таныпскому позднедевонскому атоллу, расположенному в северо-восточной части Башкирского свода [2]. В скважине 119 проведен комплекс геофизических исследований, включающих в себя: стандартный каротаж (А2М0.5N), боковое каротажное зондирование, микрозондирование, боковой, микробоковой, индукционный, акустический, радиоактивный каротажи и кавернометрию. Также в скважине с отбором керна пройдено 172,4 м пород, вынос керна составил 170,1 м (98,7 %). Определены коэффициенты пористости, проницаемости, нефтенасыщенности и объемной плотности керна. Выделены продуктивные пласты в верейском, верхневизейском и верхнедевонско-турнейском нефтегазоносных комплексах. Для исследования принимались результаты интерпретации ГИС и исследования кернового материала верхневизейского нефтегазоносного комплекса [3, 4]. Важно, что определение характеристик пород проводилось не только в коллекторах, выделенных по геофизическим данным, но и в неколлекторах. Объем полученных данных и их разнообразность позволили составить базу данных с представительной выборкой. Анализ выделения пород-коллекторов Разделение на проницаемую и непроницаемую части проводилось по комплексу промыслово-геофизических исследований по общепринятым методикам [5], также учитывались следующие качественные признаки кривых ГИС [6-11]: - низкие и средние показания относительно вмещающих пород на кривых радиоактивного каротажа. - повышенные значения сопротивлений относительно вмещающих плотных и глинистых пород. Используем критерий Стьюдента t [12-14] для анализа средних значений классов проницаемых и непроницаемых пород, полученных по результатам исследования кернового материала и по геофизическим методам (табл. 1, 2). Из геофизических данных статистически значимыми являются все параметры, кроме бокового каротажа и значений, полученных с помощью микропотенциал-зонда. Значения параметров значительно различаются: от 0,7-0,8 по данным микроградиент-зонда до 781,9-849,2 по данным акустики. По керновым данным статистически значимыми являются все параметры, кроме минералогической плотности. Значения параметров различаются (2,1-2,3 г/см3) по плотности образца до 68,7-904,1 мД по проницаемости породы. Такой разброс средних значений различных параметров не позволяет проводить их количественное сравнение. Для проверки показателей, наиболее сильно влияющих на дифференциацию пород, построены диаграммы рассеяния для данных, полученных в результате исследований кернового материала и по геофизическим параметрам (рис. 1, а, б). Наибольшее полное разделение на проницаемые и непроницаемые породы по керновым данным наблюдается при сопоставлении параметров коэффициента пористости и остаточной водонасыщенности (см. рис. 1, а). Отсюда видно, что породы, относящиеся к коллекторам, обладают высокой пористостью и низкой остаточной водонасыщенностью, и наоборот, породы-покрышки обладают низкой пористостью и высокой остаточной водонасыщенностью. Уравнение регрессии для проницаемых пород имеет вид Kв = 214,9 - 19,8Kп + 0,5Kп2, для непроницаемых пород Кв = 63,2 + 3Kп - 0,3Кп2. Анализ значений коэффициентов Z и свободных членов показывает, что они значительно различаются. Таблица 1 Значения t-критерия Стьюдента по данным ГИС Метод Среднее Значение критерия Стью-дента Число степеней свободы Уровень значи-мости Количество наблюдений Стандартное отклонение Относи-тельная дисперсия Дис-персия непрони-цаемый прони-цаемый непрони-цаемый прони-цаемый непрони-цаемый прони-цаемый BK 7,52 7,57 -0,04 334 0,97 147 189 15,68 5,35 8,59 0,00 BMK 4,59 11,19 -5,41 334 0,00 147 189 6,85 13,49 3,88 0,00 DTP 239,99 272,59 -10,66 334 0,00 147 189 9,43 36,12 14,67 0,00 GK 4,04 13,94 -32,82 334 0,00 147 189 1,63 3,36 4,25 0,00 IK 606,69 194,00 13,40 334 0,00 147 189 385,82 153,80 6,29 0,00 MGZ 0,84 0,71 2,61 334 0,01 147 189 0,16 0,58 12,71 0,00 MPZ 0,81 0,79 0,58 334 0,56 147 189 0,24 0,51 4,64 0,00 NKTB 4,75 2,58 17,22 334 0,00 147 189 0,83 1,33 2,55 0,00 NKTS 1,77 1,16 18,26 334 0,00 147 189 0,16 0,38 5,41 0,00 TP1 662,26 712,83 -9,86 334 0,00 147 189 33,61 54,66 2,65 0,00 TP2 781,98 849,22 -12,26 334 0,00 147 189 32,30 60,04 3,46 0,00 dGK 0,14 0,74 -32,82 334 0,00 147 189 0,10 0,20 4,26 0,00 W 21,04 32,64 -11,02 334 0,00 147 189 3,02 12,47 17,10 0,00 DS 0,22 0,24 -7,85 334 0,00 147 189 0,00 0,04 1961,62 0,00 Кгл 8,88 46,44 -32,82 334 0,00 147 189 6,18 12,76 4,26 0,00 КпАК 20,48 28,34 -10,66 334 0,00 147 189 2,27 8,70 14,67 0,00 Примечание: BK - боковой каротаж, Ом·м; BMK - зонд бокового микрокаратожа, Ом·м; DTP - интервальное время пробега Р-волны по базе широкополостного акустического каротажа, мкс/м; GK - интенсивность гамма-излучения по гамма-каротажу, мкР/ч; IK - индукционный каротаж, проводимость, мСм/м; MGZ - градиент-микрозонд, Ом·м; MPZ - микропотенциал-зонд, Ом·м; NKTB - нейтрон-нейтронный каротаж по тепловым нейтронам (ННКт), большой зонд, усл. ед.; NKTS - нейтрон-нейтронный каротаж по тепловым нейтронам (ННКт), малый зонд, усл. ед.; TP1 - время пробега Р-волны по короткому зонду АКШ, мкс; TP2 - время пробега Р-волны по длинному зонду АКШ, мкс; dGK - разностный параметр по ГК; W - водородосодержание, %; DS - измеренный диаметр скважины, м; Kгл - коэффициент глинистости, % (по ГИС); KпАК - коэффициент пористости по АК, %. Здесь и в табл. 2, 5, 6 красным цветом показаны значения, имеющие уровень вероятности ошибки менее 5 %. Таблица 2 Значения t-критерия Стьюдента по данным исследования керна Метод Среднее Значение критерия Стью-дента Число степеней свободы Уровень значи-мости Количество наблюдений Стандартное отклонение Относи-тельная дисперсия Дис-персия непрони-цаемый прони-цаемый непрони-цаемый прони-цаемый непрони-цаемый прони-цаемый Kп, % 10,75 18,68 -13,08 161 0 71 92 4,6 3,2 2,1 0,0 Kпр, мД 68,71 904,06 -6,75 140 0 50 92 206,0 860,9 17,5 0,0 Kво, % 45,25 9,68 6,68 54 0 20 36 29,2 10,0 8,5 0,0 ПлО 2,36 2,14 10,84 161 0 71 92 0,2 0,1 3,3 0,0 ПлН 2,47 2,36 7,74 161 0 71 92 0,1 0,1 6,3 0,0 ПлМ 2,64 2,63 0,79 161 0,43 71 92 0,1 0,0 76,6 0,0 Примечание: Kп - коэффициент пористости, %; Kпр - коэффициент проницаемости, мД; Kво - коэффициент остаточной водонасыщенности, %; ПлО - объемная плотность керна, г/см3; Плн - плотность насыщенного керна, г/см3; Плм - минералогическая плотность керна, г/см3. а б Рис. 1. Точечная диаграмма: а - построенная по керновым данным для Kв и Kп; б - для параметров GK и W Анализируя геофизические параметры, отметим, что наибольшая степень дифференциации наблюдается по показателям собственной радиоактивности (GK) и водородосодержанию (W) (см. рис. 1, б). Рисунок показывает, что проницаемые породы имеют низкое водородосодержание и низкие значения естественной радиоактивности, а непроницаемые породы - высокое водородосодержание и высокие показатели радиоактивности [15]. Уравнение регрессии для проницаемых пород имеет вид W = 22,9 - 0,5GK, для непроницаемых W = 10,7 + 1,6GK. Значения коэффициентов Z и свободных членов показывают, что они значительно различаются. Вероятностно-статистическая проверка разделения разреза на проницаемые и непроницаемые части Для определения качества разделения толщи пород составлена база по результатам интерпретации геофизических исследований и данным кернового материала. Из керновых данных использовались: коэффициент пористости, проницаемости, остаточной водонасыщенности, объемная плотность образца, плотность насыщенного керна и минералогическая плотность керна. Из геофизических методов учитывались данные, полученные с использованием электрических, электромагнитных, радиоактивных и акустических методов исследований и измеренный диаметр скважины. Во внимание принимались и параметры, полученные в результате вычислений: двойной разностный параметр по гамма-каротажу, коэффициент глинистости и пористость, рассчитанная по акустическому и нейтронному каротажам [16-18]. Сравнение влияния показателей на разделение разреза на проницаемые и непроницаемые части затруднительно из-за различной размерности показателей, учитываемых при дифференциации пород. Для приведения показателей к единой размерности необходимо построить линейные вероятностные модели для расчленения разреза по отдельным показателям. Построенные модели должны обладать следующими свойствами: · средние значения в группах должны находиться: для покрышек - в диапазоне значений от 0,5 до 1, для коллекторов - в диапазоне от 0 до 0,5. Граничным значением для разделения принималось значение 0,5; · средние значения должны быть расположены симметрично относительно значения 0,5 [19-21]. По керновым данным для всех параметров, кроме минералогической плотности, удалось построить модели, разделяющие проницаемые и непроницаемые части разреза. Модель, описывающая разделение разреза на проницаемые и непроницаемые породы по данным, полученным при исследовании пористости образцов керна, имеет обратную зависимость, что обозначает отнесение образцов с высокой пористостью к классу коллекторов; образцы, имеющие низкую пористость, относятся к классу покрышек. Таблица 3 Модели для нормировки данных, полученных при исследовании керна Показатель Линейные модели вероятности принадлежности к пластам-покрышкам Показатель Линейные модели вероятности принадлежности к пластам-покрышкам Kп P(Kп) = 0,955794 - 0,0311Kп Пло P(Пло) = -1,178 + 0,75Пло Kпр P(Kпр) = 0,544 - 0,00009Kпр Плн P(Плн) = -0,22 + 0,3109Плн Kв P(Kв) = 0,3412 + 0,0062Kв Плм Модель построить нельзя Аналогичная ситуация и с проницаемостью: коллекторы имеют высокие значения проницаемости, а покрышки - низкие. Модель, описывающая коэффициент водонасыщенности, имеет прямую зависимость, из чего следует, что покрышки имеют более высокую степень насыщения водой, чем коллекторы. Это связано с тем, что покрышками в терригенной части разреза являются в основном аргиллиты темно-серые, имеющие высокий объем закрытой пористости, с высокой степенью поглощения жидкости, благодаря которой и сохраняется высокая остаточная водонасыщенность. Прямая модель распределения объемной плотности образца и плотности насыщенного керна показывает, что покрышками являются более плотные породы, в то время как более разуплотненные являются коллекторами. Исходя из коэффициентов, полученных при построении статистических моделей плотности образца, видно, что при наполнении образца жидкостью разница в плотности между коллекторами и покрышками снижается относительно сухого образца. И все же испытания образцов керна проводятся точечно и в основном в интервалах продуктивных пластов, поэтому для получения более полной информации привлекались данные геофизических исследований. Аналогично с керновыми данными вероятностно-статистические модели построены и для геофизических параметров. Таблица 4 Модели для нормировки данных, полученных по ГИС Показатель Линейные модели вероятности принадлежности к пластам-покрышкам Показатель Линейные модели вероятности принадлежности к пластам-покрышкам BK P(BK) = 0,457 + 0,005BK NKTS P(NKTS) = 0,701 - 0,1369NKTS BMK P(BMK) = 0,465 + 0,004BMK TP1 P(TP1) = -0,194 + 0,00102TP1 DTP P(DTP) = -0,153 + 0,00259DTP TP2 P(TP2) = -0,215 + 0,00089TP2 GK P(GK) = 0,194 + 0,0353GK dGK P(dGK) = 0,425 + 0,17774dGK IK P(IK) = 0,769 - 0,0007IK W P(W) = 0,332 + 0,0069W MGZ P(MGZ) = 0,65 - 0,1914MGZ Ds P(Ds) = -0,23 + 3,1646Ds MPZ P(MPZ) = 0,379 + 0,1467MPZ Kгл P(Kгл) = 0,254 + 0,0096Kгл NKTB P(NKTB) = 0,589 - 0,0221NKTB KпАК P(KпАК) = 0,331 + 0,0072KпАК Электрические методы. Методы бокового и микробокового каротажа работают одинаково для выявления проницаемых и непроницаемых пропластков: чем выше сопротивление, тем выше вероятность выделения непроницаемого прослоя. Микроградиент- и микропотенциал-зонды работают в противофазе, т.е. при исследовании меньшего радиуса (радиус исследования микроградиент-зондом примерно равен 3,75 см) высокое сопротивление показывает наличие глинистой корки и присутствие в этом интервале коллектора. При исследовании более удаленной зоны (радиус исследования микропотенциал-зондом в 2,0-2,5 раза больше) на показания не влияет наличие глинистой корки, и с увеличением сопротивления растет вероятность прогнозирования непроницаемых пропластков в этом интервале. Индукционный каротаж имеет обратную зависимость: при увеличении значений снижается вероятность нахождения в этом интервале непроницаемых пород. Радиоактивные методы. Показания гамма-каротажа имеют прямую зависимость от вероятности определения непроницаемых пород. Это объясняется тем, что в основном непроницаемые пропластки связаны с наличием в них глин, которые имеют высокую естественную радиоактивность. Также это доказывается наличием обратной связи при определении вероятности прогнозирования непроницаемых пропластков по нейтронному каротажу. Акустические методы. Наличие прямой связи при определении вероятности нахождения непроницаемых прослоев показывает, что более плотные породы с большей долей вероятности относятся к непроницаемым, чем разуплотненные. Диаметр скважины в коллекторах уменьшается из-за наличия глинистой корки и увеличивается вследствие промытия глинистых интервалов, что демонстрирует прямая зависимость вероятностной оценки нахождения непроницаемых пропластков от диаметра скважины. Водородосодержание показывает, что в непроницаемых пропластках содержание водорода выше, чем в коллекторах. Разностный параметр ГК, как и сам метод гамма-каротажа, свидетельствует о прямой зависимости нахождения непроницаемых пропластков от значений параметра. Коэффициент глинистости подтверждает связь между наличием непроницаемых пропластков и глинистостью породы. При определении пористости по акустическому каротажу наблюдается обратная зависимость. Это связано с тем, что непроницаемая часть терригенного разреза состоит в основном из аргиллитов, имеющих высокий объем пустотного пространства при низком сообщающемся. Использование вероятностно-статистических моделей позволило привести все параметры к единой системе измерения. С помощью коэффициента Стьюдента проводилась проверка правильности выбора математических моделей и анализа степени влияния каждого из параметров на выделение проницаемой и непроницаемой частей (табл. 5). Из данных табл. 5 видно, что средние значения классов находятся в допустимых интервалах коллектора 0,0-0,5; покрышки - в интервале 0,5-1,0. Исключение составляет параметр P[e1] (ПлНас) по керновым данным, что объясняется способностью глинистых пород впитывать жидкость насыщения. Анализ степени влияния на разделение пород Таблица 5 Проверка правильности выбора математических моделей и анализа степени влияния каждого из параметров на выделение проницаемой и непроницаемой частей Параметр Среднее Значение критерия Стьюдента Число степеней свободы p Количество наблюдений Стандартное отклонение Относи-тельная дисперсия p класс 2 класс 1 класс 2 класс 1 класс 2 класс 1 Параметры, определенные по керну Р(Kп) 0,63 0,37 13,17 162 0 72 92 0,14 0,10 2,12 0 Р(ПлО) 0,58 0,43 8,66 162 0 72 92 0,15 0,07 5,04 0 Р(Kпр) 0,53 0,46 6,75 140 0 50 92 0,02 0,08 17,46 0 Р(Kв) 0,62 0,40 6,68 54 0 20 36 0,18 0,06 8,53 0 Р(ПлН) 0,54 0,51 5,06 162 0 72 92 0,06 0,02 12,22 0 Геофизические параметры Р(GK) 0,66 0,34 26,83 353 0 208 147 0,14 0,06 5,82 0 Р(dGK) 0,55 0,45 26,83 353 0 208 147 0,04 0,02 5,82 0 Р(Кгл) 0,67 0,34 26,83 353 0 208 147 0,14 0,06 5,82 0 Р(IK) 0,64 0,34 13,90 345 0 200 147 0,10 0,27 6,68 0 Р(NKTS) 0,53 0,46 11,38 353 0 208 147 0,07 0,02 10,05 0 Р(W) 0,55 0,48 8,93 353 0 208 147 0,09 0,02 19,07 0 Р(DS) 0,53 0,46 7,38 353 0 208 147 0,12 0,00 1845,21 0 Р(DTP) 0,53 0,47 7,18 353 0 208 147 0,11 0,02 19,86 0 Р(КпАК) 0,52 0,48 7,18 353 0 208 147 0,07 0,02 19,86 0 Р(NKTB) 0,52 0,48 7,17 353 0 208 147 0,06 0,02 9,22 0 Р(TP2) 0,52 0,48 6,43 353 0 208 147 0,08 0,03 6,96 0 Р(BMK) 0,52 0,48 6,20 353 0 208 147 0,07 0,03 6,51 0 Р(TP1) 0,52 0,48 5,32 353 0 208 147 0,07 0,03 4,74 0 Р(MGZ) 0,51 0,49 2,05 353 0,04 208 147 0,11 0,03 11,96 0 Р(BK) 0,50 0,49 1,50 353 0,13 208 147 0,05 0,08 2,31 0 Р(MPZ) 0,50 0,50 0,55 353 0,58 208 147 0,08 0,03 4,72 0 Примечание. Класс 1 - проницаемые горные породы; класс 2 - непроницаемые горные породы. показал, что из геофизических параметров наибольшее влияние оказывает значение естественной радиоактивности пород, а наименьшее - показания микропотенциал-зонда. Из параметров, определяемых на керне, наибольшее влияние оказывает коэффициент пористости, а наименьшее - плотность насыщенной породы. Используя результаты, полученные с помощью статистических моделей, приведенных в табл. 3 и 4, вычислим условную комплексную вероятность по следующей формуле: где P1 … Pi - индивидуальные вероятности. Для комплексного анализа построим многомерные модели при различных сочетаниях m. Первую вероятность построим при m = 2. При Таблица 6 Коэффициент Стьюдента t при включении в расчет комплексной вероятности различных параметров Параметр m = 2 m = 3 m = 4 m = 5 m = 6 m = 7 m = 8 m = 9 m = 10 m = 11 m = 12 m = 13 m = 14 m = 15 m = 16 Параметры, определенные по керну Р(Кп) + + + + Р(ПлО) + + + + Р(Кпр) + + + Р(Кв) + + Р(ПлН) + Среднее, класс 1 0,67 0,72 0,74 0,67 Среднее, класс 2 0,29 0,25 0,26 0,29 t 12,07 8,60 8,65 12,07 p 0,33 0,01 0,03 0,33 Геофизические параметры Р(GK) + + + + + + + + + + + + + + + Р(dGK) + + + + + + + + + + + + + + + Р(Кгл) + + + + + + + + + + + + + + Р(IK) + + + + + + + + + + + + + Р(NKTS) + + + + + + + + + + + + Р(W) + + + + + + + + + + + Р(DS) + + + + + + + + + + Р(DTP) + + + + + + + + + Р(КпАК) + + + + + + + + Р(NKTB) + + + + + + + Р(TP2) + + + + + + Р(BMK) + + + + + Р(TP1) + + + + Р(MGZ) + + + Р(BK) + + Р(MPZ) + Среднее, класс 1 0,69 0,77 0,82 0,82 0,82 0,81 0,80 0,79 0,79 0,79 0,80 0,79 0,78 0,78 0,79 Среднее, класс 2 0,30 0,19 0,17 0,15 0,14 0,13 0,12 0,11 0,10 0,10 0,10 0,10 0,09 0,10 0,11 t 27,75 29,66 29,58 28,38 27,68 27,57 25,63 24,77 24,43 24,11 24,63 24,12 24,01 23,07 23,25 p 0,00 0,00 0,69 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 расчетах с различным значением m используется такое сочетание, при котором средние значения в классах коллекторов и покрышек максимально различаются. Степень различия определялась по критерию Стьюдента t (табл. 3), результаты расчетов приведены в табл. 6. По керновым данным наибольшее различие наблюдается при m = 5, при котором в расчете участвуют параметры пористости, проницаемости, остаточной водонасыщенности, плотности образца и плотности насыщенного образца. Таким образом, все показатели, рассмотренные в работе, увеличивают степень дифференциации пород. По геофизическим данным наибольшее различие наблюдается при m = 12, при котором в расчете участвуют практически все изучаемые параметры, кроме TP1, MGZ, BK и MPZ. При включении этих параметров в расчеты степень дифференциации уменьшается (см. табл. 5). Получив значения, при которых распределение по геофизическим и геологическим Рис. 2. Точечная диаграмма для геологических и геофизических параметров параметрам будет максимальным, построим график распределения параметров геологии в зависимости от геофизических данных (рис. 2). Анализ диаграммы показывает, что значения, принадлежащие к классу проницаемых пород, находятся в основном в диапазоне от 0 до 0,4 по геофизическим параметрам и в диапазоне - от 0 до 0,6 по геологическим параметрам, за исключением единичных образцов. Класс непроницаемых пород находится в пределах значений от 0,4 до 1,0 по геологическим параметрам и от 0,5 до 1,0 по геофизическим параметрам, за исключением единичных образцов. Поля точек, относящихся к различным классам, не перекрываются, что позволяет сделать вывод о действенности метода для разделения толщи на проницаемую и непроницаемую части. Заключение При выполнении работы обработан материал по скважине Софьинского месторождения. Проведен анализ методики выделения пород-коллекторов. Выявлены наиболее значимые параметры, влияющие на дифференциацию разреза. Построены линейные модели для приведения различных размерностей геофизических методов и данных исследования керна к единому знаменателю. Проведена градация геофизических методов и геологических данных по степени влияния на расчленение разреза. Рассчитан комплексный вероятностный параметр при вовлечении различного количества переменных. Использование керновых данных увеличивает степень дифференциации пород, геофизические параметры в большинстве своем также увеличивают дифференциацию пород, однако некоторые показатели, наоборот, не позволяют разделить породы на проницаемые и непроницаемы. Благодаря данному методу можно наиболее полно провести разделение толщи пород на проницаемую и непроницаемую части.

About the authors

Aleksandr V. Shcherbenev

Author for correspondence.
Email: Aleksandr.Shcherbenev@pnn.lukoil.com
9 Geroev Khasana st., Building A, Perm, 614066, Russian Federation

2nd category engineer of the Department of Reserves Estimation

References

Statistics

Views

Abstract - 302

PDF (Russian) - 38

PDF (English) - 79

Refbacks

• There are currently no refbacks.