Математические модели местности

Математические связи между исходными точками цифровых моделей описываются линейными либо нелинейными (степенными) зависимостями. В первом случае связь между смежными точками модели описывается уравнениями плоскостей, проходящими через каждые три смежные точки модели, во втором — криволинейными поверхностями разного порядка, и, таким образом, рельеф местности задается либо множеством пересекающихся между собой плоскостей, либо поверхностей различной кривизны.

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

Подавляющее число регулярных и нерегулярных ЦММ предполагают при последующем математическом моделировании линейную интерполяцию высот (глубин) между смежными точками модели.

Задача определения высот точек, уровней воды и поверхностей геологических напластований сводится к нахождению в каждом случае тех трех смежных исходных точек модели, между которыми попадет соответствующая искомая точка, в нахождении коэффициентов уравнения плоскости, проходящей через эти три точки, и, наконец, в определении по полученному уравнению искомой высоты (рис. 10.2).

Рис. 10.2. Линейное математическое моделирование рельефа и инженерно-геологического строения местности.

Если искомая точка (например, 20 — рис. 10.2) попадает между смежными исходными точками ЦММ с номерами /, к и /, то уравнение искомой плоскости в общем виде может быть представлено:

В уравнении (10.3.1) известны проектные координаты X и У искомой точки (например, 20 — рис. 10.2), высоту которой нужно определить, но не известны коэффициенты А, В и С уравнения плоскости, проходящей через исходные точки /', к и / цифровой модели.

Если в уравнение (10.3.1) подставить известные координаты трех исходных точек цифровой модели, то получим три уравнения, в которых не известны только три коэффициента — А, В

и С:

Система уравнений (10.3.2) решается в матричной форме или методом «прогонки», в результате чего определяют неизвестные коэффициенты А, В и С уравнения (10.3.1), подставив в которое проектные координаты X и У искомой точки, определяют ее высоту Н.

Наиболее универсальными являются статистические ЦММ (10.1.6), математическая реализация которых заключается в использовании метода «плавающего квадрата» или «плавающего круга», в пределах которого строится криволинейная поверхность п — го порядка (рис. 10.3).

Рис. 10.3. Математическое моделирование рельефа «плавающей» криволинейной поверхностью: 1 — точки статистической ЦММ; 2 — искомые точки; 3 — проектная трасса линейного сооружения.

Наиболее часто для математического моделирования рельефа используются уравнения поверхности 2-го порядка:

где X, У — известные проектные координаты точки, высоту которой требуется определить; А, В, С, D, Е, F — коэффициенты уравнения аппроксимирующей поверхности 2-го порядка.

Основная идея «плавающей» аппроксимации заключается в том, что по аппроксимируемой поверхности от точки к точке перемещается круг или квадрат таким образом, что каждая точка поверхности, высоту которой требуется определить, размещается в его центре (например, точка 21 на рис. 10.3). Радиус круга или размеры стороны квадрата автоматически устанавливаются такими, чтобы в их пределы попало не менее десяти исходных точек модели. Поскольку радиус круга или размеры сторон квадрата меняются с дискретным шагом соответственно Дг и ДЬ, то в пределах выделяемых ими площадей может оказаться и более десяти точек модели (например, 11,12,13 и т. д.).

Поскольку коэффициенты А, В, С, D, Е, F в аппроксимирующем уравнении (10.3.3) не известны, то для кахздой точки модели, попавшей в пределы круга или квадрата, записываются уравнения:

где А, В, С, D, Е, F — неизвестные коэффициенты уравнения аппроксимирующей поверхности; H_jr Ху, У),…, Н", Х", У_п — известные координаты точек модели, попавших в пределы круга или квадрата; п — число точек ЦММ, попавших в пределы круга или квадрата.

Поскольку число неизвестных в системе (10.3.4) меньше числа уравнений (которых не менее 10), то система решается методом «наименьших квадратов». Таким образом, определяют неизвестные коэффициенты аппроксимирующего уравнения (10.3.3), подставив в которое известные проектные координаты X и Уточки (например, 21), определяют ее высоту или глубину Н.

Далее круг или квадрат перемещают в центр очередной искомой точки, и процедура повторяется. При этом если плотность исходных точек модели в районе очередной определяемой точки уменьшилась, то размеры круга или квадрата автоматически возрастут, а если плотность возросла — то наоборот уменьшатся.

Для математического описания ситуационных, почвенногрунтовых, гидрогеологических и других условий местности используют контурную индексацию объектов местности с перечнем номеров точек вдоль каждого такого контура (урез воды, подошва подводной части пляжа, газопровод и т. д.). Для замкнутых контуров (здания, причалы, волноломы и т. д.) точки замыкания повторяются.