Файл: Сохранов Н.Н. Машинные методы обработки и интерпретации результатов геофизических исследований скважин.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 11.04.2024
Просмотров: 137
Скачиваний: 0
Начало |
Ввод Ванных |
|
|
|
|
|
Печать |
^)~*] |
Коней |
|
|
|
Вычисление £ |
|
ParPc Ра?Я ß W
Рп'ИРпМ
Определение Ajuß:
Определение AJUBJ
X
Ц1+бг)-а,о55]р°И-
Нет |
n(0) |
|
PnfPcA>PnfPcV |
||
|
Рп={(Рп,+Яг)
(Рэп/Яш»*® h
|
W-№^WPc\ |
|
РпТРс^РпгРсВ'» |
Определение Aj и Bj |
Определение КЗП [*, |
|
|
1 |
Нет |
Определение AjuBj
Рис. 63. Блок-схема универсальной программы определения удельного сопротивления пластов н зоны проникновения.
V — номер пласта; JV — число пластов; I — число зондов.
граничное значение а) указывает на соответствие условных модулей
сопротивлений |
фактическим |
модулям. |
Результат |
интерпретации |
|
(рп и |
р п 1 -грпг) |
уточняют при меньших |
значениях |
о;. |
|
3. |
Проверка |
пересечения |
интервалов |
[А-, Bj], определенных по |
|
трехслойным матрицам (при небольшой мощности пласта для зондов
с £ > Я |
применяют матрицы |
ЭКЗ): |
|
|
|
||
а) выбор матрицы по величине РзП /рс , найденной с помощью дву |
|||||||
слойной |
матрицы; |
[A-, Bj] |
|
|
|
||
б) определение интервалов |
для |
строк |
выбранной ма |
||||
трицы |
при Oj — ог; |
|
|
|
|
||
в) |
проверка |
пересечения интервалов; |
|
|
|
||
г) |
уточнение |
р п при меньших значениях |
о;.; |
|
|||
д) |
уточнение |
p j n . |
|
удельного |
сопротивления |
||
Рассмотрим |
алгоритм определения |
||||||
пластов и зоны проникновения на ЭВМ по комплексу зондов (в том числе разнотипных), показанный на рис. 63.
1. Ввод исходных данных:
а) палеточные данные •— двуслойные и трехслойные кривые Б К З , правые ветви экстремальных кривых зондирования и поправочные коэффициенты на мощность, представляемые в виде описанных выше
таблиц |
(матриц); |
|
|
|
|
|
|
|
|
||
б) данные об интерпретируемом пласте (d, |
H, р в м , рс и исправлен |
||||||||||
ные на |
влияние |
мощности значения рк /рс )- |
|
|
|
||||||
2. Определение интервалов модулей кажущегося удельного со |
|||||||||||
противления (1 — 0Г) р к у/рс |
и (1 + |
ог ) р к |
у/рс для зондов /; стг обычно |
||||||||
равно |
|
0,2. |
|
|
|
|
|
|
|
|
|
3. Определение интервалов условных модулей удельного сопро |
|||||||||||
тивления руу -/рс |
в предположении, |
что среда двуслойная, и анализ |
|||||||||
их. Границы Aj |
И B J интервалов |
определяют линейной интерполя |
|||||||||
цией |
по параметру |
р п /р с |
между строками таблицы |
двуслойной |
па |
||||||
летки |
|
для (1 — ог ) |
р к / / р с |
и |
(1 + |
огг) рк //Рс- |
W/, Bj] |
условных |
мо- |
||
4. |
Проверка |
пересечения |
интервалов |
f] |
|||||||
|
|
|
|
|
|
|
;=i,...,fc |
|
|
|
|
дулей. Для этого определяют максимальное значение А1-0* из верх
них границ Aj |
и минимальное Б<о) из нижних границ Bj |
интервалов |
||||||
условных модулей удельного |
сопротивления [23]. Удельное сопро |
|||||||
тивление пластов будет находиться в интервале рсп°і — рп °2 , где рп і |
= |
|||||||
= |
РсАМ; Р п ° 2 |
= |
Р с Б(о). |
интервалов |
условных модулей (рпг |
> |
||
|
5. В |
случае |
пересечения |
|||||
> |
Рп°і) |
уточняют |
интервал |
рп1 — р п 2 . |
С этой целью |
повторяют |
||
5 раз операции по пунктам 2—4, уменьшая каждый раз интервалы
кажущегося сопротивления для |
всех зондов; границы |
интервалов |
ßy и bj в этом случае вычисляют |
по формулам |
|
ау = [(1 —crr) + 0,05-s] р к / ; |
Ъ, = [(1 +<тг ) - 0,05 |
• s] р к ( 8 9 ) |
6. Определяют интервал вероятных удельных сопротивлений пласта р п 1 = pc ^.( s \ рП 2 = РсВ^', за удельное сопротивление при нимают среднее значение, определенное по формуле (82). Если при
157
определении удельного сопротивления был использован |
хотя |
бы |
|||||||
один зонд с большим радиусом исследования, то оценивают |
погреш |
||||||||
ность |
е определения |
р п по формуле (88). |
|
|
|
|
|
||
7. |
Запоминают и печатают номер пласта ѵ и значения р п |
1 , р п |
2 , р п |
||||||
и е (или признак пересечения). |
|
|
|
|
|
|
|||
8. Отсутствие пересечения указывает на трехслойный |
характер |
||||||||
кривой БК З или на увеличение диаметра скважины |
против |
пласта |
|||||||
(d>dK). |
|
|
|
|
|
|
|
|
|
9. |
Определение |
приближенного |
значения |
р з п / р С і |
если |
d |
> |
dn. |
|
За это значение принимают интервал пересечения условных |
модулей |
||||||||
удельного |
сопротивления, полученных для малых и средних |
зондов |
|||||||
(например, |
зонды БК З с АО = 0,45; 1,05 и 2,25 м). |
|
|
|
|
||||
10. Выбор трехслойной палетки (матрицы) для интерпретации. |
|||||||||
При глубоком проникновении (пп > |
3, где пп |
— число зондов с не |
|||||||
большим |
радпусом |
исследования, |
интервалы |
р у ; |
- / р с |
которых пере |
|
||||||||
секаются |
при определении |
р з п |
/ р с п 0 |
двуслойной |
матрице) шифр па |
|
|||||||||
летки определяют из условий (85) или (83) и (84). В других |
случаях |
|
|||||||||||||
шифр палетки определяют из условий (86) или (83) и (84). |
|
|
|
||||||||||||
11. |
Определение интервалов условных модулей для кривых (строк |
|
|||||||||||||
выбранной матрицы) с различными значениями |
Did. |
|
|
|
|
|
|||||||||
12. |
Проверка |
пересечения |
интервалов |
|
условных |
моду- |
|
||||||||
лей |
Л [А,, |
В;]. |
|
|
|
|
|
|
|
|
|
|
|
|
|
13. |
В случае пересечения уточняют интервал |
р п 1 |
— р п 2 |
и |
р |
||||||||||
так же, как и в случае двуслойной |
среды (см. пункт 5). |
|
|
|
|
||||||||||
14. |
При-D/c? < |
4 определяют удельное сопротивление р п |
1 |
= p c ^ l ( s ) |
|
||||||||||
и р п 2 |
= |
р с 5 М пласта и параметры Did |
= (D/d)^ |
и |
р з п |
= р с |
( р з п / р с ) ш |
|
|
||||||
зоны |
проникновения. |
|
|
|
(Did >> 4) вычисляют |
|
|
|
|||||||
15. |
При глубоком проникновении |
отноше |
|
||||||||||||
ние р ; . п |
/ р с |
п определяют удельное сопротивление зоны проникновения |
|
||||||||||||
по формуле (87). Диаметр |
зоны проникновения |
в |
обоих |
случаях |
|
||||||||||
(см. пункт 14) определяют |
по параметру кривой |
(строка |
матрицы), |
|
|||||||||||
для |
которой зафиксировано |
пересечение. |
|
|
|
|
|
|
|
||||||
16. |
Определяют |
р п по формуле (82) и погрешность е по формуле (88). |
|
||||||||||||
17. |
Запоминают и печатают номер пласта ѵ и значения |
р п 1 , |
р |
п |
|||||||||||
рп , р з п , Did и е (или признак пересечения).
Программа определения удельного сопротивления пластов ши роко опробована в различных районах. Оценка результатов машин ной интерпретации данных Б К З в скважинах различных месторожде ний производилась путем сопоставления их с данными квалифициро ванной ручной интерпретации (рис. 64 и 65). Определения величины р п обоими способами считались достаточно близкими, если выпол нялось одно из следующих условий:
'S € [рп X. Р п 2І> |
(90) |
ss 0,2. |
(91) |
158
ЭВМ
ßn, Он-м
10000
/ООО
Рис. 64. Сопоста вление результа тов определения р п
на ЭВМ ( р ^ в м ) и ручным способом рР по скв. 102 од ного из районов Украины (по Д. Е.
Коваленко).
1000 р р ' о н м
Удельные сопротивления зоны проникновения, полученные на ЭВМ и ручным способом, считались приемлемыми в том случае, если они удовлетворяли условию (91) или расхождения их могли быть объяснены эквивалентностью решений.
Сопоставление показало, что результаты машинной ЭВМ и ручной интерпретации, как правило, хорошо согла суются между собой. Число расхождений между вели чинами Рп, полученными обоими способами, не пре вышает 18% от числа сопо ставленных пластов. Поло вина расхождений связана
Рис. 65. Сопоставление результатов определения р з п на ЭВМ (pf^M) и ручным способом (рР ).
159
с ошибками ручной интерпретации. Другая половина |
объясняется |
|||
плохим качеством данных Б К З |
(разброс |
точек |
кривой |
зондирова |
ния), влиянием ограниченной |
мощности |
пласта, экранированием |
||
соседним пластом н Г л у б о к и м |
проникновением |
раствора, которые |
||
трудно учитывать как при ручной, так и при машинной интерпре тации [40].
Определенные на ЭВМ удельные сопротивления зоны проникнове ния проницаемых пластов отличаются от результатов ручной интер
претации значительно больше и чаще, чем значения р п . |
Однако, как |
|
видно из рис. 65, определенные |
на ЭВМ значения pf B M |
равномерно |
отклоняются от значений р | п , |
полученных ручным способом, как |
|
в сторону уменьшения, так и в сторону увеличения. Это свидетель
ствует о том, что алгоритм выбора |
палетки удовлетворителен. |
|||||
Универсальная программа определения удельного сопротивления |
||||||
проверена также при интерпретации данных комплекса |
разнотипных |
|||||
зондов: а) |
малые градиент-зонды с АО = 0,25 и 0,45 м, |
зонд |
АБК-3 |
|||
и зонд индукционного каротажа 4ФО,75; |
б) неполное |
Б К З |
и з о н д |
|||
ИК-2; в) малый градиент-зонд с АО |
= 0,55 |
м, зонды АБК-3 и ИК-2; |
||||
г) малый градиент-зонд с АО |
= 0,55 |
м, зонд АБК-3 и большой гра |
||||
диент-зонд |
с АО = 4,25 м. |
Результаты интерпретации показаний |
||||
этих зондов на ЭВМ сопоставлены с результатами ручной обработки данных Б К З . Удельные сопротивления, полученные при интерпрета ции данных этих зондов на ЭВМ, хорошо согласуются с определе ниями пх по Б К З , если влияние зоны проникновения невелико.
Глава IX
ЛИТОЛОГИЧЕСКОЕ РАСЧЛЕНЕНИЕ РАЗРЕЗА
ИВЫДЕЛЕНИЕ КОЛЛЕКТОРОВ
22.П Е Т Р О Ф И З И Ч Е С К И Е ОСНОВЫ ЛИТОЛОГИЧЕСКОГО Р А С Ч Л Е Н Е Н И Я
ПОРОД I I В Ы Д Е Л Е Н И Я КОЛЛЕКТОРОВ
Расчленение разреза скважин с оценкой литологического харак тера пород и выделение коллекторов являются одними из важней ших" этапов интерпретации материалов каротажа. Получаемые при этом данные служат основой для дальнейшей геологической интер претации и выдачи окончательных заключений о характере пластов.
Литологическое расчленение разреза скважин и выделение кол лекторов в нем основаны на различии физических свойств горных пород [14, 20]. На рис. 66 приведены интервалы изменения физи ческих и геофизических свойств осадочных пород, составленные по данным работ [10, 14]. Видно, что отдельные литологические разно видности пород, в том числе и песчаники, характеризуются опреде ленными диапазонами физических и геофизических величин: ДС,
Д£/"пс, р п , А / т , V (или At) и А / ѵ ѵ (или плотность о), где AUnc, AJy и А / ѵ ѵ — приведенные существенные значения, вычисленные для кри-
160
вых ПС, гамма-каротажа (ГК) и гамма-гамма-каротажа (ГГК) по формуле (35). Но вместе с тем нн одна из рассматриваемых пород не может быть определена (опознана) по величине одного геофизиче ского параметра, так как интервалы этих параметров пересекаются. Поэтому однозначное определение породы требует получения гео физической характеристики по комплексу исследований.
Для получения геофизической характеристики горной породы (пласта) можно вместо численных значений геофизических величин
|
rVrjg |
|-J,7l7l- 7 |
') ') 'J |
в |
|
Рис. СО. Физические и геофизические свойства осадочных |
пород. |
|
1 — песчаник;' 2 — глинистый |
песчаник; 3 — глина; 4 — антрацит; 5 — каменный уголь; |
|
s — известняк; 7 — глинистый известняк; s — мергель; 9 — каменная соль; 10 — ангидрит; 11 — интервал изменения свойств.
пользоваться их приблизительными оценками. С этой целью диапазон значений каждой из геофизических величин ajt разбивают на отдельные градации (поддиапазоны) граничными значениями Ац, А,-2,..., Ajn, где ) —номер параметра; ?г —число градаций, которые обозначают двоичными кодами Kj [45 ]. Совокупность кодов К; отдельных гео физических величин (комплексный код К) дает геофизическую харак теристику породы.
При описании показанных на рис. 66 пород диапазоны значений геофизических величин целесообразно разбить на четыре градации (малые, средние, большие и очень большие) за исключением данных каверномера ДС и ПС, которые можно разбить на три градации
И Заказ Ш |
161 |
|
|
|
|
|
|
|
ï |
a б л к ц a 14 |
|
Диагностические |
коды для |
определения |
литологии пород |
|
|||||
|
|
|
|
|
Номер з и обозначение параметра и |
||||
|
|
|
|
|
|
кривой |
|
|
|
Индекс |
Название |
породы |
1 |
2 |
3 |
6 |
8 |
10 |
|
литологии I |
|||||||||
|
|
|
|
ДС |
ПС |
Рп |
г к |
АК |
ггк |
|
|
|
|
|
|
|
|
|
|
1 |
Песнанпк |
|
|
00 |
00 |
00 |
00 |
01 |
00 |
|
|
|
|
00 |
00 |
01 |
00 |
01 |
00 |
|
|
|
|
00 |
00 |
10 |
00 |
01 |
00 |
|
|
|
|
01 |
00 |
00 |
00 |
01 |
00 |
|
|
|
|
01 |
00 |
01 |
00 |
01 |
00 |
|
|
|
|
01 |
00 |
10 |
00 |
01 |
00 |
4 |
Глинистый |
песчаник |
00 |
01 |
00 |
01 |
01 |
00 |
|
|
|
|
|
00 |
01 |
01 |
01 |
01 |
00 |
|
|
|
|
01 |
01 |
00 |
01 |
01 |
00 |
|
|
|
|
01 |
01 |
01 |
01 |
01 |
00 |
|
|
|
|
00 |
01 |
00 |
10 |
01 |
00 |
|
|
|
|
00 |
01 |
01 |
10 |
01 |
00 |
|
|
|
|
01 |
01 |
00 |
10 |
01 |
00 |
|
|
|
|
01 |
01 |
01 |
10 |
01 |
00 |
2 |
Глина |
|
|
10 |
10 |
00 |
11 |
00 |
00 |
|
|
|
|
01 |
10 |
00 |
11 |
00 |
00 |
|
|
|
|
01 |
10 |
00 |
11 |
01 |
00 |
6' |
Антрацит |
|
|
01 |
— |
00 |
00 |
01 |
10 |
6 |
Каменный |
уголь |
01 |
— |
10 |
00 |
01 |
10 |
|
|
|
|
|
01 |
11 |
00 |
01 |
10 |
|
|
|
|
|
10 |
— |
10 |
00 |
01 |
10 |
|
|
|
|
10 |
— |
11 |
00 |
01 |
10 |
5 |
Известняк |
|
|
01 |
01 |
01 |
00 |
10 |
00 |
|
|
|
|
01 |
10 |
01 |
00 |
10 |
00 |
|
|
|
|
01 |
01 |
10 |
00 |
10 |
00 |
|
|
|
|
01 |
01 |
11 |
00 |
10 |
00 |
|
|
|
|
01 |
10 |
10 |
00 |
10 |
00 |
|
|
|
|
01 |
10 |
11 |
00 |
10 |
00 |
5' |
Известняк |
глинистый |
01 |
01 |
01 |
01 |
01 |
00 |
|
|
|
|
|
01 |
10 |
01 |
01 |
01 |
00 |
|
|
|
|
01 |
01 |
10 |
01 |
01 |
00 |
|
|
|
|
01 |
01 |
11 |
01 |
01 |
00 |
7 |
Мергель |
|
|
01 |
10 |
10 |
01 |
01 |
00 |
|
|
|
|
01 |
10 |
11 |
01 |
01 |
00 |
162
