Оптимизация темпов разработки газовых месторождений с учетом ресурсосбережения: Производственно-практическое издание, страница 98

Примененный при моделировании пластов алгоритм триан­гуляции заключается в следующем. Все точки (скважины) в рассматриваемой области соединяются между собой отрезками. Затем отрезки попарно проверяются на пересечение. Если от­резки пересекаются, то выбирается тот из них, который обес­печивает деление четырехугольника, образованного двумя па­рами вершин, на более равносторонние по указанному выше критерию, другой отрезок становится запрещенным и в даль­нейшем сравнении не участвует. Такая выбраковка отрезков позволяет свести зависимость затрат времени Г на триангу­ляцию от числа точек К к квадратичной форме Г, К% а не факториальной 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   - шаг сетки по оси х; - шаг сетки по оси у; Zip У -

коэффициент сверхсжймаемости; ро - атмосферное давление.

Суммирование проводится для всех узлов разностной сетки внутри контура газоносности.

Проведя суммирование по выбранной группе узлов, можно получить значение начальных запасов любой зоны газовой за­лежи, например, одной из зон УКПГ или периферийной зоны.

При появлении промысловых данных по истории разработки становится возможным применение методов уточнения запасов с учетом фактической информации по падению пластового давле­ния по скважинам.