Способ контроля лесопожарной опасностиПатент на изобретение №: 2147253 Автор: Давыдов В.Ф., Новоселов О.Н., Харин О.А., Азметов Р.Р., Мещерякова И.А., Щербаков А.А. Патентообладатель: Московский государственный университет леса Дата публикации: 10 Апреля, 2000 Начало действия патента: 28 Мая, 1999 Адрес для переписки: 141005, Московская обл., Мытищи, МГУЛ, Патентный отдел ИзображенияИспользование: лесное хозяйство, в дистанционной оперативной оценке состояния лесных горючих материалов (ЛГМ) на обширных площадях. Сущность изобретения: при контроле лесопожарной опасности регистрируют собственное излучение подстилающей поверхности, преобразуют зарегистрированную функцию электрического сигнала в цифровые матрицы отсчетов зависимости амплитуды от координат, калибруют тракт зондирования по измерениям эталонных участков, производят вычисление характеристик матриц и расчет результирующего показателя. При этом собственное излучение регистрируют в диапазоне, соответствующем максимуму теплового ИК-излучения подстилающей поверхности. В качестве результирующего показателя используют величину влажности ЛГМ. Величину влажности ЛГМ внутри выделенных контуров рассчитывают по регрессионной зависимости. Перечисленные признаки обеспечивают оперативное, дистанционное отслеживание влажности лесных горючих материалов. 5 ил., 3 табл. , , , , , , ОПИСАНИЕ ИЗОБРЕТЕНИЯ К ПАТЕНТУИзобретение относится к лесному хозяйству, в частности к дистанционной оперативной оценке состояния лесных горючих материалов (ЛГМ) на обширных площадях. Лесные пожары были и остаются доминирующим фактором, определяющим структуру и динамику бореальных лесов. По масштабу воздействия на лесной фонд они существенно превышают все остальные вместе взятые факторы. Регулярное слежение за лесопожарной обстановкой в лесах, осуществляемое ограниченными авиационными и наземными средствами, не обеспечивает требуемой периодичности и достоверности лесопожарного мониторинга. Известен способ пирологической оценки лесов по комплексному показателю (КПО) горимости В.Г. Нестерова (см., например, "Лесопожарные показатели засухи", в сборнике "Пирологическое районирование в таежной зоне", Софронов М. А. , Волокитина А.В., АН СССР, Сибирское отделение. Наука, Новосибирск, 1990 г., стр. 37...38 - аналог). В способе-аналоге расчет КПО проводится по каждому району, области, лесхозу, лесничеству на основе данных наземных метеостанций и метеопунктов по следующей зависимости: где To - температура воздуха на 12 часов местного времени, oC; tp - температура точки росы, oC; - сумма разности температур за все дни i сухого периода со дня схода снежного покрова до дня выпадения осадков в количестве более 3 мм. После выпадения осадков в количестве более 3 мм счет обнуляется, а расчет КПО начинается заново со дня установления бездождевой погоды. По полученной сумме КПО устанавливаются классы лесопожарной опасности по В.Г. Нестерову: I класс до 300, II класс 300...1000, III класс 1000...4000, IV класс 4000...8000, V класс чрезвычайной опасности - более 8000. Выделенные участки различных классов наносят на пожарную карту лесхоза (лесничества) по укрупненным выделам и закрашивают красным цветом различной насыщенности. Недостатками известного аналога являются: - существенные ошибки результирующих оценок при неравномерном выпадении осадков на площади наблюдении; - невысокая достоверность при редкой сети метеопунктов в таежной зоне; - неадекватность (косвенность) самого показателя, характеризующего засушливое состояние погоды, а не состояние горючих материалов. Восприимчивость лесного отпада и подстилки к огню напрямую зависит от их влажности. Самым надежным способом определения пожарной опасности являлся бы способ непосредственного измерения влажности самих ЛГМ. Поэтому применительно к задачам оценки и прогнозирования лесопожарной опасности наиболее важной проблемой является повышение оперативности и точности определения влагосодержания ЛГМ. Наземные методы непосредственного расчета влагосодержания ЛГМ путем взятия множества точечных проб неприемлемы из-за большой трудоемкости и длительности анализа. В большинстве национальных систем оценки и прогноза лесопожарной опасности наряду с наземными измерениями используют данные дистанционного зондирования подстилающей поверхности (см., например, Г.Н. Коровин и др. "О возможности регистрации полей осадков дистанционными методами", в научном сборнике "Аэрокосмические методы и геоинформационные системы в лесоведении и лесном хозяйстве" (ЦЭПЛ РАН, Москва, издательство МГУЛ, 1998 г., стр. 187-190). Ближайшим аналогом по технической сущности с заявляемым техническим решением является дистанционный способ определения пожарной опасности в лесу по а.с. СССР N 1648505, кл. A 62 C 3/02, опубл. 15.09.1991). В способе ближайшего аналога осуществляют регистрацию собственного излучения подстилающей поверхности в диапазоне, соответствующем максимуму теплового инфракрасного (ИК) излучения подстилающей поверхности, калибровку тракта зондирования по измерениям эталонных участков, вычисление числовых характеристик и расчет результирующего показателя. Недостатками ближайшего аналога являются: - невозможность непосредственного применения для измерения влажности из-за разнородности измеряемых величин (в известном решении это произведение радиационной и радиояркостной температур, а в заявляемом способе - влажность); - неадекватность результирующего показателя в виде произведения радиояркостной и радиационной температур физическому процессу. Задача, решаемая заявляемым способом, заключается в обеспечении оперативного, дистанционного отслеживания влажности лесных горючих материалов на основе получения изображения лесов в инфракрасном диапазоне электромагнитного спектра. Поставленная задача решается тем, что в способе контроля лесопожарной опасности, включающем регистрацию собственного излучения подстилающей поверхности в диапазоне, соответствующем максимуму теплового ИК-излучения подстилающей поверхности, калибровку тракта зондирования по измерениям эталонных участков, вычисление числовых характеристик и расчет результирующего показателя, в качестве результирующего показателя используют непосредственно величину влажности W лесных горючих материалов, преобразуют зарегистрированную функцию электрического сигнала в цифровые матрицы отсчетов зависимости амплитуды от координат, выделяют методами пространственного дифференцирования контуры на изображениях при установленных пороговых значениях градиента, а величину влажности лесных горючих материалов внутри выделенных контуров рассчитывают по регрессионной зависимости: где a - поправочный коэффициент, учитывающий параметры тракта зондирования, географическую зону и тип лесов; 2 - мощность переменной составляющей сигнала участка изображения внутри контура; мощность постоянной составляющей сигнала участка изображения внутри контура. Изобретение поясняется чертежами, где: фиг. 1 - график численного решения для функции влажности W (КПО); фиг. 2 - иллюстрация зависимости изменения параметров ИК-сигнала от возрастания КПО; фиг. 3 - результаты программной обработки ИК-изображения лесного региона, полученного со спутника NOAA (США); фиг. 4 - функциональная схема устройства; фиг. 5 - график функции регрессии. Вновь введенные операции, образующие совокупность существенных признаков, обеспечивают достижения таких качественных свойств способа, как: - адекватность контроля физического процесса, когда контролируется состояние самих ЛГМ, а не воздушной среды; - большой пространственный размах контроля при одновременной возможности точной привязки результатов по координатам; - достоверность и статистическая устойчивость результатов за счет использования эффективной мощности регистрируемого сигнала при расчете интегрального показателя. Техническая сущность изобретения заключается в следующем. В соответствии с законом Стефана-Больцмана интегральная светимость тел L(T) пропорциональна четвертой степени температуры ToK и зависит от их излучательной способности (T) (см. , например, Левитин И.Б. "Инфракрасная техника", Энергия, Ленинградское отделение, 1973 г., стр. 15). L(T) = (T)0T4, (Вт/м2), где 0 - постоянная, равная 5,6710-8 Вт/м2K-4. Максимум мощности излучения приходится на длину волны max, определяемую по закону Вина: max = b/To, b = 2897,8 мкмK; Следовательно, при пожароопасной температуре воздуха от 15 до 35oC, максимальная мощность собственного ИК-излучения лесов приходится на диапазон 8. . .10 мкм. На этот же диапазон приходятся окна прозрачности атмосферы при коэффициенте прохождения ИК-излучения до 70%. Принимаемое ИК-излучение элементов ландшафта (растительности, почв, водной поверхности) образуется суммированием излучения от элементарных площадок, попадающих в элемент разрешения сканирующих устройств зондирования. Разрешающая способность космических радиометров в ИК-диапазоне, таких как МСУ-СК (Россия), AVHRR (спутник NOAA, США) составляет от 0,6 до 1,1 км. Следует отметить, что совокупность участков, относящихся к одному типу леса, всегда отличается значительной пирологической неоднородностью за счет возрастных различий в древостое, просветов в пологе, чередования возвышенности и низин. Радиационная температура различных материалов под пологом и на открытых участках, измеренная в натурных условиях в период с 12 до 14 часов местного времени, существенно различается (см., например, "Лесная пирология", учебное пособие под редакцией Л.Б. Калинина, из-во СПб. Лесотехнической академии, 1993 г., стр. 23-27). Радиационная температура лесных материалов показана в табл.1. Основные проводники горения относятся к классу капиллярно-пористых тел, обладающих высокой гигроскопичностью. Лишайники, мхи, подстилка на открытых участках получают большой приток радиации и быстро подсыхают. Каждому сочетанию температуры воздуха соответствует определенная влажность ЛГМ, называемая равновесной. Последняя является тем пределом, к которому стремится влажность ЛГМ при их высыхания и увлажнения в воздушной среде. Реальная влажность ЛГМ находится в интервале от 7 до 80% (см., например, Софронов М. А. , Волокитина А.В. "Пирологическое районирование в таежной зоне", АН СССР, Сибирское отделение, Наука, Новосибирск, 1990 г., стр. 106-117, "Увлажнение и высыхание ЛГМ"). В табл. 2 представлены экспериментальные данные зависимости между комплексным показателем В.Г. Нестерова, влажностью основных ЛГМ и их горимостью (см., например, аналог, стр. 109). Чем больше разность между текущей влажностью (W) ЛГМ и их равновесной влажностью (Wp), тем больше скорость процесса высыхания (см. например, аналог, стр. 34) Из решения дифференциального уравнения высыхания получено: Заменяя независимую переменную из выражения для КПО как и используя в качестве начальных условий для решения дифференциального уравнения данные табл. 2, получили численное решение для функции влажности W (КПО), представленное графиком фиг. 1. В интервале пожароопасных значений влажности ЛГМ (80...10%), среднестатистическая зависимость W (КПО) при начальных условиях табл. 2 аппроксимируется выражением: где Wo 1,2; Wp 0,1; 0 500. В зависимости от влажности ЛГМ изменяется и мощность восходящего ИК-излучения лесных участков. Влажные ЛГМ экранируют ИК-излучение земного покрова, в результате кажущаяся радиационная температура увлажненных участков оказывается ниже истинной температуры объектов. Общее аналитическое выражение для мощности сигнала от элемента разрешения радиометра представляется в виде: L(To) = 1/So i(Ti)0T4iSi, где So - площадь элемента разрешения радиометра, м2; Si - площади разнородных участков внутри элемента разрешения, имеющих каждый различную влажность (Wi), излучательную способность i(Ti) и локальную температуру Ti. При большой влажности перепад температур между открытыми участками и в тени достигает десятка градусов, в результате наблюдаются большие флюктуации амплитуды сигнала от одного элемента разрешения к другому. По мере роста КПО влажность уменьшается, увеличивается средняя радиояркостная температура подстилающей поверхности и уменьшается размах флюктуаций амплитуды сигнала между соседними элементами разрешения. Таким образом, селектируемыми признаками нарастания лесопожарной опасности по параметрам регистрируемого ИК-радиометром сигнала являются две характеристики: - рост постоянной составляющей (тренд математического ожидания - уменьшение флюктуаций амплитуды сигнала (уменьшение среднеквадратического отклонения ). Количественной мерой скорости флюктуаций случайных процессов являются их автокорреляционные функции. По определению (см., например, Заездный А.М. "Основы расчетов по статистической радиотехнике", М., Связь-издат, 1969 г., стр. 94) автокорреляционная функция B() процесса вычисляется как обратное Фурье-преобразование от его энергетического спектра S (F): где - полярный радиус пространственных координат изображения Значение автокорреляционной функции в нуле B( = 0) равно средней мощности процесса, т.е. сумме мощностей постоянной и переменной составляющих 2. За меру флюктуации процесса по пространственным координатам принято отношение средней мощности процесса к мощности переменной составляющей: Чем больше влажность ЛГМ, тем меньше средняя радиояркостная температура и больше флюктуации амплитуды ИК-сигнала, следовательно, отношение мало. И наоборот, чем меньше влажность, тем больше значение и меньше (2), а упомянутое отношение велико. Чувствительность таким образом введенного относительного показателя к измеряемому параметру (W) оказывается более высокого порядка, поскольку учитывает оба селективных признака ИК-сигнала, а также позволяет ограничиться расчетом одного значения автокорреляционной функции B( = 0). Чем выше КПО, тем меньше скорость изменения мощности восходящего ИК-излучения. На фиг. 2 иллюстрируется зависимость изменения параметров ИК-сигнала от возрастания КПО. В дифференциальном виде эта зависимость представляется в виде: , 2 - коэффициент пропорциональности. Общее решение дифференциального уравнения: L = C0+2ln KПO; Co - постоянная интегрирования, определяемая из начальных условий. Совместное решение уравнений W (КПО); L (КПО): позволяет исключить промежуточную переменную (КПО) и получить закон изменения функции влажности (W) ЛГМ непосредственно от параметров регистрируемого сигнала, т.е. от его относительного показателя: Верификация изложенных аналитических соотношений осуществлялась на основе обработки реальных ИК-сигналов изображений лесных массивов. Процедура расчетов и регрессионная зависимость влажности ЛГМ от относительного показателя изложена в примере реализации способа. Конечной целью заявляемого способа является зонирование территории лесных массивов по градациям влажности ЛГМ. Последнее достигается путем выделения контуров (контурного рисунка) на изображениях, осуществляемое методами пространственного дифференцирования (см. , например, Дуда Р.О., Харт П.Е. "Распознавание образов и анализ сцен", перевод с англ., из-во Мир, М., 1976 г. "Пространственное дифференцирование", стр. 287-288). Существует несколько известных стандартных операторов (Робертса, Лапласа, Собела), позволяющих вычислять контура на двумерных изображениях. В частности, например, перекрестный оператор Робертса рассчитывают в окне (маске) 2х2 элемента, для каждой дискретной точки изображения A(x,y) из соотношения: Вычисленное значение оператора сравнивают с установленным порогом R (i, j) > П, выводят на экран точки, для которых R (i, j) превышает порог. Выделение контуров на изображении с использованием масок различных операторов представляется стандартной математической операцией, входящей в комплект специализированного программного обеспечения (см., например, МАТН САД. 6.0. PLVS, издание второе, стереотипное, информационно-издательский дом "Филинъ", 1997 г., стр. 50-68). На фиг. 3 представлены результаты программной обработки ИК-изображения лесного региона, полученного со спутника NOAA (США) (при рекомендуемых значениях пороговой величины градиента 6,3). В оконтуренных таким образом участках изображения вычисляются моменты электрического сигнала и на основе регрессионной зависимости рассчитывают уже реальную влажность ЛГМ участков на момент съемки. Пример реализации способа. Заявляемый способ может быть реализован на базе устройства по схеме фиг. 4. Функциональная схема устройства фиг. 4 содержит систему орбитальных спутников наблюдения 1 типа NOAA (США) с установленными на них сканирующими радиометрами 2 типа AVHRR. Радиометры осуществляют прием выходящего ИК-излучения подстилающей поверхности в полосе сканирования 3. Зарегистрированный радиометром сигнал с сопутствующей информацией (время приема, эфемириды) в режиме открытого доступа принимается наземными пунктами приема 4, расположенными в г. Красноярске, г. Москве (ИКИ), где записывается на видеомагнитофон 5 типа "Арктур". Из принятых изображений подстилающей поверхности в пункте приема 4 на основе сопутствующей информации формируется база данных изображений лесных массивов, которые в виде файлов помещаются в Интернет 6 (см. , например, сервер, www http:, smis. IkI. RSSI. RV.). По запросу потребителей информации в виде кадров изображений лесных массивов перекачивается из Интернет в индивидуальную систему накопления, хранения и ввода информации 7, организуемую на базе стриммеров типа FT-120. Обработка массивов данных ИК-изображений осуществляется на ПЭВМ 8 типа IBM PC 486/487 в комплекте: процессора-вычислителя 9, винчестера 10, оперативного ЗУ 11, клавиатуры 12, дисплея 13, принтера 14, графопостроителя 15. Гидрометеоцентр РФ ежедневно предоставляет службе охраны лесов информационную карту значений метеорологического показателя КПО на всю территорию лесов России. Федеральной службой лесов России ведется карта лесопожарной опасности, которая позволяет выборочно отслеживать наиболее опасные районы, осуществлять маневр силами и средствами, разрабатывать регламент авиапатрулирования (см., например, "База данных лесопожарной опасности" по отчетам Центральной авиабазы", г. Пушкино, 1997-1998 г., сервер "Интернет", http: smis. IkI. RSSI. RV.). Изображения лесных массивов, полученные из "Интернет" с известным временем и координатами съемки привязывают к местным шкалам КПО контрольных площадок. В табл. 3 представлены данные контрольных замеров и результаты программной обработки изображений со спутников NOAA. Для этого предварительно специализированный комплекс программы МАТН САД записывают на винчестер 10 и осуществляют программное выделение контуров на изображениях и расчет числовых характеристик. Поскольку мощность регистрируемого радиометром AVHRR излучения зависит от характеристик тракта зондирования и дальности, то при расчетах используют не абсолютную амплитуду сигнала, а значения квантованных уровней. Отечественные пункты приема регистрируют сигналы со спутника NOAA в шкале 0...256 уровней. Мощности постоянной и переменной составляющих вычисляют как квадраты первого и второго моментов сигнала изображения. Совместное решение уравнений W (КПО) и L (КПО) при данных контрольных замеров (табл. 3) позволило установить начальные условия и вычислить постоянные интегрирования. Получена функция регрессии влажности ЛГМ от параметров ИК-сигнала изображения следующего вида: Коэффициент a для радиолинии AVHRR и лесов Красноярского края России составил величину a1,2. График функции регрессии иллюстрируется фиг. 5. Осуществляя ежедневное обновление информации со спутников и ее автоматизированную обработку с распечаткой на графопостроителе 15, представляется возможным оперативно, достоверно, с точной координатной привязкой отслеживать состояние ЛГМ на обширных площадях.ФОРМУЛА ИЗОБРЕТЕНИЯСпособ контроля лесопожарной опасности, включающий регистрацию собственного излучения подстилающей поверхности в диапазоне, соответствующем максимуму теплового ИК-излучения подстилающей поверхности, калибровку тракта зондирования по измерениям эталонных участков, вычисление числовых характеристик и расчет результирующего показателя, отличающийся тем, что в качестве результирующего показателя используют непосредственно величину влажности W лесных горючих материалов, преобразуют зарегистрированную функцию электрического сигнала в цифровые матрицы отсчетов зависимости амплитуды от координат, выделяют методами пространственно дифференцирования контуры на двумерных изображениях лесных массивов контролируемого региона при установленных пороговых значениях градиента, а величину влажности лесных горючих материалов внутри выделенных контуров рассчитывают по регрессионной зависимости где a - поправочный коэффициент, учитывающий параметры тракта зондирования, географическую зону, тип лесов; 2 - мощность переменной составляющей сигнала участка изображения внутри контура, Вт; мощность постоянной составляющей сигнала участка изображения внутри контура, Вт.Популярные патенты: 2175189 Способ регенерации растений сорго в культуре in vitro ... активность культур при пассировании (Эльконин Л.А., Тырнов В.С., Папазян Н.Д., Ишин А.Г. Морфогенез и стабильная регенерация растений в каллусных культурах, полученных от зрелых зародышей Sorghum sp. (Poaceae) // Ботан. журн. - 1989. - Т. 74, N 12. - С. 1740-1746). Наиболее близким к предлагаемому является способ регенерации растений сорго из эмбриогенных каллусных культур (Elkonin L.A., Lopushanskaya R.F., Pakhomova N. V. Initiation and maintenance of friable embryogenic callus of sorghum (Sorghum bicolor L. Moench) by amino acids // Maydica, 1995. - V. 40. - P. 153-157), заключающийся в культивировании первичных эксплантов (незрелых зародышей, фрагментов молодых метелок) на ... 2127256 Замещенные простые оксимовые эфиры и фунгицидное, инсектицидное, арахноицидное средство ... I.017, I.058, I.086, I.096 и I.029 показала пороги действия от 200 до 1000 частей на миллион. B.2. Nephotettix cincticeps (зеленые рисовые цикады), контактное действие Круглые фильтры обрабатывают водным составом активного вещества и затем накрывают ими 5 взрослых цикад. Через 24 часа определяют процент смертности. В этом испытании соединения I.007, I.011, I.014, I.015, I.003, I.002, I. 004, I.017, I.117, I.307, I.192, I.193, I.195 и I.201 показали порог действия от 0,4 до 0,1 мг. B. 3. Prodenia litura (Египетский хлопковый червь), опыт по культивированию Пять личинок в стадии развития L3 (10-12 мм) помещают в стандартную питательную среду (3,1 л воды, 80 г агара, 137 г пивных ... 2438304 Улей ... ульи-лежаки [5, 6]. В них в зависимости от конструкции передняя и задняя стенки выполнены двойными с утепляющей засыпкой между внутренней и наружной обшивками, а также предусмотрены утеплительные подушки, холстики, разделительные диафрагмы, сетки вентиляционные и т.д. Наиболее близким аналогом является улей, содержащий дно, по меньшей мере, один корпус с рамками, подкрышник, крышу и отъемные наружные передние, задние, боковые и донный утеплительные блоки [7]. Технология производства приведенных примеров ульев 3, 4, 5 достаточно сложная. В практике за счет утеплителей они имеют большие размеры и большой вес (порядка 30-35 кг). Эти ульи громоздки, тяжелы и неудобны. Все это ... 2166252 Способ удаления костного мозга из губчатых костных трансплантатов ... круговыми движениями волновода диаметром 6 мм, равномерно озвучивая весь объем спонгиозы, погруженный в жидкость. Частота ультразвуковых колебаний составляет 24,5 - 28,5 кГц, продолжительность ультразвукового воздействия составляет 1 минуту. Далее проводят в течение 0,5 часа промывку губчатых костных трансплантатов проточной водой, затем помещают последние в спирто-эфирный раствор, выдерживая их в нем в течение 0,5-1,5 часов. Ультразвуковые колебания низкой частоты, воздействуя на биосистемы, вызывают кроме кавитации возникновение акустических течений, микротоков и капиллярных волн. Акустические течения и микротоки способны порождать вихревые движения, приводящие к разрыву ... 2253227 Устройство для регулирования температуры в улье ... внешней температуры, задатчик внешней температуры, датчик внутриульевой температуры, задатчик внутриульевой температуры, силовой коммутирующий элемент, светоизлучатель, фотоприемник, нагревательный элемент, клеммы для подключения источника питания, первый вход нагревательного элемента соединен с первым входом силового коммутирующего элемента, а второй вход - с первой клеммой для подключения источника питания, второй вход силового коммутирующего элемента соединен со второй клеммой для подключения источника питания, первые входы датчиков внутриульевой температуры и внешней температуры объединены и соединены с первым входом светоизлучателя, второй вход датчика внешней температуры ... |
Еще из этого раздела: 2059362 Установка для выращивания мидий 2196418 Устройство для укладки, сушки и хранения прессованного сена и соломы в рулонах 2472336 Соломорезка и оснащенная такой соломорезкой уборочная машина 2264065 Способ возделывания сельскохозяйственных культур на корм 2281637 Способ производства зеленого корма при возделывании в орошаемом земледелии и устройство для его осуществления 2019090 Самонапорная оросительная система 2413409 Способ и устройство для уплотнения убранной массы для получения силоса 2154940 Способ получения, содержания и хранения живого корма для биологических объектов птиц и рыб 2389173 Способ выращивания земляники садовой 2057432 Биологический состав кузнецова для подсочки деревьев, в том числе каучуконосов (варианты), и способ его приготовления |