§ 119. ПРЕДВАРИТЕЛЬНАЯ ОБРАБОТКА ИСХОДНЫХ ДАННЫХ

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

1T . Подготовка исходной информации, заключающаяся в просмотре ее с целью дальнейшего подбора геофизических, геологических и других данных, необходимых для решения поставленной задачи (в частности, сведения о результатах исследования керна, опробования коллекторов, об испытании скважин и вероятном положении газожидкостных контактов, техническом состоянии скважины и другие). Эти данные вводятся в информационно-поисковую систему (ИПС) машины на магнитные барабаны (МБ) или ленты (МЛ) (банк данных).T

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

§ 120. ИНДИВИДУАЛЬНАЯ ИНТЕРПРЕТАЦИЯ КРИВЫХ ГЕОФИЗИЧЕСКИХ ПАРАМЕТРОВ

ПриT точечной системе интерпретации кривых геофизических параметров после корректировки всех материалов осуществляют переход от геофизических величин ρT кTB B (ρT эфTB ),B UспB ,B IγTB ,TB InB γT TB измеряемых в скважинах, к определяющим их физическим свойствам пород ρT пBT B ,ρT ппBT ,B ∆τT пBT ,B δT ЦBT B и др. или

к производным геофизическим параметрам αT TBсп,B ∆T JT γT T ,∆T TJnγT T и др., по которым в последующем выполняют комплексную интерпретацию геофизических данных.T

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

следующими критериями:T

(T критерий целесообразности выделения объекта по различию физических свойств его вмещающих пород);T

(T критерий целесообразности выделения объекта по его мощности). В последних формулах αT BTmaxB и αT BTmlnB

— максимальное и минимальное значения изучаемого параметра в рассматриваемой области zaB max,B ZaB minB — глубины залегания соответствующих точек.T

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

сT заранее принятой предельной величиной λT 2TP кPB ,B грB B (здесь ρT кTB , iB — экстремальные значения ρT кTB B в пропластках, слагающих пачку; H=ΣT iBT hB iB )B В дальнейшем для количественного изучения используются только те пропластки, для которых λT 2PT кBP , HB >λT 2PT кBP , прB . При λT 2PT кBP , НB < Tλ2PT кBP , прB (λT 2PT кBP , пр —B кажущееся предельное значение) вся пачка рассматривается как объект, для которого среднее значение TρкBT B находится в соответствии с его продольным сопротивлением TρnT .T

ДляT определения точного положения границ по симметричным кривым (кривые UCB ПB ,TρэфBT ,B ∆τT T и др.) машина, считывая последовательно значения изучаемого параметра аiP jBP B (j — регистрируемый

параметр UCB ПB ,IγT ,T InB γT ,TB ρT кTB ,B ρT эфTB B и др.), вычисляет разности ∆T TaiP jPB B = ajP i+1PB -B ajP iPB B . Далее в области монотонного изменения (возрастания или убывания) ∆T aT iP jBP B определяется глубина точки, в которой ∆T aT iP jBP B максимально. Затем то же осуществляется в области монотонного изменения ∆T aT iP BPjB противоположного знака

279

(убывания или возрастания); устанавливается и выдается на печать глубина залегания следующей

точки |∆T a|T iP jBP max.B Разностью глубин указанных точек определяется кажущаяся мощность пласта.T

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

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

здесьT aT ijB TB = ln ( ρT кTB /B ρT рTB B )ijB ;B j — определяет порядковый размер используемогоT зонда, выраженный в единицах диаметра скважины (L/dcB );B т — общее число зондов (отношений L/dcB B , заложенных в память

машины); i —шифр палеточной кривой — искомые значения ρT кTB /B ρT рTB ;TB п — общее число теоретических кривых, образующих матрицу. Шифром матриц являются модули ρT зпBT /B ρT рBT .B > и D/dcB .B В машину вводится вектор интерпретируемой кривой зондирования — последовательность соотношений ln (ρT кTB /B ρT рTB )TB для используемых L/dcB .B Критерием решения задачи является минимум среднеквадратичного расхождения между логарифмами ординат интерпретируемой и теоретической двухслойных кривыхT

