СПОСОБ КОЛИЧЕСТВЕННОГО РАСЧЕТА НАСЫЩЕННОСТИ ТРЕЩИННОГО КОЛЛЕКТОРА УГЛЕВОДОРОДАМИ

1. Способ количественного расчета насыщенности трещинного коллектора углеводородами, отличающийся тем, что включает следующие стадии:
С1) на основе данных керна полного диаметра определяют зависимость удельного сопротивления (I) от водонасыщенности (Sw) матрицы; на основе анализа данных герметизированого керна определяют реальную насыщенность водой (нефтью/газом) исходной формации трещинного коллектора;
С2) определяют трещинную пористость на различных глубинах трещинного коллектора с использованием графического отображения известных данных каротажного зондирования; рассчитывают показатели удельного сопротивления на различных глубинах с использованием данных двойного бокового каротажного зондирования в сочетании с данными анализа керна;
С3) создают модель перколяционной сетки, включающей матрицу и трещину, в соответствии с известными региональными особенностями структуры пор, пористостью и проницаемостью исследованной области, и определяют зависимость между показателем удельного сопротивления (I) и водонасыщенностью (Sw) трещинного коллектора посредством дренажного моделирования;
С4) осуществляют калибровку результатов численного моделирования в соответствии с моделью перколяционной сетки на основе данных анализа керна и результатов анализа герметизированого керна;
С5) после калибровки определяют вариации зависимости между показателем удельного сопротивления (I) и водонасыщенностью (Sw) при различной трещинной пористости посредством численного моделирования;
С6) определяют функцию зависимости между показателем удельного сопротивления и водонасыщенностью при различной трещинной пористости посредством интерполяции;
С7) рассчитывают насыщенность трещинного коллектора углеводородами в соответствии с пористостью трещин, определенной на стадии С2, посредством выбора интерполяционной модели на основе функций зависимости между показателем удельного сопротивления и водонасыщенностью, определенных на стадии С6.
2. Способ по п.1, отличающийся тем, что показатели удельного сопротивления на различных глубинах рассчитывают следующим образом: удельное сопротивление нетронутой формации Rt определяют по данным глубинного бокового каротажного зондирования, а удельное сопротивление водонасыщенной формации R0 рассчитывают по формуле Арчи на основе данных каротажного зондирования пористости в сочетании с данными анализа керна.
3. Способ по п.1, отличающийся тем, что в качестве используемых на стадии С2 данных каротажного зондирования используют отображения данных полнопроходного сканирования микросопротивления или электроакустического каротажа.
4. Способ по п.1, отличающийся тем, что модель перколяционной сетки, включающей матрицу и трещину, на стадии С3 создают посредством определения радиусов и распределения пор и определения радиусов и распределения каналов, определения соединений между порами и создания трехмерной модели структуры пор.
5. Способ по п.1, отличающийся тем, что региональные особенности структуры пор на стадии С3 могут быть определены методами ртутной порометрии, ядерного магнитного резонанса и электронной микроскопии.
6. Способ по п.1, отличающийся тем, что при создании модели на стадии С3 используемые в модели параметры, главным образом радиусы и распределение пор в матрице, пористость, проницаемость и аналогичные параметры, соответствуют результатам анализа керна.
7. Способ по п.1, отличающийся тем, что результаты численного моделирования на стадии С4 калибруют следующим образом:
результаты численного моделирования при трещинной пористости, равной 0, полученные на стадии С3, калибруют на основе результатов измерения электрических свойств керна (зависимость I-Sw) в заданном интервале, в частности результаты численного моделирования при трещинной пористости, равной 0, приводят в соответствие с результатами измерения электрических свойств керна посредством подгонки параметров численного моделирования;
затем результаты моделирования, полученные на стадии С3, калибруют с использованием значений водонасыщенности, определенных по результатам анализа герметизированого керна, показателя удельного сопротивления на соответствующей глубине, рассчитанного на стадии С2, и трещинной пористости на соответствующей глубине, рассчитанной на стадии С1; в частности результаты численного моделирования при трещинной пористости, соответствующей указанной глубине, приводят в соответствие с результатами анализа герметизированого керна посредством подгонки параметров численного моделирования.
8. Способ по любому из пп.1-7, отличающийся тем, что параметры численного моделирования включают параметры структуры пор и параметры свойств флюида.
9. Способ по п.1, отличающийся тем, что значения пористости на стадии С5 находятся в пределах от 0,1% до 0,5%.
10. Способ по п.1, отличающийся тем, что оптимальная интерполяционная функция, определяемая на стадии С6, имеет вид:
<mrow><mi>I</mi><mo>=</mo><mi>a</mi><msup><mi>e</mi><mrow><mi>b</mi><msubsup><mi>S</mi><mi>w</mi><mrow><mo>−</mo><mi>c</mi></mrow></msubsup></mrow></msup><mo>,</mo></mrow>
w
a
где I - показатель удельного сопротивления; S - водонасыщенность; , b и c - константы,
a
при матричной пористости 3% и трещинной пористости 0,0%: =0,03, b=3,44, c=0,52;
a
при матричной пористости 3% и трещинной пористости 0,1%: =0,06, b=2,70, c=0,52;
a
при матричной пористости 3% и трещинной пористости 0,2%: =0,20, b=1,47, c=0,67;
a
при матричной пористости 3% и трещинной пористости 0,3%: =0,38, b=0,85, c=0,77;
при матричной пористости 3% и трещинной пористости 0,4%: a=0,52, b=0,53, c=0,91.