|
(21), (22) Заявка: 2002122145/09, 14.08.2002
(24) Дата начала отсчета срока действия патента:
14.08.2002
(43) Дата публикации заявки: 27.02.2004
(45) Опубликовано: 10.08.2004
(56) Список документов, цитированных в отчете о поиске:
ТОНКОНОГОВ М.П., ВЕКСЛЕР В.А., БИРЖАНОВ К.Ж. Диэлектрическая релаксация в водных растворах и суспензиях. Изв. вузов. Физика, 1975, №2, с.81-84. RU 2134425 C1, 10.08.1999. RU 2089889 C1, 10.09.1997. RU 3829 U1, 16.03.1997. WO 9221983, 29.06.1993. ЕР 0478815 В1, 08.04.1992. US 3927369 А, 16.12.1975.
Адрес для переписки:
625003, г.Тюмень, ул. Семакова, 10, Тюменский государственный университет
|
(72) Автор(ы):
Семихина Л.П. (RU)
(73) Патентообладатель(и):
Тюменский государственный университет (RU)
|
(54) СПОСОБ ОПРЕДЕЛЕНИЯ ДИЭЛЕКТРИЧЕСКИХ ПАРАМЕТРОВ ВОДЫ И ЕЕ РАСТВОРОВ В НИЗКОЧАСТОТНОЙ ОБЛАСТИ С ПОМОЩЬЮ L-ЯЧЕЙКИ
(57) Реферат:
Изобретение относится к физическим методам исследования состояния воды и ее растворов в различных объектах и может использоваться при решении фундаментальных и прикладных проблем водных систем. Предлагается способ определения диэлектрических параметров воды и ее растворов с удельной электропроводностью < пред.=20 мСм/см в диапазоне частот 10 кГц-100 МГц с помощью комплекта соленоидальных катушек индуктивности идентичного размера (L-ячеек), подключаемых к колебательному контуру куметра. Благодаря тому, что исследуемый объект вводится в L-ячейки в диэлектрическом сосуде так, что зазор между объектом и обмоткой измерительной ячейки составляет 4-6 мм, реализуется индуктивный принцип работы L-ячеек за счет возникновения внутри исследуемого объекта токов смещения, зависящих от диэлектрических параметров объекта. Диэлектрические параметры объекта находятся в слабом вихревом электрическом поле L-ячейки, причем значение тангенса угла диэлектрических потерь tg водного объекта на данной частоте рассчитывается по определенному соотношению для tg , а диэлектрическая проницаемость характеризуется величиной отн и рассчитывается по определенному соотношению для отн. Найденные значения tg и отн не зависят от добротности и индуктивности измерительной L-ячейки, если длина 1 и радиус r ее обмотки удовлетворяют условию 1/r 7, r 20 мм, а максимально возможная величина напряженности вихревого электрического поля внутри ее равна E*max=5-300 мкВ/см. Способ измерения диэлектрических параметров водных объектов является методом неразрушающего контроля их состояния, поскольку напряженность электрического поля E*max внутри L-ячейки на несколько порядков меньше напряженности поля внутри измерительных конденсаторов в общеизвестном емкостном способе измерения этих же параметров. В результате предложенным L-способом выявляются принципиально новые сведения о свойствах водных систем, например, в диапазоне частот 10 кГц-40 МГц фиксируется низкочастотная дисперсия диэлектрической проницаемости и связанный с ней низкочастотный максимум tg . Технический результат заявленного изобретения заключается в возможности кроме tg находить и диэлектрическую проницаемость жидкости , причем значение tg и не зависят от параметров измерительной ячейки. 6 ил., 2 табл.
Изобретение относится к физическим методам исследования состояния воды и ее растворов в различных объектах и может использоваться при решении фундаментальных и прикладных проблем водных систем. Учитывая, что воду содержат все биообъекты, предлагаемый способ представляет также интерес для биологии и биофизики.
Известен емкостной или С-способ определения диэлектрических параметров веществ, в котором исследуемый объект помещается в измерительный конденсатор (С-ячейку) [1]. Для исследования водных растворов на частотах менее 1 МГц этот способ неприменим из-за увеличения на таких частотах потерь обусловленных проводимостью.
Известен также способ бесконтактной кондуктометрии водных растворов, в котором измерительной ячейкой является соленоидальная катушка индуктивности (L-ячейка) [2-6]. Несмотря на то, что данный способ используется уже более 50 лет, даже принцип работы L-ячеек остался спорным. Отмечается плохое согласие эксперимента с теоретическими моделями, неоднозначность получаемых таким способом значений проводимостей жидкостей. Большинство авторов предполагают, что изменения индуктивного и активного сопротивления L-ячейки после введения жидкости обусловлены возникновением в ней вихревых токов.
Наиболее близким аналогом предлагаемого способа является способ определения тангенса угла диэлектрических потерь tg м водных растворов солей в диапазоне частот 100 кГц – 10 МГц, предложенный в [7], в котором исследуемый водный раствор помещается в соленоидальные катушки индуктивности (L-ячейки), подключенные к колебательному контуру куметра, а величина tg м рассчитывается из соотношения
где Q1; C1; Q2; C2– добротность и емкость колебательного контура куметра при резонансе без жидкости и с жидкостью внутри измерительной катушки индуктивности. Найденные данным способом потери в [7] были условно названы “магнитными” в отличие от обычных диэлектрических потерь tg , определяемых С-способом. Исследуемая жидкость в [7] с помощью насоса последовательно прокачивалась через измерительные ячейки, каркас которых сечением S=1,74 см2 изготавливался из тонкого стекла пирекс, поэтому зазор между обмоткой катушки и жидкостью равнялся толщине тонкой стеклянной трубки. Как отмечается в [7, стр.82], значения tg м характеризуют некие структурные релаксационные процессы в водных растворах и не зависят от добротности и индуктивности измерительной катушки, а также от сечения и формы пробы исследуемой жидкости лишь начиная с некоторой частоты, зависящей от проводимости и объема сердечника. Влияние сечения и проводимости пробы на величину tg м объясняется возникновением внутри жидкости вихревых токов, для снижения которых рекомендовано использовать L-ячейки с малым сечением (менее 7 мм). Однако даже в таких ячейках значения tg м водных растворов в [7] существенно отличались от значений tg . Недостаточная теоретическая и экспериментальная разработка предложенного в [7] способа измерения диэлектрических потерь водных растворов, а также зависимость tg м от параметров L-ячейки резко снижают его практическую и научную значимость.
Предлагаемый способ измерения диэлектрических параметров водных объектов, как и в прототипе, основан на измерении изменений параметров L-ячейки после введения в нее водных объектов методом куметра. Однако в отличие от прототипа он позволяет кроме tg находить и диэлектрическую проницаемость жидкости , причем значения tg и не зависят от параметров измерительной ячейки.
Прежде всего, в предлагаемом способе по сравнению с прототипом изменяется способ введения исследуемой жидкости в измерительную ячейку: жидкость не прокачивается через них, а вводится в цилиндрическом диэлектрическом сосуде (например, стеклянной или тефлоновой пробирке). Это не только упрощает методику эксперимента, но и за счет увеличения расстояния между обмоткой катушки и исследуемой жидкостью меняет принцип изменения параметров L-ячейки.
Суть предлагаемого способа иллюстрирует фиг.1, согласно которому после внесения водного раствора в диэлектрическом сосуде внутрь катушки индуктивности, подключенной к колебательному контуру куметра, изменяется величина емкости С переменного калиброванного конденсатора, при которой наблюдается резонанс в колебательном контуре куметра, а также существенно (в несколько раз) уменьшается добротность Q при резонансе. (Эффект от внесения пустого сухого сосуда пренебрежимо мал).
Поскольку при резонансе емкость С и индуктивность L колебательного контура с активным сопротивлением R и добротностью Q= L/R удовлетворяют соотношению
то при добротности колебательного контура Q>20 с точностью не ниже 0,1% влиянием величины активной составляющей контура и ее изменением при внесении жидкости на резонансную емкость С можно пренебречь.
Если учесть наличие у катушки индуктивности собственной емкости С0 и пренебречь вкладом активного сопротивления контура, то условие резонанса в контуре будет иметь вид
Из соотношения (3) следует, что емкость С при резонансе после внесения жидкости в L-ячейку может измениться лишь по двум причинам:
во-первых, из-за изменения собственной емкости ячейки С0 обмотка катушки в этом случае выполняет функцию внешних электродов подобно электродам обычной С-ячейки (емкостной принцип работы ячейки);
во-вторых, из-за изменения индуктивности катушки (индуктивный принцип работы ячейки).
Оценить вклад емкостного и индуктивного принципа изменения параметров L-ячеек можно путем сопоставления их с модельными CL-ячейками с тем же числом витков, нанесенных на идентичный корпус, например цилиндрическую тонкостенную стеклянную мензурку. Электродами емкостной СL-ячейки являются 2 незамкнутые между собой провода, намотанные параллельно друг другу по спирали. Эквивалентом методики ввода жидкости в измерительную ячейку в прототипе является заливка водного раствора непосредственно внутрь ее корпуса. Поскольку при этом параметры CL– и L-ячейки меняются практически идентично, то в прототипе преобладает емкостной принцип работы L-ячеек. Так как рассчитать изменение распределенной собственной емкости катушки с многослойной обмоткой при внесении в нее жидкости практически невозможно, то L-ячейки, работающие по емкостному принципу, не представляют интереса. Между тем, все предыдущие попытки использования L-ячеек для измерения диэлектрических параметров основывались именно на емкостном принципе ее работы [2-6].
При удалении жидкости от обмотки L-ячейки емкостной механизм изменения ее параметров резко снижается. В предлагаемом способе жидкость вводится в L-ячейку так, что между обмоткой катушки и жидкостью находится примерно 2-2,5 мм материала корпуса катушки, затем 1-1,5 мм воздушного зазора и 1,5-2 мм стенок пробирки, итого 4-6мм. При таком зазоре между обмоткой и водным раствором чувствительность СL-ячейки к вводу жидкости оказывается существенно ниже, чем у L-ячейки – фиг.1. При указанном выше расстоянии между обмоткой ячейки и жидкостью емкость СL-ячейки из-за пробирки с водным раствором возрастает всего на 5% от собственной емкости ячейки. Поэтому для L-ячеек с собственной емкостью С0=(2-6)пФ емкостной эффект примерно на порядок меньше наблюдаемого в эксперименте сдвига резонансной емкости и им можно пренебречь.
Итак, в отличие от прототипа в предлагаемом способе реализуется индуктивный принцип работы L-ячеек. Следовательно, после введения жидкости происходит изменение активного R и индуктивного  L сопротивления измерительной L-ячейки, отношение которых определяет величину тангенса угла диэлектрических потерь жидкости tg = R/ L.
Значения R и  L при индуктивном принципе работы L-ячейки с собственной емкостью С0 можно найти из соотношений
где R0 и L0 – активное и индуктивное сопротивление контура до введения жидкости, C1, С2; Q1, Q2 – значения емкости калиброванного конденсатора и добротности контура при резонансе до и после введения жидкости в L-ячейку. Добротности контура Q1, Q2 соотношениями (6, 7) связаны со значениями Q 1, Q 2 так называемыми множителями вольтажа колебательного контура, отсчитываемыми при резонансе по шкале вольтметра, измеряющего напряжение на калиброванном конденсаторе куметра. При С0=0, Q1,2=Q 1,2.
Из соотношений (4-7) получаем
Соотношение (10) оказывается идентичным (1) лишь при С0=0, причем различия в значениях числителя и знаменателя выражений (10) и (1) при С0=(5-10)пФ достигают 20%-30%. Однако различие значений tg , рассчитанных по соотношениям (10) и (1) при С0<10 пФ, не превышают 2 10-3, что в несколько раз ниже погрешности эксперимента – таблица 1.
Поэтому при С0<10пФ для расчета tg нет смысла пользоваться более сложным соотношением (10), вполне применимо и упрощенное соотношение (1). Независимость значений tg от величины С0 является дополнительным подтверждением отсутствия заметного вклада емкостного принципа изменения параметров L-ячейки в предлагаемом способе. Поскольку соотношение (1) справедливо лишь при индуктивном принципе работы L-ячейки, то его использование в прототипе некорректно.
Значения тангенса угла потерь tg жидкостей, найденные предложенным способом, в отличие от прототипа не зависят ни от индуктивности, ни от добротности измерительной ячейки. В качестве примера в таблице 1 приведены экспериментальные данные для бидистиллированной воды, полученные данным способом на одной и той же частоте 100 кГц на трех L-ячейках одинакового размера, но с разной индуктивностью (диаметр пробы 20,5 мм).
Возможность нахождения с помощью L-ячейки не только значений tg , но и диэлектрической проницаемости жидкости, основывается, во-первых, на обнаруженном экспериментальном факте, что для немагнитных жидкостей величина C=C1-С2>0, равная изменению резонансной емкости колебательного контура куметра после ввода жидкости в L-ячейку, определяется лишь исследуемой жидкостью и не зависит от индуктивности измерительной ячейки – таблица 1. И, во-вторых, на выяснении природы изменения индуктивности L-ячейки с немагнитными жидкостями и теоретическом доказательстве, что в случае индуктивного принципа работы L-ячейки величина С при одинаковом размере проб жидкости определяется лишь их диэлектрической проницаемостью.
В прототипе и в других работах по L-ячейке полагается, что причиной изменения индуктивности ячейки с немагнитными жидкостями является возникновение внутри жидкости вихревых токов [2-7]. Между тем, учитывая, что индуктивность катушки Z, индукция магнитного поля В и магнитный поток Ф внутри ее связаны между собой соотношением Ф=BSN=LI, где S и N – площадь сечения и число ее витков, I – ток через катушку, а поле вихревых токов по правилу Ленца может только уменьшать магнитное поле катушки, то возникновение вихревых токов внутри жидкости должно приводить к уменьшению индуктивности L-ячейки.
В то же время из (4-5) имеем
Откуда при C=C1-C2>0 получаем, что L=L-L0>0. Следовательно, введение немагнитной жидкости в L-ячейку приводит к увеличению индуктивности ячейки. Т.о. гипотеза об определяющей роли вихревых токов в принципе работы L-ячейки не согласуется со знаком наблюдаемых в эксперименте изменений ее параметров.
То, что в случае водных систем изменение индуктивности L-ячейки обусловлено не вихревыми, а токами смещения, выявляется после нахождения магнитной индукции В=В0J0(аr)/J0(аr0) и магнитного потока Ф=2 В0r0J1(аr0)/аJ0(аr0) внутри цилиндрической бесконечно длинной ячейки, заполненной жидкостью с магнитной проницаемостью , диэлектрической проницаемостью и удельной электропроводностью на частоте , где В0 – поле внутри ячейки без жидкости, a2= ( -i ), J0(ar), J1(ar)- функции Бесселя, r0 – радиус пробы жидкости. Разложив функции Бесселя в ряд и избавившись от комплексных величин, получаем, что Ф=Ф0+ Фсмещ.– Фвихр., где Ф0, Фсмещ. и Фвихр.– магнитный поток в отсутствие жидкости и его изменения за счет токов смещения и вихревых токов в жидкости. Знак изменения L, наблюдаемый в эксперименте, совпадает с эффектом от токов смещения, а от вихревых токов -противоположен. Вкладом от вихревых токов можно пренебречь, если
Для проб водных растворов с r0=1,5 см условие (12) будет выполняться при
Соотношения (12-13) определяют верхний предел удельных электропроводностей водных растворов, для которых применим предлагаемый L-способ измерения их диэлектрических свойств. Пренебрегая вкладом от вихревых токов для таких растворов, получаем

