Примененный при моделировании пластов алгоритм триангуляции заключается в следующем. Все точки (скважины) в рассматриваемой области соединяются между собой отрезками. Затем отрезки попарно проверяются на пересечение. Если отрезки пересекаются, то выбирается тот из них, который обеспечивает деление четырехугольника, образованного двумя парами вершин, на более равносторонние по указанному выше критерию, другой отрезок становится запрещенным и в дальнейшем сравнении не участвует. Такая выбраковка отрезков позволяет свести зависимость затрат времени Г на триангуляцию от числа точек К к квадратичной форме Г, К% а не факториальной 7\ К.
На рис. 9 приведена иллюстрация примера выбора одного из пересекающихся отрезков, который является ключевым моментом данного алгоритма триангуляции. Из отрезков АВ и CD выбирается АВ и бракуется CD> так как, исходя из критерия равно-сторонности, пара треугольников ABD и ABC вкупе лучше, чем пара треугольников CDA и CDB, хотя последний треугольник среди рассматриваемых четырех является наиболее равносторонним. В дальнейшем процессе выбраковки при сравнении со следующим пересекающимся отрезком может быть таким же образом запрещен и сам отрезок АВ.
Программа отдельно обрабатывает более сложные случаи выбраковки таких, как, например, при соосности отрезков или их параллельности осям прямоугольной области. Итогом триангуляции является полученный набор разрешенных отрезков, образующих оптимальную по выбранному критерию треугольную сетку, ячейками которой являются произвольные треугольники. Далее на прямоугольной области рассматриваемого пласта
158
Рис. 9. Пример выбора одного из пересекающихся отрезков
Рис 10. Пример интерполяции внутри треугольной ячейки
строится прямоугольная сетка с числом узлов по вертикали и по горизонтали соответственно 93 и 27. Те узлы, которые попадают в одну треугольную ячейку, соответствующую проекции наклонного треугольника в трехмерном пространстве на плоскость прямоугольной области, интерполируются этим треугольником.
На рис. 10 дан пример интерполяции внутри одной треугольной ячейки. Значения параметра пласта (например, отметки кровли) в точках 1-5, соответствующие узлам прямоугольной сетки 1'- 5' внутри ячейки А\В\Си берутся на пересечении с наклонным треугольником ABC. После определения во всех узлах прямоугольной сетки значений параметра пласта соответствующая сеточная функция считается построенной и записывается в базу данных. Таким образом строились все поверхности пластовых параметров.
Примеры объемных изображений построенных поверхностей показаны на рис. 11.
Для начальной поверхности ГВК подобрана плоскость с учетом наклона с юга на север от абсолютной отметки - 1187 (скв. 47, 12) до 1199 м (скв. 32, 45). На этом фоне отмечается изменение направления наклона поверхности ГВК на северо-восток [38].
При расчете начальных запасов по сеточным данным предлагается использовать формулу
a. . М.Л..Х
н ° х
(50)
где а - коэффициент газонасыщенности; i, / - узел разност-
159
Рис П. Примеры объемных изображений построенных поверхностей
ной сетки; Н - эффективная газонасыщенная толщина; М -коэффициент пористости; п - число узлов сетки; ри - начальное пластовое давление; 7о - стандартная температура; d - шаг сетки по оси х; d - шаг сетки по оси у; Zip У -
коэффициент сверхсжймаемости; ро - атмосферное давление.
Суммирование проводится для всех узлов разностной сетки внутри контура газоносности.
Проведя суммирование по выбранной группе узлов, можно получить значение начальных запасов любой зоны газовой залежи, например, одной из зон УКПГ или периферийной зоны.
При появлении промысловых данных по истории разработки становится возможным применение методов уточнения запасов с учетом фактической информации по падению пластового давления по скважинам.
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.