ADG3D

3. Данные, используемые программой.

   Исходные данные задачи содержат в себе информацию о геометрии модели, о физических характеристиках среды, о вариантах и условиях счета, а кроме того, некоторые паспортные сведения и информацию о скважинах. Элементам исходных данных даны обозначения. Они используются в названиях таблиц и в сообщениях справочной системы.

   Исходные данные состоят из следующих элементов

Скаляры.

Все эти параметры суть простые скалярные величины.

·    Решаемая задача - вариант (тип) решаемой задачи; Он может принимать восемь значений:

o        Нет - задача не выбрана;

o        Прямая - решается прямая задача;

o        Обратная - решается обратная задача и вслед за ней, по вычисленным значениям ro решается прямая задача.

o        Контактная - решается контактная задача и вслед за ней, по вычисленным значениям z решается прямая задача.

·    Количество:

o       слоев - количество слоев в разрезе (не более 5).

o       профилей - количество блоков по оси X (не более 100).

o       точек в профиле - количество блоков по оси Y (не более 100).

o       итераций - количество итераций для решения обратной или контактной задачи.

o       кривых - количество кривых наблюдения, по которым решаются обратная и контактная задачи.

·    Параметры модели

o       Расстояние между профилями - шаг между точками наблюдения по оси X.

o       Расстояние между точками - шаг между точками наблюдения по оси Y.

o       Радиус внешний -

o       Радиус внутренний -

o       Эти константы определяют размеры блоков в плане. Задаются в метрах, не должны быть более 55500 м.

o        

·    ly,lx- вокруг каждой точки наблюдения выделяется воображаемое прямоугольное окно : по ly блоков в обе стороны по оси Y и по lx блоков в обе стороны по оси X. Это окно используется для расчета поля в точке наблюдения и в нем же итерационно уточняются неизвестные плотности (или границы).

o       ly < ky , lx < kx.

Описание данных

Эта группа исходных данных не связана с основным вычислительным алгоритмом, а используется в оформлении. Эти данные не имеют обозначений.

·      Описание объекта - любая последовательность не свыше 60-ти символов.

·      Масштаб   - любая последовательность не свыше 11-ти символов.

Входные массивы.

Следующие ниже данные задаются в виде таблиц размером ку*kx:

·      z - таблицы границ. Отметки задаются в метрах. Отметки ниже уровня моря задаются отрицательными числами. Наибольшая отметка : 55500 м.Задается ks+1 таблица по числу границ. Первая граница - поверхность наблюдения, последняя граница - подошва нижнего слоя. z задается для прямой и обратной задач.

·      z1 - начальные значения внутренних границ. Задаются аналогично z, для контактной и смешанной задач.

·      sz - начальные погрешности оценки границ. Задаются для контактной и смешанной задач.

·      ro - таблицы плотностей. ks таблиц, по числу слоев. Задаются избыточные или абсолютные значения плотности. Для прямой и контактной задач.

·      ro1 - начальные значения плотностей. Задаются аналогично ro. Для обратной и смешанной задач.

·      sr - начальные погрешности оценки плотностей. Задаются для обратной и смешанной задач.

·      dg - значения наблюдаемой аномалии силы тяжести. Задаются в мГал, для обратной и контактной задач.

·      s0 - среднеквадратическая погрешность значений dg.

·    u - уровень (в километрах) для каждой кривой наблюдения, на который будут пересчитаны значения поля dg, от модели, подобранной в результате решения обратной или контактной задач. Если u не задано, то поле рассчитывается для поверхности наблюдения. Отметки ниже уровня моря задаются отрицательными числами.

·    fij - фоновые составляющие плотности. Они задаются для каждого слоя, если предполагается вводить и получать абсолютные, а не избыточные значения плотностей.

·    b - вектор из ks констант ( задается в м/сек на единицу изменения плотности ).

·    v - вектор из ks констант ( задается в м/сек ).

 Для вычисления статических поправок, нулевых времен и средних скоростей используются коэффициенты bk и vk, связывающие избыточные плотности ro и скорости V в слое k и блоке i,j по формуле: V[k,i,j] = b[k]*ro[k,i,j]+v[k].

Данные по скважинам сведены в таблицу. Считается, что

 каждая скважина привязана к определенной точке наблюдения на

 поверхности и имеет свой регистрационный номер. В скважине

 может быть до 6 пластов, каждый из которых описывается отмет-

 кой нижней границы и значением плотности. Границы и плотнос-

 ти для скважин задаются по тем же правилам, что и для слоев.

 В эту же таблицу заносятся данные по шурфам и пикетам на пло-

 щади. Эти объекты характеризуются только собственным номером

 и координатами точки привязки. Максимальное суммарное число

 объектов в таблице скважин - 24.

Выходные массивы.

Входные массивы отделяются от выходных массивов буферным массивом. В списке он обозначается ‘----------’ и может использоваться для хранения результатов пересылки массивов.

·      zx - таблицы границ, результат решения контактной задачи.

·      rox - таблицы плотностей, результат решения обратной задачи.

·      nevz -  невязки - вектор среднеквадратических невязок (в мГал) между фактическими значениями dg и вычисленными по моделям, формировавшимся в итерациях при решении обратной или контактной задачи. Он содержит ki значений

·      dg - значения поля dg, полученные при решении прямой задачи, а также полученные при решении прямой задачи по модели, полученной при решении обратной или контактной задачи.

·      tau - вектор, характеризующий задержки продольной волны, проходящей от поверхности наблюдения до подошвы последнего слоя, вызванные изменением плотности. Другими словами - cтатические поправки за пункт приема отраженных волн.

·      t0 - нулевые времена прихода волн, отраженных от нижней границы.

·      vср - средние скорости прохождения волн от поверхности до подошвы