Индуктивный способ электромагнитного мониторинга процесса оттайки грунта
Реферат
Изобретение относится к электромагнитным методам исследования вещества земли. Оно может быть использовано при геокриологических и других исследованиях. Область преимущественного применения - определение электропроводности и поляризационной характеристики мерзлых грунтов в естественных условиях. Технический результат - повышение точности определения сопротивления геологической среды и дополнительного определения поляризационных характеристик исследуемой среды. Сущность изобретения: возбуждают горизонтальной рамкой переменное электромагнитное поле, измеряют горизонтальную (Нr) и вертикальную (Нz) составляющую поля на нескольких частотах в диапазоне десятки - первые сотни килогерц, определяют по отношению Нz/Нr электросопротивление земли. Измерение Нr на всех рабочих частотах f1, f2, ... , fi, ... , fN выполняют после компенсации составляющей Нr,o первичного магнитного поля на нижней рабочей частоте f1, находят приращение (lg) = lgi+1-lgi и по изменению во времени i и коэффициента Kf= (lg)/(lgf) судят о состоянии процесса оттайки грунта. 4 ил.
Изобретение относится к электромагнитным методам исследования вещества Земли. Оно может быть использовано при геокриологических, инженерно-геологических, геоэкологических и других исследованиях, когда требуется знать электросопротивление земли в интервале глубин от земной поверхности до нескольких метров в диапазоне частот от нескольких десятков килогерц до первых сотен килогерц. Область преимущественного применения предлагаемого изобретения - определение электропроводности и поляризационной характеристики мерзлых грунтов в естественных условиях, контроль за изменением указанных характеристик грунтов во времени.
Известен индуктивный способ определения электропроводности массивов с использованием разнесенных по земной поверхности излучателя поля и магнитоприемника. В наиболее близком нам способе-аналоге [1] источником поля является горизонтальная рамка или стержневая самоотвешивающая антенна (аппаратура типа ДЭМП), а измеряются взаимно перпендикулярные компоненты магнитного поля - радиальная Hr и вертикальная Hz; электросопротивление среды определяют через отношение Hz/Hr при данной рабочей частоте f и разное r установки. Недостаток этого способа - низкая точность определения поляризационной характеристики мерзлых грунтов из-за сильного мешающего влияния скин-эффекта в более электропроводной, смежной (например, подстилающей, не мерзлой) геологической среде, а также влияния других электромагнитных эффектов (интерференция, дифракция и т.п.) на посторонних геоэлектрических неоднородностях. Для уменьшения влияния указанных нежелательных эффектов в этом способе уменьшают разнос установки для нескольких метров или снижают рабочую частоту. При этом измерения переходят в область малых параметров, где информативно только поле Hr. С применением же известного способа в области малых параметров он имеет недостаток, т.к. регистрируемое поле Hr кроме полезной информации (отклик земли на электромагнитное возбуждение) содержит и составляющую Hr,o первичного поля (помеху). Это происходит вследствие неизбежного нарушения требуемой взаимной ортогональности в системе излучатель (вертикальный магнитный диполь) - датчик поля Hn (например, из-за земной неровной поверхности, неодинаковой высоты излучателя и приемного датчика и т. п. ). По указанной причине при применении способа-аналога с малыми разносами или низкими частотами резко возрастает погрешность в определении геологической среды, и он становится совершенно не применимым для изучения частотной дисперсии электропроводности грунтов. Известен также способ дистанционных зондирований с возбуждением вертикальным магнитным диполем, в котором датчик поля Hr устанавливают по малой оси эллипса поляризации магнитного поля на частоте, соответствующей малым параметрам поля [2]. Но он применим лишь при отношении Hz/Hr, не превосходящем 10 - 30, если же среда имеет высокое сопротивление (мерзлота), а разносы составляют несколько метров, то отношение Hz/Hr может достигать 100 и более; в этом случае для отыскания малой оси эллипса поляризации измерительный датчик Hr необходимо поворачивать в пределах около 1o с точностью в несколько угловых минут, что превосходит возможности применяемых в способе [2] угломерно-поворотных устройств. Из-за этого резко снижается точность измерения поля Hr и сильно возрастает погрешность определения в нужном для изучения поляризуемости мерзлых грунтов диапазоне частот и разносов установки. Для достижения высокой точности измерения отклика земли на электромагнитное возбуждение при малых параметрах поля генераторную рамку и приемные датчики объединяют в одну жесткую конструкцию, переносимую по рабочему профилю как одно целое устройство, а для выделения сигнала, пропорционального только вторичному полю, используют различные вспомогательные схемно-технические приспособления, включаемые в это устройство, как его неотъемлемая часть. Цель предлагаемого изобретения - повышение точности определения электросопротивления и поляризационных характеристик геологической среды при использовании индуктивного метода возбуждения и регистрации поля. Поставленная цель достигается тем, что в индуктивном способе электромагнитного мониторинга процесса оттайки грунта возбуждают горизонтальной рамкой переменное электромагнитное поле в диапазоне частот десятки - первые сотни килогерц, измеряют вертикальную Hz и горизонтальную Hr составляющие на всех рабочих частотах f1, 2,......, fi,...., fN после компенсации составляющей Hr,o первичного магнитного поля на каждой рабочей частоте f, определяют по отношениям Hz/Hr электросопротивление земли , находят приращение (lg) = lgi+1-lgi и по изменению во времени i и коэффициента Kf = (lg)/(lgf) судят о состоянии процесса оттайки грунта, где (lg) и (lgf) соответственно приращения логарифмов от электросопротивления и частоты, Kf - коэффициент частотной дисперсии. На фиг. 1 приведена структурная схема устройства для осуществления предлагаемого способа. В устройство входят: 1 - приемный датчик магнитного поля, способный поворачиваться в вертикальной плоскости и занимать два фиксированных положения - по оси установки и в ортогональном положении , датчик подключается к селективному вольтметру 6; 2 - генераторная рамка с магнитным моментом перпендикулярным оси установки и параллельным плоскости вращения датчика 1; к рамке 2 подключается генератор переменного тока 7; 3 - комплексационная рамка с магнитным моментом направленным по оси установки 4 - коммутатор компенсирующего тока, он изменяет на противоположное направление магнитного поля момента 5 - размыкатель тока компенсационной рамки; в положении А ток в компенсационной рамке отсутствует, а в генераторной рамке сохраняется, в положении Б обе рамки включены последовательно; 6 - жесткая несущая рама. На фиг. 2 показаны результаты применения предлагаемого способа при мониторинге процесса оттайки сезонной мерзлоты в г.Екатеринбурге с 4 марта по 26 июня 1997 г.; на фиг. 2 приведены графики изменения максимального и минимального электросопротивления грунта в диапазоне частот 2,5 - 160 кГц, на фиг. 3 - графики поинтервального коэффициента частотной дисперсии Kf, цифры у графиков соответствуют интервалу частот (в килогерцах), на котором определяли Kf, на фиг. 3 приведен также усредненный график изменения температуры воздуха за время наблюдений - в градусах Цельсия (ToC), на фиг. 4 изображены графики кажущегося сопротивления при изопараметрическом зондировании, выполненном в августе 1997 г. на пункте, где проводилось опробование предлагаемого способа в марте-июне, цифры у точек графика - значение в килогерцах. Предлагаемый способ заключается в следующем (фиг. 1). При установке приемного датчика 1 по оси регистрируемое им поле Hr состоит из трех слагаемых Hr = ImHrHr,0Hk, (1) где ImHr - вторичное поле индуцированных в земле токов; при малых параметрах поля он находится в квадратуре с первичным полем H0, создаваемым генераторной рамкой; Hr,0 - доля первичного поля H0 по направлению возникающая за счет ошибки в перпендикулярности генераторного магнитного момента продольной оси установки и неточного совпадения магнитной оси датчика Hr с осью Hk - поле компенсационной рамки 3. Поскольку в схеме фиг. 1 через генераторную и компенсационную рамки протекает один и тот же ток, то поля Hr,0 и Hk синфазны друг другу и оба одинаково сдвинуты по фазе на /2 по сравнению с ImHr. Следовательно, чтобы регистрируемый сигнал содержал только информативную часть, зависящую от электропроводности земли, должно выполняться векторное равенство. В схеме на фиг. 1 противоположность направления полей и обеспечивается коммутатором 5 направления тока в компенсационной рамке, а равенства их амплитуд добиваются изменением напряженности поля путем перемещения компенсационной рамки по продольной оси установки. В предлагаемом способе продольной осью установки и одновременно жесткой основой, на которой размещены генераторная, компенсационная рамки и приемный датчик поля с поворотным антеннодержателем, является неметаллический брус (рама) длиной 3,4 м. Генераторная рамка и антеннодержатель для удобства при транспортировке выполнены съемными, но имеют на концах рамы посадочные гнезда, чтобы обеспечивать строгое постоянство расстояния между центрами приемного датчика и генераторной рамки, равного 2,8 м. Средняя часть опорного бруса, где размещена компенсационная рамка, имеет размеченную в сантиметрах продольную шкалу. Сама компенсационная рамка для плавности перемещения по брусу выполнена в виде отрезка трубы из фольгированного медью гетинакса длиной 15 см с сечением немного большим, чем у бруска. Размеры компенсационной рамки подобраны так, чтобы при ее перемещении на 1 см шкалы создаваемое ей компенсационное поле в центре приемного датчика было в пределах 10-4 - 10-5 За счет этого в предлагаемом способе достигается высокая точность исключения влияния на измерения нарушений строгой геометрии системы излучатель поля - магнитоприемник. Измерения по предлагаемому способу выполняются следующим образом. 1. Измерения начинают с нижней из выбранных рабочих частот. На этой частоте приемный датчик устанавливают в положение для регистрации поля Hr и при разомкнутом размыкателе 4 измеряют поле Hr. Это соответствует по (1) сумме вторичного поля ImHr и добавки Hr,0. 2. Включают переключателем 4 ток в компенсационной рамке 3. 3. Перемещая компенсационную рамку по несущему брусу, находят ее положение, при котором регистрируемый сигнал минимален. 4. Если отсчет по измерительному прибору не уменьшается, то коммутатором 5 меняют полярность компенсирующего тока и повторяют операцию по п.4. 5. В положении, соответствующем минимуму регистрируемого сигнала, компенсационную рамку закрепляют. В этом положении выполняется равенство (2) и регистрируемый сигнал Urmin пропорционален только вторичному полю ImHr. 6. Измеряют сигнал Urmin. 7. Поворачивают приемный датчик на фиксированный угол 90o и измеряют сигнал Uz. Этот сигнал пропорционален полю Hz, а отношение Uz/Urmin равно искомому отношению полей Hz/Hr. 8. При измерениях на все более высоких частотах операции по пп. 4.5.6 опускают, а измеряют только сигналы Uz и Ur и находят их отношение. Электросопротивление исследуемой среды определяют по известной формуле [1] . Коэффициент частотной дисперсии вычисляют для каждой пары смежных рабочих частот по формуле где индексом i и i+1 обозначены величины, относящиеся к более низкой и более высокой частоте соответственно. Опробование предлагаемого способа проведено при мониторинге процесса оттайки сезонной мерзлоты. Мониторинг проводили в южном лесопарке г. Екатеринбурга в период с 4 марта по 26 июня 1997 г. Возбуждение и измерение поля выполняли на 7 фиксированных частях: 160, 80, 40, 20, 10, 5, 2,5 кГц. В качестве генератора тока и частотно-избирательного микровольтметра использовали аппаратуру малоглубинного частотного зондирования типа МЧЗ-10. На фиг. 2 приведены графики изменения во времени максимального и минимального сопротивления в рабочем диапазоне частот. Максимальная величина уменьшилась от 247-230 Омм (4-11 марта) до 65 Омм 26 июня (примерно в 3,7 раза), минимальное уменьшилось на это же время от 115-125 Омм до 47 Омм (2,6 раза). Очевидно, что такие сильные изменения могут быть вызваны только изменением физического состояния содержащейся в грунте влаги - переходом ее из твердой фазы (лед) в жидкую (вода). На фиг. 3 приведены графики изменения во времени коэффициента частотной дисперсии Kf для шести пар частот (в кГц): 160-80, 80-40, 40-20, 20-10, 10-5, 5-2,5. В нижней части фиг. 3 показан усредненный график изменения температуры атмосферного воздуха. Из сравнения температурного графика и графика Kf видно, что процесс оттайки характеризуется появлением отрицательных значений коэффициента Kf в интервале частот 160-20 кГц. Максимальная величина Kf в отрицательной области, регистрируемая с 3 по 11 апреля на частотах 160-80 кГц, достигает -0,57 - 0,87, что соответствует коэффициенту поляризуемости 50-80%. В период 4 марта по 11 апреля изменения коэффициента Kf отображает неравномерность процесса оттайки мерзлоты, заметную и по изменению погодных условий: начавшееся после 4 марта потепление и быстрое таяние снега сменилось во второй половине марта резким похолоданием с временным установлением снегового покрова; только около 1 апреля вновь установилась теплая погода с устойчиво возрастающей положительной температурой воздуха. Из сравнения графиков Kf, температурного графика (фиг. 3) с графиками на фиг. 2 видим, что изменение коэффициента при частотах 160 - 40 кГц лучше, чем изменение отображает особенности процесса оттайки грунта. По графику Kf на частотах 160-40 кГц можно также заключить, что 4 марта процесс оттайки еще практически не начался (Kf = 0,08), а к 26 июня уже закончился (Kf= 0 - 0,08), в изменении же сохраняется его тенденция к уменьшению и после 26 июня, что, вероятно, связано с обводнением грунта прошедшими дождями. Таким образом, использование при мониторинге дополнительной электрофизической характеристики - коэффициента Kf, получаемого с применением предлагаемого способа, дает важную дополнительную информацию о состоянии процесса оттайки мерзлоты. На фиг. 3 также хорошо видно, что на частотах ниже 20 кГц коэффициент Kf переходит в положительную область, причем положительный знак сохраняется независимо от времени наблюдений. Это показывает, что почастотные измерения при f 20 кГц характеризуют уже не временные изменения состояния вещества грунта, а особенности геоэлектрического разреза в статическом состоянии, а именно, наличие подземного непромерзающего водоносного слоя повышенной электропроводности. В доказательство этого на фиг. 4 приведен график кажущегося сопротивления при изопараметрическом зондировании, выполненном нами по способу [2, стр. 46-49] на том же самом пункте, в августе 1997 г. Цифры у точек графика указывают частоту тока в генераторной рамке в килогерцах. На графике фиг. 4 при разносах 9,8, 15,7 и 20 метров (частота 20,10 и 5 кГц соответственно) наблюдается отчетливый минимум Судя по графику глубина до электропроводного слоя равна 5-7 м, его удельное сопротивление 30 Омм, удельное сопротивление перекрывающего слоя равняется 43 Омм. Таким образом, результаты опробования подтвердили возможность применения предлагаемого способа для определения электросопротивления и коэффициента частотной дисперсии электропроводности мерзлого грунта в процессе его оттаивания. Нам известны другие индуктивные способы, которые позволили бы в естественных условиях определять обе эти характеристики мерзлых грунтов. Обычно сведения о частотной дисперсии электропроводности мерзлых грунтов получают из данных электромагнитных зондирований как изменения не соответствующие теоретической модели горизонтально-слоистой земли. Но, очевидно, что такие изменения могут вызываться и любыми геоэлектрическими неоднородностями, не имеющими к мерзлоте никакого отношения. Соотнесение, связь с мерзлотой не очевидна, и возможна лишь при известной априорной информации о наличии мерзлоты и ее состояния. Достоинство предлагаемого способа состоит в том, что в нем не требуется априорная информация о мерзлоте, он сам дает эту информацию, причем из непосредственных измерений. А это исключает субъективный фактор из получаемого результата. Одним из важных достоинств предлагаемого способа является то, что в нем можно выполнять измерения вторичного поля Hr в области малых параметров при очень малой величине полезного сигнала по сравнению с возбуждающим полем. Так в описанном выше опробовании отношение полезного сигнала (ImHr) к полю на нижней рабочей частоте 2,5 кГц составляло 0,0005 - 0,001. За счет этого достигнута высокая чувствительность предлагаемого способа к почастотным изменениям , вызванным повышенной поляризуемостью мерзлого грунта и повышенной электропроводностью водоносного слоя ниже зоны промерзания. Источники, использованные при составлении заявки: 1. Вешев А.В., Любцева Е.Ф. и др. Временное руководство по методу электромагнитного зондирования с вертикальным магнитным диполем. - МЦМ СССР, М., 1978, с. 19 - 25. 2. Титлинов В.С., Журавлева Р.Б. - Технология дистанционных индуктивных зондирований. - Екатеринбург, УИФ "Наука", 1995, с. 29 - 30.Формула изобретения
Индуктивный способ электромагнитного мониторинга процесса оттайки грунта, отличающийся тем, что в нем возбуждают горизонтальной рамкой переменное электромагнитное поле в диапазоне частот десятки - первые сотни килогерц, измеряют вертикальную Hz и горизонтальную Hr составляющие на всех рабочих частотах f1, f2, ..., fi, ... fn после компенсации составляющей Hr,o первичного магнитного поля на нижней рабочей частоте f1, определяют по отношению Hz/Hr электросопротивление земли, находят приращение электросопротивления (lg) = lgi+1-lgi и частоты (lgf) = lgfi+1-lgfi на смежных частотах и по изменению во времени fi и коэффициента частотной дисперсии kf= (lg)/(lgf) судят о состоянии процесса оттайки грунта.РИСУНКИ
Рисунок 1, Рисунок 2, Рисунок 3, Рисунок 4