Обработка результатов измерений гравитационного поля по участку
Дипломная работа - Геодезия и Геология
Другие дипломы по предмету Геодезия и Геология
медианный фильтр (середина чисел, расп. в порядке возрастания)
Для расчетов я выбрал m=2 для того чтобы не потерять в точности полученных данных и получить большее количество значений.
результат фильтрации по каждому профилю:
Таким образом, с помощью применения медианного фильтра и фильтра аппроксимации мы выделили локальную составляющую поля из значений наблюденного поля. Фильтр аппроксимации позволяет выделить региональную составляющую сигнала. Вычитая ее из значений наблюденного поля, мы получили сигнал локальной составляющей и помехи. Сигнал помехи мы исключили медианным фильтром.
Графики автокорреляционных функций
Для изучения корреляционных свойств сигналов используется автокорреляционная функция:
j = 0,1,2,…, m, (2)
, (3)
Процедура вычисления АКФ
По графикам АКФ видим, что прослеживается закономерность в значениях аномальной составляющей.
По графикам АКФ случайной составляющей определяем, что перед нами действительно помеха, так как все графики распределены по нормальному закону с нулевым средним
Выделение объектов с различной плотностью
Выделение объектов с различной плотностью производится с помощью программного модуля prizma-g.exe, куда заносятся значения отфильтрованной аномальной составляющей гравитационного поля по каждому профилю в отдельности. После ввода значений получаем график измеренных значений приращения силы тяжести. Внизу графика поле для построения объектов различной палитры допустимых цветов, которые отличаются избыточной плотностью. Значение избыточной плотности мы вводим сами исходя из собственных соображений и имеющегося опыта.
После вычислений на экран выводится графики измеренного (сплошная линия) и вычисленного (пунктир) поля. Сравнив эти два графика, принимаем решение относительно дальнейших действий. Можем произвести корректировку разреза и повторить вычисления, если результаты моделирования сильно отличаются от полевых данных. Если поле над моделью удовлетворительным образом совпадает с результатами полевых измерений, принимается решение о прекращении процесса подбора.
Примеры интерпретации:
Профиль 134
Профиль 128
На основе полученных данных об объектах с избыточной плотностью по каждому профилю, строим модель исследуемого участка.
Оценка плотности горных пород на площади исследования
Модель представлена объектами с плотностями =2.53, =2.59, =2.69, =2.73.если это класс магматических г.п., то это кислые горные породы, например, гранит; из минералов это может быть: ортоклаз (=2.53), альбит (=2.60), кварц (=2.67), кальцит (=2.70).
Оценивая по статистическим параметрам, данная группа модели подходит под известняки, граниты и порфириты.
Физико-геологическая модель исследуемого участка
регрессия автокорреляционный плотность порода
Исходные данные:
Значения ?g в каждой точке измерения (пикете) по всей системе профилей
Профиль138Профиль136Профиль134Профиль132Профиль130Пикет?gПикет?gПикет?gПикет?gПикет?g60-16.6660-16.7360-16.8060-16.7160-16.4862-16.7662-16.6262-16.8062-16.6562-16.3664-16.6964-16.4964-16.6864-16.5664-16.2366-16.6566-16.4766-16.5166-16.4666-16.2068-16.6268-16.4468-16.4168-16.3668-16.1370-16.5670-16.3970-16.2770-16.2170-15.9772-16.5172-16.3272-16.1172-16.0572-15.9074-16.4074-16.1074-15.9174-15.8474-15.8176-16.2676-15.8976-15.8676-15.9176-15.6178-16.1978-15.8378-15.7578-15.8478-15.6380-16.4980-15.6580-15.6080-15.7980-15.5682-16.2782-15.6982-15.5082-15.7282-15.5284-16.1584-15.7384-15.4684-15.5084-15.6686-15.6586-15.7486-15.4386-15.4586-15.3888-15.3788-15.8288-15.4388-15.3288-15.2990-15.2590-15.4890-15.5190-15.2390-15.1392-15.1792-15.0992-15.6192-15.2192-15.0394-14.9694-14.8494-15.3094-15.1094-14.9796-14.8196-14.6596-15.2796-15.0996-14.9498-14.7098-14.6598-15.1498-14.9498-14.84100-14.55100-14.63100-14.94100-14.94100-14.68102-14.44102-14.54102-14.54102-14.93102-14.55104-14.37104-14.47104-14.51104-14.82104-14.64106-14.32106-14.36106-14.39106-14.51106-14.56108-14.18108-14.35108-14.27108-14.41108-14.23110-14.13110-14.19110-14.07110-14.23110-14.20112-13.66112-14.14112-13.97112-14.09112-14.12114-13.47114-13.88114-13.84114-14.07114-14.10116-13.36116-13.70116-13.77116-13.93116-14.08118-13.35118-13.60118-14.03118-13.88118-13.87120-13.36120-13.68120-13.97120-13.97120-13.84122-13.30122-13.43122-13.62122-14.15122-13.83124-13.24124-13.34124-13.43124-13.90124-13.86126-13.17126-13.21126-13.26126-13.66126-14.02128-13.12128-13.12128-13.20128-13.51128-14.33130-13.18130-13.23130-13.03130-13.46130-14.09132-13.08132-13.11132-12.98132-13.30132-13.77134-12.96134-13.09134-12.94134-13.20134-13.58136-12.70136-12.90136-13.00136-13.02136-13.21138-12.79138-12.81138-12.87138-12.95138-12.79140-12.61140-12.71140-12.73140-12.89140-12.68142-12.61142-12.56142-12.65142-12.92142-12.66144-12.50144-12.56144-12.76144-12.69144-12.68146-12.51146-12.51146-12.65146-12.56146-12.68148-12.40148-12.50148-12.55148-12.36148-12.71150-12.34150-12.46150-12.43150-12.34150-12.49152-12.25152-12.35152-12.38152-12.06152-12.25154-12.06154-12.25154-12.25154-12.18154-12.32156-12.01156-12.10156-12.17156-12.21156-12.20158-12.14158-12.13158-12.04158-12.21158-12.15160-12.03160-12.12160-12.09160-12.12160-12.03162-11.99162-12.01162-12.11162-12.13162-11.98164-11.93164-11.94164-12.16164-12.27164-11.93166-11.80166-11.87166-11.99166-11.94166-11.79168-11.85168-11.80168-11.98168-11.84168-11.74170-11.69170-11.80170-11.88170-11.78170-11.61172-11.63172-11.68172-11.83172-11.68172-11.54174-11.47174-11.74174-11.82174-11.70174-11.55176-11.34176-11.52176-11.72176-11.60176-11.53178-11.25178-11.41178-11.60178-11.52178-11.52180-11.09180-11.35180-11.52180-11.62180-11.46182-11.05182-11.22182-11.45182-11.57182-11.41184-11.00184-11.17184-11.51184-11.44184-11.39186-10.89186-11.07186-11.37186-11.33186-11.24188-10.81188-10.98188-11.18188-11.17188-11.17190-10.63190-10.82190-11.06190-11.12190-11.04192-10.47192-10.71192-10.98192-11.04192-11.02194-10.36194-10.71194-10.90194-10.96194-10.77196-10.13196-10.65196-10.79196-10.85196-10.64198-9.99198-10.42198-10.70198-10.65198-10.50200-9.94200-10.00200-10.52200-10.48200-10.47202-9.97202-10.10202-10.08202-10.28202-10.16204-9.97204-9.98204-10.01204-10.19204-10.00Профиль128Профиль126Профиль124Профиль122Профиль120Пикет?gПикет?gПикет?gПикет?gПикет?g60-16.1460-15.8160-15.6160-15.3960-15.1862-16.1962-15.8362-15.7562-15.5362-15.0564-16.1864-15.8164-15.8164-15.5864-15.1166-15.9266-15.6566-15.8666-15.7766-15.1868-15.8068-15.6568-15.8568-16.0268-15.2770-15.7070-15.6470-15.8870-15.9970-15.3672-15.5972-15.6772-15.7872-15.8972-15.1374-15.5574-15.7174-15.5474-15.5574-14.7376-15.6176-15.4676-15.3876-15.3676-14.8378-15.6178-15.3278-15.3078-15.3478-15.0380-15.4980-15.3380-15.0280-15.2480-15.0382-15.3882-15.2482-14.9682-15.2182-14.9884-15.3384-15.1784-14.8684-15.0484-14.9286-15.2686-15.0386-14.8686-14.8686-14.9188-15.1288-14.6788-14.8088-14.7588-14.8990-14.8790-14.6590-14.7390-14.5190-14.7892-14.7392-14.6592-14.6892-14.4492-14.7394-14.