Численно-аналитическая модель теплового поля поверхностных слоев земли

Гаврилов В.В. Статья в формате PDF 113 KB Необходимость постановки геотермических исследований определяется практическими задачами, связанными с проблемами поисков, разведки и разработки глубокозалегающих месторождений полезных ископаемых, а также использования внутриземного тепла.

Распределение теплового поля в сечении описывается двумерным уравнением стационарной теплопроводности в неоднородной среде в прямоугольной области

где T( x, y) - температура, λ( x, y) - коэффициент удельной теплопроводности, Q -мощность внутренних источников и стоков тепла (теплогенерация пород). Для задания граничных условий используется средняя температура по  верхности Земли T0 и тепловой поток qm x) на  границе литосферы и астеносферы (зоны частичного плавления пород).

Для определения функции T( x, y) распределения λ( x, y) и Q( x, y) должны быть известными. Для получения этих распределений используются опубликованные данные для разных типов пород региона, а также геологическая интерпретация сейсмических исследований.

Уравнение (1) преобразуем, получив линейное дифференциальное уравнение второго порядка эллиптического типа. Итак, задача решается конечно-разностным методом. Для этого строится сетка с постоянным или переменным шагом.

В работе для итерационных расчетов используются метод Гаусса-Зейделя и метод релаксаций. При чем Метод Гаусса-Зейделя применительно к эллиптическим разностным уравнениям называется методом Либмана или методом последовательных смещений.

Таким образом, в данной работе с помощью численного эксперимента и сравнения с точным решением удалось установить, что предложенная разностная схема при разрывной функции λ(x,y) обладает первым порядком точности. Применение обобщенного правила Рунге для оценки погрешности численного решения показало хорошее соответствие полученных оценок фактическим погрешностям. Применение повторной экстраполяции позволило на несколько порядков уточнить полученные численные данные.

Таблица 1. Сравнение численного результата с аналитическим

y

T yi )

T ( y)

ε, %

0

5

5

-

1000

19.55

19.39

0.8

2000

33.93

33.77

0.5

3000

48.30

48.13

0.3

4000

62.63

62.45

0.3

Программный продукт реализован в среде Microsoft Visual C++ версии 5.0 (под Windows), и зарегистрирован в РОСПАТЕНТЕ. Кроме того получены справки о практическом и научном использовании данной работы.






В ОПЕК думают о Калуге

25 ноября в Калуге многое изменится. Рассматривая вопрос о вступлении России в ОПЕК, страны картеля решили открыть свое представительство именно в Калуге, на Калугинском бульваре. Анатолий Васильевич Квашнин, как официальный представитель Администрации Президента, высоко оценил этот шаг, а Вячеслав Краснов , будущий глава представительства, пообещал как можно более выгодные результаты для обеих сторон.