ПриT работе с палетками U-эквивалентных кривых аналогичным способом находят параметр U-эквивалентности.T

ПриT определении электрического удельного сопротивления по данным индивидуальных

кривых ρT кTB TB и TρэфTB ,B зарегистрированных с зондами различных типов, предварительно в значения TρкTB TB (ρT эфTB )B вводят поправки за ограниченную мощность пласта и сопротивление вмещающих пород и по

исправленным значениям в дальнейшем вычисляют искомые параметры TρлBT ,B ρT ппBT B или ρT зиBT .B T

ПриT автоматической интерпретации кривых потенциалов собственной поляризации с помощью машины вводят поправки в ∆T UT спB заB различие удельных сопротивлений ρT лBT B, ρT зпBT ,B ρT рBT ,B ρT вмBT иB ограниченную мощность пласта. Возможен также вывод разности T∆AT даB B = AдаB , пB - AдаB , вмB ,

пропорциональной ∆T UT спB .B T

ИнтерпретируяT данные гамма-метода, в память машины вводят сведения о гамма-активности глинистого раствора, диаметре скважины, толщине глинистой корки, мощности пласта, произведения

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

ПриT интерпретации нейтронных методов по значениям hгкB B и данным о водородосодержании глинистого раствора машина вносит соответствующие поправки в InB γT ,BT In,B TB и In,B HB и рассчитывает

исправленные значения InB γT (ITB n,B T),B выраженные в единицах разностного параметра, в вероятностных единицах или в величинах коэффициента kпB , nB нейтронной пористости.T

§ 121. КОМПЛЕКСНАЯ ИНТЕРПРЕТАЦИЯ РЕЗУЛЬТАТОВ ГЕОФИЗИЧЕСКИХ ИССЛЕДОВАНИЙ СКВАЖИН

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

280

ЛитологическоеT расчленение разрезов скважин выполняют приближенно способом градаций и более точно — решением уравнений петрофизических связей (284) и (285).T

ВT способе градаций разрез скважины расчленяют по комплексу значений физических свойств исследуемых пород или, что менее точно, по комплексу числовых значений измеряемых в скважине геофизических параметров. Эти данные либо вводятся в машину на основании обобщения всех имеющихся петрофизических определений с учетом вероятного влияния на физические свойства пород термобарической обстановки, в которой находитсяT изучаемаяT порода в условиях ее естественного залегания, либо выдаются как результат предшествовавшего этапа индивидуальной интерпретации. Каждому из физических свойств или геофизических параметров отводится несколько разрядов двоичного кода, соответственно которым тот или иной физический параметр или свойство могут быть изучены по ряду диапазонов — градаций. Таким образом, при изучении литологического разреза скважины каждый из обрабатываемых параметров X разбивается на классы с граничными значениями АхB ,B iB -B (х — порядковый номер параметра, i — номер границы). В программе «Каротаж» две границы выделяют три градации, которые кодированы: малые значения X — 00, средние значения X — 01, большие значения X — 10.T

ГраничныеT значения Ax,B iiB предварительно устанавливают, исходя из дифференцированности

