Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Электроразведка при поисках нефти и газа - курс...doc
Скачиваний:
28
Добавлен:
24.11.2019
Размер:
11.46 Mб
Скачать

5. Основные способы вычисления кривых кажущегося сопротивления в зст.

Первичные материалы при выполнении зондирования становлением поля представляются в двух основных формах: кривые кажущегося сопротивления и кривые кажущейся продольной проводимости. При вычислении кажущегося сопротивления возможно использование приближенных соотношений типа (15) и (17). Вместе с тем, используется расчет кажущегося сопротивления и по точным формулам. Принцип такого вычисления рассмотрим на примере азимутальной составляющей напряженности электрического поля, регистрируемого в поле вертикального магнитного диполя. Величина, обратная параметру становления u, может быть записана в виде: 1/u = τ /2πt. Следовательно,

= = =

Преобразуем коэффициент в (12):

= = = .

Теперь уравнение (12) можно записать в виде:

. (18)

Здесь функция еφ(u) зависит только от параметра становления u и ее можно предварительно представить в виде таблицы или графика. Таким образом, если по измеренным значениям Eφ(t) для каждого времени регистрации, пользуясь (18) вычислить значения еφ(u), то, пользуясь заранее составленными таблицами или графиками, можно найти соответствующие значения u и далее вычислить кажущееся сопротивление:

= r20/2tu2 (19)

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

F(u) = .

С помощью таблиц и номограмм по F(u) находятся значения u и далее ρm. Приближенная формула для вычисления кажущегося сопротивление имеет вид:

.

Если измерения выполняются совмещенными установками, то вначале вычисляется, f(u) =

Далее находится u2 и затем вычисляется кажущееся сопротивление:

ρτ = .

Здесь Rрадиус круглой петли. При другой форме петли в качестве R рассматривается величина R = .

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

. (20)

здесь К = - геометрический коэффициент:

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

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

. (21) Здесь hi Si - глубины до центров отдельных слоев и их продольные проводимости. Так как при каждом времени становления достигается определенная глубина исследования, толще пород, в которой при данном времени проходит этот процесс, может быть поставлена в соответствие эквивалентная плоскость, расположенная на глубине hτ и имеющая продольную проводимость Sτ. Принцип вычисления параметров эквивалентной плоскости проследим на примере становления азимутальной составляющей электрического поля вертикального магнитного диполя. По В.А.Сидорову

, (22) где , K = πR3/3I, (23)

Отсюда величина S может быть определена как

S = F(m)/KE . (24)

Для определения F(m) авторами предложена следующая схема.

  1. По Еφ проводится численное дифференцирование, т.е. для выбранных моментов времени t считается значение производной .

  2. Для каждого выбранного момента времени считается функция φ(m) = . (25)

3. По заранее составленным для φ(m) таблицам или номограммам определяется параметр m и далее F(m).

4. По формуле (24) вычисляется Sτ, а так же

hr = mr - . (26)

Выражения типа (23) - (26) являются стандартными при расчетах кривых Sτ . Следует лишь иметь в виду, что вместо Еφ при расчетах следует использовать конкретный измеряемый сигнал; при этом для каждой разновидности установки будет различен и коэффициент К.