Изобретение относится к электромагнитным методам геофизических исследований земной коры и может быть использовано при глубинных зондированиях при поисках и разведке месторождений нефти и газа.
Известен способ геоэлектроразведки, в котором возбуждают импульсное электромагнитное поле в горных породах и измеряют неустановившееся поле на исследуемом участке и в нормальном поле, при этом о строении исследуемого участка судят по разности зависимостей сигналов становления поля, полученных на исследуемом участке и в нормальном поле [1]
Недостатком способа является слабая дифференциация разреза на отдельные слои, т. е. невозможность расчленения всего разреза на отдельные слои при зондировании мощных геологических толщ вследствие того, что в промежуточной зоне сигнал становления формируется полем, обусловленным лишь той частью разреза, которая имеет глубину, соизмеримую с разносом, а также выявление только региональных аномалий при зондировании мощных геологических толщ.
Известен также способ геоэлектроразведки, заключающийся в том, что возбуждают электромагнитное поле в земле с помощью питающей линии, измеряют неустановившееся поле несколькими приемными рамками, расположенными на разных разносах от питающей линии, полученные результаты интегрируют по разносу для каждого фиксированного времени регистрации, затем увеличивают диапазон разносов и снова находят интегральное значение поля для каждого из диапазонов, а геоэлектрические параметры определяют в таком диапазоне разносов, при котором величина интегрального поля выходит на асимптотическое значение [2]
Недостатками способа являются невысокая точность и разрешающая способность из-за усиливающийся осредненности геоэлектрических параметров в пределах выделяемых толщ из-за больших баз интегрирования сигналов в пределах расстановок приемников поля.
Наиболее близким к предлагаемому является способ геоэлектроразведки, основанный на возбуждении и регистрации неустановившегося электромагнитного поля в виде производной вертикальной компоненты вектора магнитной индукции с помощью установки, включающей систему источников поля и группу приемников поля, расположенных в ближней зоне источника на определенной площади, обработке сигналов становления поля и суждении о строении исследуемого участка разреза [3]
Недостатком способа является слабая точность определения параметров неоднородности из-за определения преимущественно интегральных параметров среды, в которых влияние неоднородности проявлялось незначительно. По-существу, способ предназначался только для определения параметров вмещающей среды.
Технический результат изобретения повышение эффективности геоэлектроразведки за счет возможности учета аномальных составляющих электрического и магнитного полей.
Технический результат достигается тем, что в способе геоэлектроразведки, основанном на возбуждении и регистрации неустановившегося электромагнитного поля в виде производной вертикальной компоненты вектора магнитной индукции с помощью установки, включающей систему источников поля и группу приемников поля, расположенных в ближней зоне источника на определенной площади, обработке сигналов становления поля и суждении о строении исследуемого участка разреза, дополнительно регистрируют первые горизонтальные компоненты электромагнитного поля с последующим определением электрических составляющих аномального поля и вторые горизонтальные компоненты электромагнитного поля с последующим определением магнитных составляющих аномального поля, при этом первые горизонтальные компоненты электромагнитного поля регистрируют путем расстановки системы источников поля в количестве не менее двух симметрично исследуемому участку по разные стороны от него и пропускания электрического тока, одинакового по величине и направлению по отношению к исследуемому участку, а вторые горизонтальные компоненты путем пропускания также одинакового по величине электрического тока, но направление которого по одну сторону от исследуемого участка сохраняют первоначальным, а по другую меняют на противоположное, по выявленным параметрам аномальных составляющих электрического и магнитного полей судят о строении исследуемого участка, причем размер площади расстановки приемников, величину шага между приемниками, максимальное и минимальное расстояния между источниками и количество источников системы выбирают из условий
Δy = 0,05h,
N=2n, где n=1, 2, 4, 6, 8, 10, 12
где Q размер площади расстановки приемников,
Δy величина шага между приемниками,
ΔLmax максимальное расстояние между источниками,
ΔLmin минимальное расстояние между источниками,
N количество источников,
Qн размер площади исследуемого участка,
h глубина залегания исследуемого участка,
α погрешность измерений.
Первые и вторые горизонтальные компоненты электромагнитного поля могут быть зарегистрированы путем расстановки дополнительной второй системы источников поля с количеством не менее двух, которую располагают перпендикулярно первой системе источников, при этом при регистрации первых горизонтальных компонент напряженности электромагнитного поля источники располагают симметрично исследуемому участку по разные стороны от него, а направление тока выбирают одинаковым, при регистрации вторых горизонтальных компонент напряженности направление тока по одну сторону от исследуемого участка сохраняют первоначальным, а по другую меняют на противоположное.
При этом горизонтальные компоненты напряженности электромагнитного поля при первом режиме регистрируют, создавая направление тока в источниках поля одинаковым по отношению к исследуемой неоднородности, а при втором режиме направление в источниках поля с одной стороны неоднородности сохраняют первоначальным, а с другой меняют на противоположное.
На фиг.1 представлен состав электромагнитного поля над геоэлектрической неоднородностью (кривая 1) в виде совокупности поля вмещающей среды (кривая 2) и добавки, связанной с присутствием неоднородности (кривая 3); на фиг.2 - структура электромагнитного поля при размещении одного источника, при размещении двух источников симметрично исследуемому участку при пропускании тока одинакового и противоположного направлений; на фиг.3 структура электромагнитного поля при размещении двух взаимно перпендикулярных систем источников поля при пропускании тока одинакового и противоположного направлений; на фиг. 4 размещение системы источников поля по профилю, проходящему через неоднородность (d расстояние между профилями наблюдений, Dy расстояние между пунктами наблюдений); на фиг.5 антисимметричные относительно положения неоднородности кривые профилирования электрических составляющих аномального поля для нескольких моментов времени: t1 (кривая 4) < t2 (кривая 5) < t3 (кривая 6) < t4 (кривая 7), характеризующих переход от начальных стадий становления поля к поздним; на фиг.6 симметричные относительно положения неоднородности кривые магнитных составляющих аномального поля для нескольких моментов времени: t1 (кривая 8) < t2 (кривая 9) < t3 (кривая 10) < t4 (кривая 11), характеризующих переход от начальных стадий становления поля к поздним; на фиг.7 размещение системы источников поля по двум взаимно перпендикулярным профилям, проходящим через неоднородность; на фиг.8 - результаты применения предлагаемой методики на начальной и завершающей стадиях обработки данных наблюдений.
Представим наблюдаемые компоненты напряженностей электромагнитного поля в виде суммы величин, соответствующих вмещающей среде и аномальных составляющих связанных с влиянием геоэлектрических неоднородностей. Если для величины источниками поля являются сторонние токи , то для величин их роль играют "избыточные" токи в пределах неоднородностей создаваемые вмещающей средой, где
σн электропроводность неоднородности,
σo электропроводность вмещающей среды.
Существуют интегральные соотношения, связывающие наблюдаемые величины напряженностей электромагнитных полей (ЭМП) с полями в глубине разреза при постоянных значениях его параметров.
Таким образом, ЭМП, возбуждаемое в среде с геоэлектрической неоднородностью объема Vн, находящейся на глубине h и отличающейся от среды на проводимость Δε = εн-εo, где напряженность электрического поля, - напряженность магнитного поля. Вектор представляется в виде суммы поля вмещающей горизонтально-слоистой среды и аномальных полей, порожденных вторичными источниками, созданными токами и зарядами в пределах неоднородности. Аномальные ЭМП разделяются по их воздействию на ЭМП гальванического происхождения и ЭМП магнитного происхождения , т.е.
На фиг.1 представлен состав электромагнитного поля над геоэлектрической неоднородностью (кривая 1) в виде совокупности поля вмещающей среды (кривая 2) и добавки, связанной с присутствием неоднородности (кривая 3).
Когда ЭМП в пределах объекта (неоднородности) мало меняется, то практически величина представляется через поле от эквивалентного электрического диполя, а от эквивалентного магнитного поля.
На начальных стадиях становления поля влияние неоднородности проявляется достаточно слабо по типу "однократных отражений", которые впервые достигают точки наблюдения через интервал времени Δt после выключения тока в источнике поля по закону Хевисайда. Величина интервала Δt определяется формулой
где ΔR1 расстояние от точки истока поля до неоднородности;
ΔR2 расстояние от неоднородности до точки наблюдения;
среднее значение сопротивления вмещающей среды;
K коэффициент, включающий магнитную проницаемость mo.
На средних стадиях происходит интенсивное взаимодействие полей от неоднородности и вмещающей среды, требующее описание их через полные выражения, содержащие различные "мультипольные моменты". На поздних стадиях влияние неоднородности на полное поле заметно ослабевает и носит достаточно простой характер.
Такой характер влияния поля от неоднородности предполагает и соответствующую систему наблюдений и реализации алгоритмов, которая на начальных стадиях позволяет получить годографы однократных отражений, а на промежуточных и поздних добиться полного описания наиболее значимых параметров неоднородности. Используя эти представления и поздние стадии становления, получим достаточно простые выражения для напряженности электромагнитного поля над простейшей локальной неоднородностью (сферической, радиуса "а"), где величины определяются через вектор потенциалы неоднородности ,
где вектор-потенциал электрического типа для аномального поля;
вектор-потенциал магнитного типа для аномального поля.
Вектор-потенциал неоднородности (с учетом применения "метода изображений") выражаются через дипольные моменты неоднородности аномальный электрический и аномальный магнитный .
Обобщая приведенные формулы на случай малой неоднородности произвольного вида (а точнее вида "бруска" длины 2а1, ширины 2а2 и высоты 2a3), приходим к возможности достаточно простых корректировок формул (4, 5), вводя обозначения Qн проекция неоднородности на дневную поверхность, Vн объем неоднородности
Переобозначив параметры перейдем к выражениям общего вида
Каждой неоднородности как бы соответствует аномальный электрический и аномальный магнитный диполи, причем
а "Методом изображений" учитывается влияние дневной поверхности как границы раздела сред. Здесь расстояние от точки наблюдения до неоднородности;
x, y координаты точки наблюдения относительно центра неоднородности;
h глубина залегания неоднородности.
Таким образом, в аномальных электромагнитных полях на поздних стадиях фактически отмечается присутствие двух независимых составляющих, т. е.
где векторные функции зависимости аномальных компонент от координат,
f2 и f4 функции зависимости аномальных компонент от времени.
Если величина U
Зная поведение какой-либо из функций от координат, нетрудно определить эпицентр и глубину залегания неоднородности. Все сходится к необходимости разделения аномальных полей, которое проведем следующим образом.
Если расположить источники ЭМП заземленные линии АВ и А1B1 параллельно и симметрично относительно предполагаемой неоднородности, то можно регулировать характер воздействия ЭМП на среду.
Пропуская через А1B1, A2B2, А3B3, А4B4 (фиг.1 и 2) равные по величине и одинаково направленные токи при первом режиме I1=I2 (фиг.1), I1=I2=I3=I4 (фиг. 2), добиваются создания в зоне предполагаемого объекта однородного электрического поля , полностью ликвидируя магнитное поле .
Таким образом, получают первые горизонтальные компоненты (электрические составляющие) электромагнитного поля
Пропуская через А1B1, A2B2, А3B3, А4B4 (фиг.1 и 2) равные по величине, но противоположно направленные токи при втором режиме, т.е. I1=-I2 (фиг.1) и I2= -I1, I3= I1, I4=-I3 (фиг.2), добиваются создания в зоне неоднородности однородного магнитного поля , практически сводя к нулю электрическое .
Таким образом, получают вторые горизонтальные компоненты (магнитные составляющие) электромагнитного поля:
Определив параметры вмещающей среды по измерениям за пределами зоны неоднородности, рассчитывают , а затем, используя формулы (12, 13) по величинам UI и UII, определяют раздельно U
Функции f2(t) и f4(t) отражают влияние на характер аномального поля полей вмещающей среды, а также описывают поведение во времени эквивалентных неоднородности дипольных моментов. Функции и содержат информацию о глубине залегания и эпицентре неоднородности.
Регистрация процессов становления поля в присутствии локальной геоэлектрической неоднородности ведется в диапазоне [tmin; tmax] где величина tmin определяется по формуле (1), а tmax≈ μSH, где S продольная проводимость разреза. Если неоднородность находится на глубине Н=5000 м, то при среднем сопротивлении разреза ее влияние наиболее значительно, начиная с tmin ≈ 1,2 с до tmax ≈ 15 c. Для более высокоомного разреза при этот диапазон сокращается до [0,24 c; 3 c] и т.д. По данным зондирований строится и уточняется геоэлектрическая модель исследуемой площади. Исследования ведутся в пределах площади Q ≈ Qн+h2/2, где h глубина залегания неоднородности, Qн ее проекция на поверхность наблюдения.
На начальных стадиях прогнозируется положение неоднородности по смещению tmin начала ее влияния с изменением расстояния до нее как от источников электромагнитного поля, так и от точек приема. Если задаться изменением tmin на определенный процент α, гарантирующий получение информации
то для минимального шага ΔLmin перемещения источника поля применяют оценочную формулу
где tmin минимальное время становления поля;
минимальное время становления при перемещении источника поля на расстояние ΔL;
расстояние от точки истока поля до неоднородности при перемещении его на расстояние ΔL.
Таким образом,
что при α ≈ 0,09 составит 0,3h. Максимальный шаг перемещения определяется соотношением , т.е. величина ΔL ограничена пределами , а количество позиций выбирается четным (К=2n).
Для площадной съемки возможны два варианта размещения точек наблюдения:
по равномерной сетке,
профильно-площадная.
При этом для оценок шага между точками наблюдения Δy важным критерием является особенность поведения ΔUa/Δy~1/hn в зонах наиболее резкого изменения электромагнитного поля, т.е. Δy~hn•ΔUa, что допускает заметное увеличение Δy с ростом глубины исследования. Поэтому зададимся шагом Δy ≈ 0,05h. Если же фактически выбранный шаг заметно меньше Δy, например при решении задач для нескольких глубин одновременно, то целесообразно группирование данных соседних точек наблюдения.
По особенности поведения различных компонент аномальных составляющих электромагнитного поля процесса становления обосновывается возможность наиболее точного определения параметров аномального объекта. Определим, например, Vz на поздних стадиях при расположении неоднородности посредине между двумя заземленными линиями (фиг.2, 3). При включении тока в линиях в одном направлении (первый режим) будет зарегистрирована величина
U
которая состоит из U
При включении тока в линиях в противоположных направлениях (второй режим) величина U
Подробные параметры геоэлектрического разреза определяют интерактивным путем через минимизацию специального функционала.
Определение начальных значений для итерактивных процедур производится последовательно. Сначала определяются параметры горизонтально слоистой среды в предположении, что разрез заменяется эквивалентными проводящим полупространством через функцию становления поля Fт(ρτ,R,t) и для "ближней волны" наблюдений эквивалентной проводящей пластиной через функцию Fт(Sτ, R,t). Потом, вычитая из наблюдаемых значений Fн(R,t) функцию влияния вмещающей среды Fт(ρτo,t) или Fт(Sτo,t) выделяют аномалийную составляющую и через изучение ее геометрической характеристики определяют эпицентр и глубину залегания неоднородности. При первом режиме изучается характеристика f1z(x,y,z,h), а при втором f3z(x, y, z,h). Четкое разделение ЗМП на симметричную и антисимметричную части позволяют точнее определить параметры неоднородности. Величина r
где y координата точки наблюдения по профилю, проходящему перпендикулярно источникам поля A1B1 и A2B2, отсчитываемая от эпицентра неоднородности. Составив сочетание
получим, что величина r2 определяется по формуле
из которой получим
для устойчивого определения величины h с применением формулы требуется данные расчетов h в различных координатах точках подвергнуть статистическому осреднению.
Независимо можно определить h с применением ΔUza/Δx используя выражения, полученные при двух режимах возбуждения тока:
Поскольку величина U
Во втором варианте такое же отношение можно получить через U
Сходные соотношения получают из формул для ЭДС с заменой f2(t) и f4(t) через пpоизводные от них. Повышение достоверности этого параметра также производится через статическое осреднение данных по ряду координатных точек.
Определение Δσ (для малоамплитудных объектов -DS) произведем, пользуясь выражением функции fi(x,y,h) и данными о вышеприведенных параметрах h и a по любой из измеренных компонент. Взяв в качестве исходных данных величины ЭДС получим два независимых выражения
которые также подвергаются статистическому осреднению.
Получив интегральные параметры , проведем дальнейшее уточнение строения неоднородности, разбивая ее на более мелкие элементы и используя данные измерений при первом режиме возбуждения тока, так как из (4 7) следует, что с уменьшением размеров неоднородности величины аномальных полей одновременно уменьшаются, но электрическое аномальное поле от малой неоднородности становится заметно больше магнитного и для него справедливо правило суммарное поле от мелких объектов равно полю от совокупности объекта.
Для получения уточненных значений параметров ищется минимум функционала Φ
где под Fн понимаются наблюдаемые значения какой-либо из компонент ЭМП, Fт значения теоретической кривой для той же компоненты при заданных традиционных параметрах вмещающей среды {ρi,Hi} и геоэлектрической неоднородности , где hj= h+Δhj, , Δσj= Δσ1+Δσ2j
при условии выполнения уравнений экстремумов
Усложняя модель среды при высокой точности полученных данных, можно дополнительно определить другие параметры геоэлектрического разреза. В частности, по величине вклада поля аномального геоэлектрического диполя, порождаемого неоднородностью, определяются параметры поляризуемости неоднородности за счет замены в формуле (4) для произведения через интеграл
и определения дополнительных экстремумов, связанных с поляризацией неоднородности
где ηн поляризуемость, а βн= 1/τн скорость затухания поляризационных процессов.
Процедура уточнения параметров происходит в несколько итераций путем задания закона изменения параметров. Обозначим Gq= ∂Φ/∂q и зададим последовательность шагов
Δ = (0,1÷0,05)•q(i),
q1= q(i),q2= q1+Δ
где g(i) значение искомого параметра на i-м шаге итерации. Остальные параметры при этом считаются фиксированными. В качестве начальных значений для итераций используются те, которые получены по асимптотическим формулам.
В качестве Fт нами используется выражение
где М число элементов, на которое условно разбита изучаемая неоднородность.
Способ осуществляется следующим образом.
По результатам предварительных исследований выявляют параметры неоднородности (протяженность, глубину залегания, ширину). Рассчитывают площадь этой неоднородности Qн при ее проекции на дневную поверхность. Например, для выявленной неоднородности протяженностью 8 км, шириной 2 км и глубиной 5 км величина площади неоднородности составит Qн=16 км2.
Рассчитывают размер площади, на которой проводят наблюдения , где h глубина залегания неоднородности. В нашем примере Q=28 км2. Рассчитывают максимальное и минимальное расстояния между источниками поля
(при погрешности измерений 10%).
Расстанавливают приемники по системе профилей с шагом между приемниками Δy = 0,05h. В нашем примере Δy составят 250 м, число пунктов наблюдения 112, расстояние между профилями до 1 км.
Располагают не менее двух источников поля (заземленных линий А1B1 и А2B2) параллельно неоднородности по разные стороны от нее, например 4 источника (фиг. 4). Пропускают электрический ток, одинаковый по величине (например, 100 А) и направлению.
Регистрируют в пунктах наблюдения производную вертикальной компоненты вектора магнитной индукции и первые горизонтальные (электрические) компоненты электромагнитного поля U
Учитывая параметры вмещающей среды из зарегистрированных величин выделяют электрические составляющие аномального поля U
На фиг. 5 представлены кривые профилирования электрических составляющих аномального поля для нескольких моментов времени становления. На рисунке видно, что с увеличением времени (t1<t4) аномальное поле сначала меняет знак, что демонстрирует переход от кривой 4 к кривым 5 7, а затем уменьшается по амплитуде, одновременно увеличивая зону влияния на профиле (переход от кривой 5 к кривой 7).
Затем через линии А1B1 и А2B2 пропускают одинаковый по величине электрический ток, но направление его по одну сторону от неоднородности оставляют первоначальным, а по другую меняют на противоположное. Регистрируют в пунктах наблюдения производную вертикальной компоненты вектора магнитной индукции и вторые горизонтальные (магнитные) компоненты электромагнитного поля U
С учетом параметров вмещающей среды из зарегистрированных величин в соответствии с формулой (13) выделяют магнитные составляющие аномального поля U
На фиг.6 представлены кривые магнитных составляющих аномального поля для нескольких моментов времени становления.
Из графиков видно, что с увеличением времени становления (t4>t1) аномальное поле также сначала меняет знак, что демонстрирует переход от кривой 8 к кривым 9 11, а затем уменьшается по амплитуде, одновременно увеличивая зону влияния на профиле (переход от кривой 9 к кривой 11).
Параметры неоднородности , рассчитывают с использованием формул (17 24), а методом минимизации, используя соотношения (25 30) уточняют распределение электропpоводности Dsj в пределах отдельных элементов неоднородности.
При необходимости, дополнительно определяют поляризуемость вмещающей среды и собственно неоднородности.
В случае более сложного строения изучаемой геоэлектрической неоднородности, что характеризуется нарушениями в характере поведения кривых профилирования (отклонениями от симметричного или антисимметричного распределений), выполняются наблюдения от дополнительной системы источников поля, размещаемой по профилю, перпендикулярному первому (фиг.7). В результате трансформаций данных получаем "куб геоэлектрических параметров".
На фиг.8 показан переход от предварительной геоэлектрической модели исследуемого участка, полученной с применением соответствующих геоэлектрических разрезов и срезов по результатам обработки ∂βz/∂t к уточненной детальной модели неоднородности (уточненному кубу геоэлектрических параметров), в которой выделен ряд элементов, характеризующих ее особенности.
название | год | авторы | номер документа |
---|---|---|---|
Способ геоэлектроразведки | 1979 |
|
SU1075831A1 |
СПОСОБ ГЕОЭЛЕКТРОРАЗВЕДКИ | 2009 |
|
RU2411549C1 |
СПОСОБ ПОИСКА ЗАЛЕЖИ УГЛЕВОДОРОДОВ | 1992 |
|
RU2065615C1 |
СПОСОБ ГЕОЭЛЕКТРОРАЗВЕДКИ | 1989 |
|
SU1679880A1 |
СПОСОБ ГЕОЭЛЕКТРОРАЗВЕДКИ В УСЛОВИЯХ ТЕХНОГЕННОЙ ИНФРАСТРУКТУРЫ | 2013 |
|
RU2528115C1 |
СПОСОБ ПРОГНОЗА ЕМКОСТНЫХ ПАРАМЕТРОВ И ТИПА ФЛЮИДОНАСЫЩЕНИЯ КОЛЛЕКТОРОВ | 2013 |
|
RU2540216C1 |
СПОСОБ ГЕОФИЗИЧЕСКОЙ РАЗВЕДКИ ЗАЛЕЖЕЙ УГЛЕВОДОРОДОВ | 2013 |
|
RU2527322C1 |
СПОСОБ МАГНИТНО-РЕЗОНАНСНОЙ ГЕОЛОГИЧЕСКОЙ РАЗВЕДКИ ЗАЛЕЖЕЙ ВОДЫ И УГЛЕВОДОРОДОВ | 1995 |
|
RU2088955C1 |
СПОСОБ МОДЕЛИРОВАНИЯ В ИНДУКТИВНОЙ ГЕОЭЛЕКТРОРАЗВЕДКЕ | 1987 |
|
SU1540513A1 |
СПОСОБ ГЕОЭЛЕКТРОРАЗВЕДКИ | 2005 |
|
RU2291472C1 |
Использование: при глубинных зондированиях для поиска и разведки месторождений нефти и газа. Сущность изобретения: По результатам предварительных исследований выявляют параметры неоднородности. Рассчитывают площадь этой неоднородности Qн при ее проекции на дневную поверхность. Рассчитывают размер площади, на которой проводят наблюдения Q=Qн+h2/2, где h - глубина залегания неоднородности. Рассчитывают максимальное и минимальное расстояния между источниками поля - погрешность измерения. Расстанавливают приемники по системе профилей с шагом Δy между приемниками. Располагают не менее двух источников поля параллельно неоднородности по разные стороны от нее. Пропускают электрический ток, одинаковый по величине и направлению. Регистрируют в пунктах наблюдения производную вертикальной компоненты вектора магнитной индукции и первые горизонтальные компоненты электромагнитного поля U
Δy = 0,05h,
N 2n, где n 1, 2, 4, 6, 8, 10, 12, 14.
где Q размер площади расстановки приемников;
Δy величина шага между приемниками;
ΔLmax максимальное расстояние между источниками;
Lmin минимальное расстояние между источниками;
N количество источников;
Qн размер площади исследуемого участка;
h глубина залегания исследуемого участка;
α погрешность измерений.
СПОСОБ ГЕОЭЛЕКТРОРАЗВЕДКИ | 0 |
|
SU324602A1 |
Печь для непрерывного получения сернистого натрия | 1921 |
|
SU1A1 |
Сб | |||
Зондирование становлением поля в ближней зоне.- М.: Недра, 1976, с | |||
Кипятильник для воды | 1921 |
|
SU5A1 |
Способ геоэлектроразведки | 1979 |
|
SU1075831A1 |
Печь для непрерывного получения сернистого натрия | 1921 |
|
SU1A1 |
Авторы
Даты
1997-03-27—Публикация
1994-05-26—Подача