изучаемого параметра и его роли в определении литологии. Например для ρT кTB B (порядковый номер параметра «3»), зарегистрированного стандартным зондом, Н. Н. Сахранов рекомендует A3B ,B B 1B = (1,2 +

1,5ρT глBT )B Ом*м и A3*3B =B (150— 1200) Ом*м (в зависимости от значений ρT вмBT B и ρT вBT B [54]. Однако для интерпретации кривых сопротивлений, учитывая огромные пределы изменения этого параметра,

целесообразнее выделять три разряда.T

ПриT использовании двухразрядного кода интервал возможного изменения AдаB B (∆T UT cuB )B может быть разделен на следующие диапазоны: АдаB B < 10 мВ (∆T UT cuB B > 0,8∆T UT cuB ,B mB ах)B — код 00 — чистые и слабоглинистые породы; АдаB B = 10-30 "мВ [∆T UT cuB =B = (0,5 — 0,8) ∆T UT cB п, maxB] — код 01 — глинистые песчано-алевритовые и карбонатные породы; АдаB B = 30-50 мВ ∆T UT cB п=B = (0,2 — 0,5) ∆T UT cB п, maxB ]— код 10

— сильно глинистые песчано-алевролитовые и карбонатные породы; АдаB B > 50 мВ (∆T TUcB п<B < 0,2∆T UT cB п,

max)B — код 11 — глины, плотные мергели.T

ДляT кавернограмм два разряда могут включать следующие диапазоны изменения dcB B и

соответствующие им коды: dcB B < dпB B ± TδdTB B - 00; dcB =dB дB ±B δT dTB —B 01; dcB B = (1 — 1,5) dBдB ± δT dTB B — 10 и dcB B > 1,54dдB B —

11 (δT dTB B —погрешность в определении dcB )B .T

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

значения kпB ,nB ,δT пTB ,B ∆τT ,T ρT пT , TρппTB B и др. за влияние условий измерений в скважине, в зависимости от числа использованных методов осуществляет частичное или полное решение системы следующих

уравнений:T

РешениеT системы уравнений (461) дает возможность:T

1)T устанавливать вероятное нефтегазонасыщение коллекторов, определять kBn,B kглB ,B kнB B и kгB ;B T

2)T в случае сложного минерального состава скелета коллектора и достаточного числа уравнений выделять основные его компоненты — кальцит, доломит, ангидрит, гипс, кварц с учетом

характерных значений δT ,T T∆τT и An,B B пB B для этих минералов;T

3)T по сопоставлению kвB B с kвB , ппB устанавливать наличие неизвлекаемых углеводородов и тем

самым уточнять коэффициент вытеснения βT вBT .B T

РезультатыT комплексной обработки данных геофизических исследований скважин электронно-вычислительные машины выдают на печать в цифровой или графической форме (рис. 206).T

281

РисT . 206. Примеры графического изображения результатов непрерывной комплексной

интерпретации данных геофизических исследований скважин на ЭВМ (по Дрессер Атласу).

а — песчаный разрез; б — карбонатный разрез

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

ИспользованиеT электронно-вычислительных машин имеет большое значение при анализе различных методик интерпретации и выборе из них наиболее обоснованной. Это особенно важно в случае глинистых коллекторов, для которых в зависимости от принятой модели распределения глинистых частиц и их влияния на электрические и другие физические свойства горных пород могут быть получены различные связи между электрическими свойствами и коэффициентами пористости и нефтегазонасыщения [4, 11, 44].T

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

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

§ 122. ИСПОЛЬЗОВАНИЕ ЭВМ ПРИ РЕГРЕССИОННОМ АНАЛИЗЕ И РЕГИОНАЛЬНЫХ ИССЛЕДОВАНИЯХ

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

проницаемости, зольностью и др.) вида kпB B = f (ρT пTB /B ρT вTB ,B ρT ппTB /B ρT фTB ,B αT спTB ,B ∆T IT γTB ,TB T∆IT n,B γT TB , ∆τT пTB )TB ;T kT пB B = f (РT нB ,B (kвB ),B αT спTB ,B ∆T IT TBγ ,TB

∆T IT n,B γT BT , ∆τT пBT )B иT др. Использование многомерных связей при достаточной их обоснованности в ряде случаев позволяет уточнить интерпретацию как при машинной, так и при неавтоматической обработке результатов геофизических исследований скважин.T

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

ОдноT из новых направлений — использование ЭВМ при o6o6щающейT интерпретации результатов геофизических исследований скважин с целью получения наиболее объективной корреляции их разрезов (см. § 125) и изучения цикличности осадконакопления методами гармонического анализа (см. § 133).