где = lE2max/2(UQ1)2– есть постоянная измерительной ячейки; Еmax= В0r/2 – максимальная напряженность электрического поля внутри L-ячейки радиусом r, U – напряжение, подаваемое на вход колебательного контура и не меняющееся во время измерений, l – длина L-ячейки. Из соотношения (15) следует, что для данной ячейки величина С зависит только от квадрата радиуса пробы жидкости r20 и диэлектрической проницаемости жидкости . Линейную зависимость С от при r0=const подтверждает таблица 2, а от r2 0– фиг.2.
Измерения диэлектрических параметров с помощью L-ячейки упрощаются из-за наличия “эффекта насыщения”, благодаря которому значения С и tg водных растворов не меняются при изменении поля внутри L-ячейки в диапазоне Е*max=5-300 мкВ/см. Независимость величин С и tg от индуктивности и добротности измерительной ячейки есть следствие данного эффекта и происходит лишь если внутри ее *mах=E*mах.
При Е*max=Е*max и r0=const величина С определяется лишь величиной диэлектрической проницаемости исследуемой жидкости, а рассчитанные по (1) значения tg – есть значения тангенса угла диэлектрических потерь этой жидкости в очень слабом электрическом поле E*max измерительной L-ячейки.
Оптимальные размеры измерительных L-ячеек выявляются из данных, представленных на фиг.2-3. Из фиг.2 следует, что рекомендованное в [7] уменьшение радиуса пробы жидкости приводит к снижению чувствительности способа. При r0<10 мм искажается и зависимость между r2 0и С. Оптимальный размер радиуса пробы в предлагаемом способе составляет примерно 15 мм, т.о. оптимальный радиус обмотки измерительной L-ячейки r равен приблизительно 20 мм. Дальнейшее увеличение r приводит к увеличению сопротивления L-ячейки и снижению поля внутри ее ниже E*max на низких частотах.
Влияние длины обмотки l соленоидальной измерительной катушки индуктивности на полученные с ее помощью значения tg и С исследуемого объекта показано на фиг.3 для L-ячейки с r=19 мм. Как видим, при длине катушки l>l0, где l0 – некоторое критическое значение длины катушки, отмеченное на фиг. 3 пунктиром, величина tg исследуемого объекта практически не зависит от l, а между С и l наблюдается четкая линейная зависимость, подтверждающая соотношение (15). Очевидно, что для практических измерений пригодны лишь ячейки с l>l0. При l0 соотношение (15) уже не выполняется, а значение tg существенно уменьшается. Критической длине катушки l0 соответствует условие
Оно совпадает с условием, при котором магнитное поле внутри соленоида можно считать равным полю бесконечно длинного соленоида. В измерительных катушках меньшей длины происходит занижение значений tg . B катушках, удовлетворяющих условию (16), возрастает область однородного магнитного поля. На необходимость высокой однородности магнитного поля внутри измерительной ячейки указывает также то, что при неравномерной обмотке L-ячейки значения tg , найденные с ее помощью, заметно снижаются, даже если ее длина удовлетворяет условию (16).
Для получения частотной зависимости tg и объектов в широком диапазоне частот необходим комплект L-ячеек одинакового размера с перекрывающимся частотным диапазоном и комплект куметров (например, куметры TESLA BM-311, TESLA ВМ-560, ВМ 409G). При выполнении указанных выше условий различия между значениями tg и С объектов с < npeд., найденные на разных ячейках, не превышают погрешности эксперимента – фиг.4.
Индикатором того, что объект имеет удельную электропроводность > npeд. и предложенный способ для его исследования непригоден из-за увеличения вклада вихревых токов, является возникновение существенных различий между значениями tg , найденными с помощью L-ячеек с разной индуктивностью – фиг.5. Эксперимент подтверждает полученную по (12-13) оценку величины пред., ограничивающей применимость L-способа измерения tg – фиг.5.
Особенностью предложенного способа является то, что с его помощью выявляется низкочастотный максимум tg и низкочастотная дисперсия диэлектрической проницаемости водных объектов, которая фиксируется по резкой частотной зависимости величины С – фиг.4. Для растворов с < пред. изменение С завершается на частоте менее 10 МГц. Поэтому частотные изменения водных объектов наиболее просто характеризовать величиной отн= C/ C , где С , – сдвиг резонансной емкости колебательного контура куметра после введения исследуемого объекта на частоте (1-10) МГц. Одновременное фиксирование значений tg и отн позволяет получить более полную информацию о свойствах водных объектов.
То, что определяемые предложенным способом потери – есть обычные диэлектрические потери, подтверждает совпадение во всем исследованном диапазоне частот частотных зависимостей tg ( ), полученных общеизвестным С- и предложенным L-способом для спиртов – фиг.6. Однако для водных объектов совпадение значений tg , полученных L- и С-способом наблюдается лишь на высоких частотах – фиг.6.
Совпадение значений tg , полученных L- и С-способом для жидкостей без сплошной сетки водородных связей, например спиртов, а также подобие частотных зависимостей tg ( ), полученных С-способом для воды и спиртов и их принципиальное различие в случае L-способа (фиг.6), указывает на искажение состояния воды внутри С-ячейки. В L-ячейке жидкость подвергается воздействию вихревого электрического поля E*max, напряженность которого на 3-5 порядков меньше, чем в С-ячейке, поэтому его влияние на состояние воды гораздо меньше. Именно этот факт является причиной различия значений tg воды и ее растворов, полученных L- и С- способом. Т.о. предлагаемый L-способ измерения диэлектрических параметров в отличие от С-способа является методом неразрушающего контроля водных систем. Поэтому он представляет особый интерес для исследования именно водных объектов.
Введение в воду примесей, увеличивающих ее удельную электропроводность, смещает область дисперсии и частоту, на которой наблюдается максимум tg , в область более высоких частот. Именно по этой причине наблюдается влияние проводимости жидкости на изменение параметров L-ячейки. Однако взаимосвязь между параметрами L-ячейки и проводимостью жидкости не прямая, не линейная и не однозначная и обусловлена влиянием ионов примесей на состояния воды. Поэтому, несмотря на полувековую историю разработки способа применения L-ячейки как способа бесконтактной кондуктометрии, приходим к выводу о его нецелесообразности.
Однако предлагаемый способ применения L-ячейки как L-способ определения диэлектрических параметров водных растворов представляет несомненный интерес, так как позволяет получить принципиально новые сведения о свойствах воды. Например, как отмечено выше, с его помощью в водных объектах выявляется низкочастотная дисперсия диэлектрической проницаемости и низкочастотный максимум tg , не фиксируемый общеизвестным емкостным или С-способом измерения диэлектрических параметров.
Предлагаемый способ обладает всеми признаками, предъявляемыми к изобретению. В доступной литературе нет технического решения с данной совокупностью признаков и вполне может быть выдан патент на него.
Литература
1. Лопатин Б.А. Теоретические основы электрохимических методов анализа. М., Высшая школа, 1975, 295 с.
2. Андреев B.C. Кондуктометрические методы и приборы в биологии и медицине. М., Медицина, 1973, 336 с.
3. Заринский В.А. Высокочастотный химический анализ. М., Наука, 1970, 200 с.
4. Лопатин Б.А. Кондуктометрия. Новосибирск, Изд-во СО АН СССР, 1964, 280 с.
5. Лопатин Б.А. Теоретические основы электрохимических методов анализа. М., Высшая школа, 1975, 295 с.
6. Лопатин Б.А. Высокочастотное титрование с многозвенными ячейками. М., Наука, 1980, 207 с.
8. Ахадов Я.Ю. Диэлектрические свойства чистых жидкостей. М., 1972, 412 с.
9. Хиппель А.Р. Диэлектрики и их применение. М., Госэнергоиздат, 1959, 336 с.