282

ЧАСТЬ ТРЕТЬЯ ОБОБЩАЮЩАЯ ИНТЕРПРЕТАЦИЯ РЕЗУЛЬТАТОВ ГЕОФИЗИЧЕСКИХ ИССЛЕДОВАНИЙ СКВАЖИН

Глава XX

НОРМАЛЬНЫЕ И СВОДНЫЕ ГЕОЛОГО-ГЕОФИЗИЧЕСКИЕ РАЗРЕЗЫ. КОРРЕЛЯЦИЯ РАЗРЕЗОВ СКВАЖИН

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

СводныйT геолого-геофизический разрез, составленный по средним истинным мощностям пластов, горизонтов и свит в их нормальном залегании, называется нормальным геологогеофизическим разрезом (рис. 207).T

Рис. 207. Нормальный геолого-геофизический разрез

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

§ 123. СОСТАВЛЕНИЕ НОРМАЛЬНЫХ И СВОДНЫХ ГЕОЛОГОГЕОФИЗИЧЕСКИХ РАЗРЕЗОВ

НормальныеT геолого-геофизические разрезы по месторождению составляют по истинным мощностям в «нормальной» по возрасту последовательности.T

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

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

углом а в азимуте θT ,T истинные мощности пород подсчитывают по значениям видимых мощностей по формуле (288) или с помощью специальных номограмм [2].T

283

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

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

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

1T . Для нормальных диаграмм кажущегося сопротивления и потенциалов собственной поляризации используют кривые, зарегистрированные в скважинах, заполненных глинистыми растворами наиболее часто встречающегося (нормального для данного района) удельного сопротивления. Для нефтеносных коллекторов на диаграммы ρT кBT B наносят две предельные кривые, соответствующие продуктивному и непродуктивному состоянию коллектора. Площадь, ограниченную этими кривыми, заштриховывают. В заголовке нормальной диаграммы указывают тип и размер зонда, диаметр dcB B и удельное сопротивление рвB ,B а у нормальной диаграммы — UcB п,B полярность диаграммы и значение ρT рBT .B T

2T . Нормальные диаграммы гамма-метода и нейтронных методов составляют по кривым, зарегистрированным в необсаженных скважинах однотипной аппаратурой. В заголовке диаграммы IγTB ,TB записанных в гаммах, указывают диаметр скважины dcB ,B тип прибора, тип и число индикатора, постоянную времени τT яTB ,B скорость регистрации, радиоактивность kγTB ,T PB (qγTB ,T P)B глинистого раствора, уровень фона. В заголовке диаграмм lnB ,B InB γT ,TB вычерченных в единицах условных, двойного разностного параметра или в нормализованных, записывают тип прибора (тип и число индикаторов), интенсивность источника значения dcB ,B τT яTB ,B v и уровень фона. На диаграмме отмечают интенсивности излучений, соответствующие опорным пластам разреза.T

3T . Нормальные диаграммы интервального времени ∆τT пBT B и коэффициента затухания упругих волн составляют по кривым, зарегистрированным однотипной аппаратурой. В заголовке указывают тип аппаратуры, размер зонда, мощность излучателя, частоту излучения и подачи импульсов, значение dcB .B T

4T . В заголовке нормальной термограммы t (или ∆τT )T отмечают диаметр скважины dcB ,B вязкость и термические свойства глинистого раствора, скорость регистрации и температурную инерцию τT tTB B термометра. На термограмму наносят точку А равенства температур tpB B и tпB .B T

5T . Для нормальной диаграммы τT прBT B используют кривые, зарегистрированные при наиболее часто применяемых типах и размерах долот, постоянных нагрузках на долото и частотах вращения. В заголовке диаграммы записывают тип и диаметр долота, нормальные нагрузки Р и частоту вращения долота, для которых составлена диаграмма τT прBT .B

На всех кривых должны быть даны масштабы их регистрации.

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

284