Формула изобретения
Способ определения диэлектрических параметров воды и ее растворов в диапазоне частот 10 кГц-100 МГц с помощью комплекта соленоидальных катушек индуктивности идентичного размера (L-ячеек), подключаемых к колебательному контуру куметра, отличающийся тем, что исследуемый объект с удельной электропроводностью < пред=20 мСм/см, где пред – верхний предел удельных электропроводностей водных растворов, вводится в L-ячейки с собственной емкостью С0<10 пФ в диэлектрическом сосуде так, что зазор между объектом и обмоткой измерительной ячейки составляет 4-6 мм, длина 1 и радиус r обмотки L-ячеек удовлетворяют условию 1/r 7, r 20 мм, а максимально возможная величина напряженности вихревого электрического поля внутри ячеек поддерживается в пределах 5-300 мкВ/см, при этом реализуется индуктивный принцип работы L-ячеек за счет возникновения внутри исследуемого объекта токов смещения, зависящих от диэлектрических параметров объекта, в таком случае значение тангенса угла диэлектрических потерь объекта рассчитывается по соотношению
где R – изменение активного сопротивления L-ячейки;
 L – изменение индуктивного сопротивления L-ячейки;
Q1, C1, Q2, C2 – значения добротности и емкости колебательного контура куметра при резонансе до и после помещения жидкости внутрь L-ячейки,
а диэлектрическая проницаемость характеризуется величиной отн= С/ С , где С и С – сдвиг резонансной емкости колебательного контура куметра после введения исследуемого объекта на частоте измерения и на частоте порядка 10 МГц.
РИСУНКИ
MM4A – Досрочное прекращение действия патента СССР или патента Российской Федерации на изобретение из-за неуплаты в установленный срок пошлины за поддержание патента в силе
Дата прекращения действия патента: 15.08.2007
Извещение опубликовано: 10.03.2009 БИ: 07/2009
|
|